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ABSTRACT 


MANLEY, MARK B. Experimental Study of Airfoil Trailing-Edge Noise: 
Instrumentation, Methodology and Initial Results. (Under the direc- 
tion of DR. THOMAS H. HODGSON.) 

The mechanism of aerodynamic noise generation at the trailing 
edge of an airfoil is investigated. Instrumentation was designed, 
tested and built and digital signal analysis techniques applied to 
gain insight into the relationship between the dynamic pressure 
close to the trailing edge and the sound in the acoustic far-field. 
Attempts are made to verify some trailing-edge noise generation 
characteristics as theoretically predicted by several contemporary 
acousticians. 

A NACA 0012 airfoil of 2.0 foot chord and 1.5 foot span was 
mounted in the laminar flow of the potential core of an open jet. 

The test was conducted ill an anechoic environment thereby providing a 
quiet-flow test capability. Two-dimensional boundary layer flow was 
established over the two surfaces of the airfoil. Means were used to 
"trip" each boundary layer into turbulence and thereby maintain 
trailing edge flow conditions such as are found on full-scale aircraft 
wings. 

Considerations are discussed which had bearing on the design of 
a miniature semiconductor strain-gauge pressure transducer and 
associated electronic amplifier circuitry. Several of the transducers 
were tested and then mounted in the "upper" and "lower" surfaces of 
the airfoil in the vicinity of the trailing edge. Signals from these 
near-field transducers and other far-field condenser microphones were 
analyzed using digital signal processing techniques. The origins of 



the computer analysis functions are briefly discussed and some errors 

encountered in their use are covered. Instrumentation of the overall 

experiment is described and the results and conclusions stemming 

from the computer analysis are presented. 

It is found that the noise detected in the far-field is comprised 

of the sum of many uncorrelateo emissions radiating from the vicinity 

of the trailing edge. These emissions appear to be the result 

acoustic energy radiation which has been converted by the trailing- 

edge noise mechanism from the dynamic fluid energy of independent 

strearawise "strips" of the turbulent boundary layer flow. It is also 

found that the far-field sound power level scales on the freestream 

2 5.3 

velocity "u" raised to a power of 5.3; i.e., p oc u * . This is in 

5.0 

good agreement with the theoretically predicted u 
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1. INTRODUCTION 


The understanding of the aerodynamic noise generation phenomena 
has been a long sought goal of acousticians. Up until the last few 
years, however, most engineering attention and research has been 
focused on the roost obvious source of aircraft noise, namely the 
engine. Recently, airframe component noise due to struts, flaps, 
wheel cavities, etc. has been studied because of its possible impor- 
tance on the aircraft landing configuration. 

Theoretical studies by many acousticians including Lighthill 
[30, 31], Curie [15], Powell [33], Ffowcs-Williaras and Hall [20], 
Crighton [14], and most recently Howe [27] have indicated that purely 
aerodynamically induced noise may be generated by the interaction of 
turbulent air flows with rigid planar surfaces. That is to say that 
the surfaces themselves need not vibrate for the noise to be radiated. 
In particular, enhancement of the noise is predicted to occur at the 
trailing edges of the planar surfaces as turbulent boundary layers 
develop as a result of forward flight. (See Howe’s review paper [27].) 

For the case of turbulent flow over an infinitely large rigid 
flat plate, theory predicts that the pressure (force) fluctuations 

of the turbulent boundary layer flow generate dipole noise sources. 

2 

The sound power level or p level of dipole sources generally 

depends on the sixth power of a characteristic flow velocity "u," 

2 6 

typically the freestream flow velocity u ; i.e., p *= ku , where k is 

a constant. However, because these dipole sources are reflected in 

the rigid surface of the plate, they give rise to quadrupole noise 

emissions as sensed in the acoustic far field. Typically, quadrupole 

2 8 

noise scales on the eighth power of the velocity; i.e., p = ku . 


2 

For the case of a semi-infinite rigid flat plate, the theories 

of both Ffowcs-Williams and Hall [20] and Howe [27] predict that the 

large extent of the plate can possibly cause enhancement of the sound 

field due to scattering at the edge. This enhancement of the 

quadrupole noise has been theorized to scale on the fifth power of 

2 5 

the velocity; i.e., p = ku . This makes the efficiency of the 
trailing-edge noise phenomena lie between that of the highly efficient 
monopole (p^ = ku 4 ) and the aforementioned dipole (p^ = ku^). Thus, 
the extensive rigid planar wing surfaces of aircraft such as a Boeing 
747 jumbo jet may actually become the dominant contributors to the 
overall radiated noise sensed on the ground as the large aircraft 
approach airport runways at low altitudes with the engines throttled 
back. 

This purely aerodynamic nonpropulsive related noise source has 
been recognized as a possible lower bound or "barrier" in preventing 
the noise floor of future aircraft from being reduced. H. G. Morgan, 
Division Chief of the NASA Langley Acoustics and Noise Reduction 
Division has stated: 

We have made enough progress in quieting jet engines 
that we are nearing a point where the airflow over the 
airframe becomes what is actually heard on the ground. 

. . . We cannot quiet engines more than another 5-8 dB 
before we run into the new noise floor. [13] 

Figure 1.1 illustrates that as engines for a given aircraft 
size are becoming quieter, increasingly larger aircraft and cor- 
respondingly larger planar surfaces areas are radiating ever increas- 
ing amounts of aerodynamic noise. 

As part of its Federal Aviation Regulations, in 1974 the Federal 
Aviation Administration issued a noise standard for aircraft, FAR-36, 


Present 


Time & Size 


Figure 1.1 Past and future trends showing (a) decreasing 
engine noise for a given size aircraft due to 
better technology and (b) increasing aerodynamic 
noise due to increasing surface area 

which states in part that it ", . . does not agree with those that 
believe the use of specific noise abatement/flight procedures . . . 
would, by itself, provide the necessary noise control and abatement" 
[18]. It is the opinion of the FAA that noise control must be under- 
taken at the source of the noise generation. 

This has led experimentalists to attempt verification of the 
different facets of the proposed acoustic theories. Experimental 
work had been accomplished by many parties including Siddon [37], 

Fink [21], Healy [26], Putnam et al. [34], Hahn [22], Yu and Tam [40], 
and others (see Hardin's review paper [24]). The most promising 
theory at present appears to be Howe's "Unified Theory of Trailing 
Edge Noise" [27] which is a variation of Lighthill's treatise on the 


subject In 1952 f 30 ] with sweeping modifications to suit recent 
findings and theories. 

However, only lately has extensive experimental work of labora- 
tory quality involving the latest signal processing techniques been 
possible, These can be used to accurately confirm the legitimacy of 
Howe's or other proposed theories. It is to this end that an experi- 
mental study was initiated under Contract No. NSG 1377 for testing 
to be performed at the quiet flow facility located at the NASA 
Langley Acoustics and Noise Reduction Laboratory, Langley, Virginia. 
It is the purpose of this thesis to discuss the preparations for and 
initial findings of the first stage of testing, known as Series I 
tests. 

A NACA 0012 airfoil was used for the tests over which turbulent 
boundary layers could independently develop on both surfaces. 

Boundary layer characteristics and flow speeds were achieved which 
approximated those actually encountered on the surface of full-size 
aircraft wings. Thus, the primary differences between model and 
actual wing are those of geometric scale, i.e., dimensions of chord 
length and trailing-edge span. 

The considerations that were involved in the design of a 
miniature strain-gauge pressure transducer and associated electronic 
amplifier circuitry are discussed. Several of the pressure transduc- 
ers were tested and then mounted in the upper and lower surfaces of 
the airfoil in the vicinity of the trailing edge. Signals from 
these near-field transducers and other far-field microphones were 
analyzed using digital signal processing techniques. The origins 


of the analysis functions are briefly discussed and some errors 
encountered in their use are covered. Instrumentation of the overa 
Series I experiment is described and the results and conclusions 
stemming from the computer analysis are presented. 
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2. TEST CONFIGURATION 
2. I Airfoil 

In order to generate the turbulent boundary layer at a trailing 
edge in a controlled manner, a symmetrical NACA 0012 airfoil was 
employed. It was machined from solid aluminum stock as shown in 
Figure 2.1. Dimensions are given in Figure 2.2 (a). Slots were 
machined to allow small strain gauge pressure transducers to be 
mounted flush with the surface. The trailing edge as manufactured 
had a bluntness of 0.10 inch with corners of 0.030 inch radius 
(see Figure 2.2 (b)) which allowed some of the pressure transducers 
to be mounted in very close proximity to the edge itself. (In later 
tests. Series II, a sharp extension was added.) 

2.2 Facility and Experiment Geometry 

The test was conducted in the recently-modified open-jet quiet 
flow anechoic test facility located at the NASA Langley Acoustics and 
Noise Reduction Laboratory (ANRL) . Figures 2.3 and 2.4 depict the 
basic facility setup and define the coordinates used throughout this 
thesis. Also the positions of the one-half inch Brilel and Kjaer 
type 4133 condenser microphones are shown. 

A rectangular nozzle with dimensions 1.0’ W x 1.5' L was mounted 
on the air supply opening and was capable of airflow velocities in 
excess of 250 feet per second. Essentially laminar flow was achieved 
at exit through the use of the nozzle profile and flow straightners 
and silencers located upstream in the air supply tunnel. The airfoil 
was mounted in the potential core of the jet so that laminar flow was 
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Figure 2.1 Unmounted NACA 0012 airfoil with Kulite strain-gauge 
pressure transducers installed 
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Figure 2.2 Dimensions of the unmodified NACA 0012 airfoil (general 
or basic configuration): (a) overall, showing boundary 

layer trip at 15% chord, (b) detail of trailing edge 





Figure 2.3 Overall view of the experimental setup and facility 
as viewed from the control room door ("upper" 
surface of the airfoil) 
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maintained in the freestream up to and beyond the trailing edge of 
the airfoil. 

2 . 3 Flow Conditions 

Two-dimensional flow was maintained over the surface of the 
airfoil by means of sideplates mounted on either side of the nozzle 
opening. The sideplates had foam attached to the rear edges outside 
of the main airflow to prevent detrimental edge- tone and diffraction 
effects from occurring (Figure 2.5). 

To obtain fully-developed turbulent boundary-layer flow within 
such a short distance as the 24-inch chord length of the airfoil, a 
boundary layer "trip” was used. It consisted of a strip of coarse 
grit 1.0 inches wide centered at 15 percent chord as measured from 
the leading edge. Both upper and lower surfaces were provided with 
such a trip to facilitate the growth of statistically equivalent but 
separate turbulent boundary layers. The turbulent boundary layer 
characteristics tten closely represented the flow characteristics 
over a full-scale wing. 

The airfoil was tested at angles of attack equal to a = 0°, 

+5 , and ±10°, but little dependence on a was detected so only 
a = 0° results are presented. 

Freestream flow velocities ranged from 50 ft/sec to 250 ft/sec. 
Most testing was done at a velocity of 225 ft/sec which is equivalent 
to 900 rpra of the air supply fan. Most results in this thesis relate 
to the preliminary study which was conducted using a flow velocity of 
225 ft/sec. 
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Figure 2.5 Airfoil mounted at a = 0° angle of attack with 
Kulites, sideplates, foam, nozzle opening and 
boundary-layer trip visible 
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Boundary layer and airflow conditions were established using 
widely accepted standardized techniques. Some data parameters per- 
taining to the boundary layer conditions when the freestream flow 
velocity was 225 ft/sec are listed in Table 2.1. 

Table 2.1 Initial test boundary-layer characteristics 


Parameter 

Value 

Description 

U 

oo 

225 ft/sec 

freestream velocity 

U 1 

2 1 A ft/sec 

velocity at edge of B.L. 

U 

T 

7.01 ft/sec 

friction velocity 

6 

1.24 in. 

boundary layer thickness 

5 i 

0.154 in. 

displacement thickness - 6/8 

0 

0.162 in. 

momentum thickness 

R e 

1.67 x 10 4 

Reynolds it based on 0 

R 

c 

2.88 x 10 6 

Reynolds it based on chord 

V 

1.56 x 10 4 ft^/sec 

kinematic viscosity 
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3. PRESSURE TRANSDUCER DESIGN 

To sense the pressure fluctuations occurring in the turbulent 
boundary layer at the trailing edge of the airfoil, a very small 
pressure transducer or microphone was designed and flush-mounted as 
near the trailing edge as possible. The considerations that went 
into the design, calibration and installation of this transducer 
follow. 

3. 1 Design Considerations 

Several competing phenomena mandated that the transducer dimen- 
sions be carefully chosen to provide the best compromise design. Of 
prime importance was that the transducer be small enough so that 
several could be clustered side-by-side on the upper and lower 
surfaces at the trailing edge. Kulite Semiconductor Products, Inc., 
of Ridgefield, New Jersey, was contacted due to their prior experience 
in manufacturing miniature transducer for this purpose. Manufacturing 
and physical limitations produced a transducer of overall outside 
dimensions of 0.375” L x 0.125” W x 0.045” D. This permitted a 
minimum center-to-center transducer diaphragm spacing of 0.125 
inches (Figure 3.1). 

It was desired that the flow of air over the airfoil surface be 
disturbed as little as possible by the transducer so that only valid 
boundary-layer pressure-field data would be collected. Therefore, 
the transducer design incorporated a pinhole opening in the upper 
surface which led to an internal cavity. The bottom of the cavity was 
formed by a diaphragm of semiconductor material on which a strain- 
guage sensing element was etched. 
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It was predicted that increasing the diameter of the pinhole 
opening would have as a result an increasingly pronounced effect on 
the frequency spectrum of the naturally occurring turbulent flow at 
the airfoil surface. This is because the hole diameter approaches 
the order of dimension of the boundary layer thickness [10]. However, 
decreasing the hole size lowers the natural frequency of the Helmholtz 
resonator formed by the cavity and pinhole opening. Then the frequen- 
cies thought to exist in the extent of the pressure-spectrum become 
prone to resonant magnification within the transducer cavity, again 
causing erroneous data readings. 


3. 2 Parameter Calculations 

Using normalized wind tunnel wall-pressure spectrum data and 
methods outlined by Bull and Thomas [10] the approximate upper bound 
of frequencies that could be measured accurately was computed. They 
conclude that measurements made with a pinhole microphone are prone 
to error for frequencies such that: 


> 0.01 

u u 

T 00 

-4 2 

where: v = 1.56 10 ft /sec, kinematic viscosity of air; 

= 7.01 ft/sec, friction velocity; 
u^ = 225 ft/sec, freestream velocity; 
u> = 2vf rad/sec, radian frequency. 


(3.1) 


The resulting value for f = 16.1 kHz. 

Inspection of curves presented in reference [10] suggested that 
for pinhole transducers, the term 
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du 

= 50 (3.2) 

v 

where d is the pinhole diameter. Substituting previous values for u t 
and v yielded a diameter d = 0.0135 inches. 

It was necessary to place the Helmholtz frequency well above the 
highest frequency to be analyzed, i.e., well above 16 kHz. Thus, 
accompanying phase and amplitude effects on the pressure spectrum were 
kept to a minimum. Keeping in mind the dimensions of the basic 
transducer model to be modified and the range of permissible pinhole 
diameters, the remaining critical dimensions for a high Helmholtz 
frequency were obtained. Repeated calculations with the following 
equations were made in the selection process. 

The basic Helmholtz resonant frequency equation employed [29] is 
given as: 



S 

( £+1 . 7a) V 


Hz 


(3.3) 


where: c = 13,500 in. /sec, speed of sound; 

2 2 

S = Tra in. , pinhole cross-sectional area; 

a = radius of the pinhole, in.; 

Z - length of the pinhole opening, in.; 

2 3 

V = Trr h in. , cavity volume; 
h - cavity height, in.; 
r = cavity or diaphragm radius, in. 

The '‘sharpness” of the resonance or quality factor Q is given 


by [29]: 
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q - 2, [ < W . ^ . ) 3 V 1 * 

S 


(3.4) 


The upper and lower "half-power" or -3dB frequencies for the 
resonant peak may be found by [7]: 


f 2 * f o + -S)>- »*• 


(3.5) 


f 1 " f o ^ " 2Q^ ’ Hz ’ 


(3.6) 


where: £. < f < f„ , Hz. 

1 o 2 

By simple algebraic manipulation the -3dB bandwidth results: 



(3.7) 


The amplification and phase effect due to the resonance is found 
using a transfer function quantity [9]: 


H(f) 




o 



} 


(3.8) 


where H(f) is a complex function of frequency. 

The finalized dimensions of the finished transducer, henceforth 
referred to as a "Kulite," are given in Figures 3. 1-3. 3. The 
accompanying manufacturer's specifications and computed values result- 
ing from the above equations are given in Table 3.1. 


TR ; 
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Figure 3.3 Internal Kulite dimensions and construction details 


Table 3.1 Specified and calculated Kulite performance parameters 


Parameter 

Value 

Description 

f . 

o,d 

100 

diaphragm resonant frequency 


25 mV/psi 

nominal sensitivity 


15 VDC 

excitation voltage 

Z 

o 

750 

electrical output impedance 


±l%/±2 psi 

combined nonlinearity & hystersis 
in percent over the given pressure 
range about ambient conditions 

f 

o 

63.2 kHz 

Helmholtz resonant frequency 

Q 

24.6 

quality factor 

Af -3dB 

2.57 kHz 

-3dB bandwidth 

f i 

61.9 kHz 

lower -3dB frequency 

f 2 

64.5 kHz 

upper -3dB frequency 
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Various calculated values for resonant magnification factor and 
phase are presented in Table 3.2 and plotted in Figure 3.4 as func- 
tions of frequency. 


Table 3.2 Absolute value and phase of the resonant magnification 
factor for selected frequencies 



Frequency, kHz 


<p(f), degrees 


10.0 


1.03 

0.378 


16.0 


1.07 

0.630 


61.9, 

< f i> 

17.7 

44.7 

- 

63.2, 

< f o> 

24.6 

90.0 


64.5, 

<f 2 > 

17.1 

134.7 


100.0 


0.663 

{77.6 


Standing wave effects in a cavity of such small dimensions were 
not considered to be at all significant. 

3.3 Pressure Coupler Design and Use 
The Kulites were to be mounted symmetrically on the "upper" and 
"lower" surfaces of the airfoil. To best match the phase, frequency 
response and sensitivity of symmetric Kulites to one another, an 
acoustic calibration pressure coupler was employed (Figure 3.5). 

It was constructed of steel and had cavity dimensions as shown in 
Figure 3.6. 

The coiled tube connected to the coupler was provided to relieve 
static pressure build-up as the microphones were inserted. It was cut 
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Figure 3.4 Calculated (a) magnitude and (b) phase of the magnifi- 
cation factor resulting from the Helmholtz resonance of 
the Kulite pressure transducer cavity 






Figure 3.5 Pressure coupler components showing (a) driver, Kulite 
and reference microphones with cavity visible, (b) unit 
in operation 
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long and had a small inner diameter to place any pipe resonances well 
below the lowest measurement frequency. 

The principle of operation for a pressure coupler/calibrator is 
straightforward. A high quality microphone with known frequency and 
phase response characteristics is used as a reference. It is inserted 
into one side of the coupler until it reaches the internal cavity. 
Another microphone which may be of lesser quality is used as a driver 
to provide a noise source. It too is inserted into the coupler until 
it reaches the cavity and is facing the reference microphone. The 
pinhole opening of the Kulite is also provided access to the cavity. 
Except for the static pressure relief tube, the cavity is now sealed. 
The very close proximity of driver, reference and Kulite microphones 
together with the small volume of the cavity insure that the acoustic 
pressure is the same throughout when the driver emits a tone. Thus, 
the dynamic pressure response of the Kulite may be accurately cali- 
brated against that of the reference microphone. 

The small dimensions involved prevented standing wave resonances 
from becoming a problem in the frequency range of interest which was 
< 100 kHz. 

A Briiel and Kjaer type 4134 one-half inch condenser microphone 
was used as the driver. It was connected to a high impedance 200 VDC 
power supply to provide the polarization voltage necessary for opera- 
tion as a transducer. A swept-sine signal generator provided the 
source signal. A Briiel and Kjaer type 4138 one-eighth inch condenser 
microphone was used as the reference. A Briiel and Kjaer type 2619 
preamplifier and type 2607 measurement amplifier connected the 
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reference microphone to a type 2971 phase meter. The output of the 
phase meter was connected to an X-Y plotter. Synchronizing the X-Y 
plotter to the signal generator allowed plots of phase to be obtained 
from 1 kHz to about 30 kHz. Figure 3.7 shows the actual measured 
values of the magnitude and phase lag of the Kulite pressure response 
as functions of frequency. 

The phase and magnitude of the frequencies near Helmholtz 

resonance were checked utilizing the coupler, signal generator, micro 

phone amplifier and a Tektronix model 3103N oscilloscope. Only one 

Kulite was measured in order to verify the accuracy of the predicted 

resonance. The measured f =62 kHz was within i 1 kHz of the pre- 

o 

dieted 63.2 kHz and the phase was within an acceptable margin 
(if < 1.0°) as observed using Lissajous figures. 

It was found during calibrations that the strain-gauge sensing 
element in the Kulites required a considerable warm-up period before 
stable operation was achieved. This period was avoided during the 
weeks of measurement by leaving the strain-gauge excitation power 
supply on at all times. 


3.4 Kulite Mounting Positions 

The Kulites were mounted flush to the surface with epoxy. The 
wiring was routed through channels machined into the airfoil surface. 
The geometric spacing of the Kulites was based on prior knowledge of 
the approximate rate of falloff of the correlation of the pressure 
field [22]. The possible failure of the transducers was also con- 
sidered so some of the spacings between pinhole centers were repeated 
(Figure 3.8). This redundancy also allowed a second measurement of a 
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Figure 3.7 


Measured values of (a) magnitude and (b) phase lag of 
the Kulite pressure response as functions of frequency 




Figure 3.8 Detail of the upper airfoil surface and trailing 
edge, showing mounted Kulites and epoxy-filled 
grooves which accommodate wiring 



mm 



32 

parameter to be made under similar dimensional conditions for 
verification purposes. Figures 3.9 and 3.10 show the positions of 
the Kulites as used in this initial study, 

3.5 Kulite Amplifier Design 

A circuit was needed to amplify the double-ended electrical 
signal originating at the corners of the semiconductor strain-gauge 
bridge network. Desirable characteristics included: 

1) Low electrical noise of operation at high gain factors. 

2) High common mode rejection ratio (CMRR) > 100 dB to exclude 
60 Hz hum and radio frequency noise pickup. 

4 

3) Suitable gain on the order of 10 V/V to boost signals 
from a level of 25 pVrms to a voltage level compatible with 
most analysis instrument inputs, about 0.5 Vrras. 

4) Wide bandwidth capability to prevent amplitude or phase 
distortions. 

5) dc coupling at the input to allow the dc excited bridge to 
operate properly. 

6) ac coupling at the output to prevent small dc voltage 
offsets due to the transducer from being amplified ten- 
thousand fold and clipping or overloading the output 
circuitry. These offsets may originate from barometric 
and/or temperature shifts acting on the bridge. 

These criteria were easily and economically met by using two 
stages of amplification. Thus, the goals were accomplished: 

1) The gain-bandwidth product for any one amplifier stage 
being fixed, wide bandwidths for each stage were achieved 
by requiring only relatively low gains. This kept ampli- 
tude and phase distortions to a minimum. 
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2) The first stage was designed with a relatively low gain of 
50 V/V (as compared with the overall gain of approximately 
5000 V/V) . This prevented small input dc offsets from 
overloading the output. The balanced input of the first 
stage operational amplifier used provided a minimum CMRR 
of 110 dB. 

3) The two-stage design allowed the ac coupling to be placed 
between stages rather than at the output of the amplifier. 

A rather small-valued capacitor provided the coupling to 
the high but constant input impedance of the second stage. 
The low frequency cutoff achieved was ; 100 Hz. 

A prototype circuit of the above general description was built 
at North Carolina State University using two operational amplifiers. 
The circuitry as actually used at NASA was constructed by Wyle 
Labs, Hampton, Virginia, and is shown schematically in Figure 3.11. 
Neff model 122 amplifiers operating at unity gain were connected to 
the outputs to provide line driving from the anechoic chamber area 
where the Kulite amplifiers were situated to the control room where 
all other amplifier circuitry was located. 

Ithaco model 455 active filters operating at voltage gains of 
25 to 45 dB and set for low frequency cutoffs of 100 Hz provided the 
remaining voltage boost required by this design. Flat frequency and 
linear phase response beyond the 100 kHz range of testing was 
achieved with this configuration. 


Figure 3.11 Schematic of the power supply and one of eight 
first-stage amplifiers for the Kulite pressure 
transducers 













4. SIGNAL ANALYSIS FUNCTIONS 


It is the purpose of this section to provide a brief review of 
the mathematical development and usefulness of the functions used 
for signal analysis in this study. 

4 . 1 Auto-Correlation Function 

Correlation as it applies to signal processing is generally under- 
stood to be a means of providing an indication of the amount and 
nature of the time dependency of one time signal x(t) on another y(t). 

For the case where the two time signals are identical, the "auto- 
correlation" results. Its calculation provides an indication of the 
influence of "past" events on a chosen "present" event time t as a 
function of the time delay or lag t between the two events. It is 
defined in its infinite-limit continuous form as: 

I +T / 2 

R (t) = limit — / x(t) x(t+x)dt (4.1) 

304 T -+oo 1 -T/2 

where T is the length in time over which the multiplication is carried 
out. Because the signal is correlated with itself, the resulting 
function is symmetric about t = ti. The maximum value for the auto- 
correlation will always occur at t = 0 which also corresponds to the 
mean-square value of the signal. 

The auto-correlation is useful for characterizing a single 
signal as to its periodic versus nonperiodic content, which includes 
echoes, periodic pseudo-random noise, etc. [6, 35]. Examples of some 
characteristic auto-correlation curves may be found in the literature 
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Physical limitations of data storage space and computational time 
prohibit infinitely long data collection; thus, finite time dependent 
functions are employed to provide approximate values in lieu of the 
exact values provided by the infinite continuous functions. 

The auto-correlation may be re-defined for continuous but finite 
data record lengths as: 

1 +T / 2 

R v (t) = - / x(t) x(t+i)dt (4.2) 

XX J _ T/2 

where the ,, '* M indicates the result of many averages (Section 7.1). 

If, however, the auto-correlation is to be computed on a digital 
computer, the continuous time signal must be sampled and quantized 
into discrete data points (see Section 6.1). Then the straightforward 
computation involves the equation: 

(N-n) 

ft xx,n (nh> = J, x(ih) x[(i+n)h] (4.3) 

1=1 

where: n = 0, 1, 2, ...,r . 

N is the total number of data points in the sample record, nh is the 
incremental time displacement, r = T max ^ is t * ie max i mum lag number, 
and h is the time interval between consecutive sample points. 

In many instances information in the frequency domain is more 
desirable for characterizing signals since amplitude and phase as 
functions of frequency are obtained. These are not directly avail- 
able when the same data are processed through a time function such as 


correlation. 
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4 . 2 Fourier Transform and Auto-Spectral Density Functions 
The most common time - to - frequency domain transform func- 
tion used is the forward Fourier transform. In its infinite-limit 
continuous form it is "double-sided” in frequency about f = 0: 

■fon 

FT{x(t) } = X(f) = / x ( t) e " 32 nft dt (4.4) 

— Ou 

where: < f < +«•> , 

x(t) is some piecewise-smooth function of time and the "~" denotes 
a complex quantity. At times, rather than presenting a complex 
quantity in terms of a magnitude with phase (which is known as 
"phasor" notation) the quantity is separated into its real and 
imaginary components. The use of Euler’s equation [39]: 

e = cos(z) - jsin(z) (4.5) 

where z is a general variable yields the real and imaginary components 
of the frequency ip»ctrum when applied to Equation (4.4): 

X(f) «= |x(f )J[cos(2Ttf t) - jsin(2Tift)] 

- |x(f)| cos(2irft) - j|x(f) |sin(2irft) . (4.6) 

Fourier transformed functions are often called' "spectral densi- 
ties" because they show how the frequency spectra are amplitude 
weighted . 

The finite Fourier transform for obtaining the spectral density 
function is defined as: 


AO 


X ( f , T) 



x(t)e^ 27rft 


dt 


(4.7) 


where x(t) is a continuous function of time but is integrated over the 
finite period T. The resulting function is no longer a continuous 
function of frequency as in Equation (4.4) but is instead represented 
by discrete frequency bands of bandwidth: 

BW = ~ , Hz (4.8) 


with an infinite number of bands in the positive and negative frequen- 
cy axis yet symmetric about zero. 

Since, as stated previously, only finite capabilities exist for 
collecting and analyzing the data, the function x(t) must be repre- 
sented by a finite time length sample record x(ih) (Section 6.1). 

Thus, the period of integration T as used in Equation (4.7) is replaced 
by the product Nh and the discrete finite Fourier transform is given 
by: 


(n-1) 

X yi (f n ,Nh) = h l x(ih) 


-_12Trni 



(4.9) 


where: n = 0, 1, 2, ..., (N-1) , 

f = Nh * t * ie i ncremen tal frequency, (4.10) 

x(ih) is the finite length time history of events comprised of N data 
points h seconds apart, and N, h, and i are as previously defined. 


A I 

For actual direct computation of the spectral density on a 
computer, Euler's relationship, Equation (A. 5) may he used on (A. 9) 
to obtain [2, 6, 28, 32]: 


(N- 1 ) 

X (f ,Nh) «= h l x(ih)cos[ 2tt(~) ] 

~n n i ^ o N 

(N-D . 

- jh' l x(ih)sin[2TT(-2i)] 
i=o 


(A. 11) 


where the variables are as defined previously. 

A . 3 Fast Fourier Transform 

A much faster method of computing the discrete finite Fourier 
transform known as the fast Fourier transform, or FFT, was developed 
by Cooley and Tukey in 1965 [12] and has since been refined by other 
mathematicians. The FFT algorithm used during this study was devel- 
oped by R. C. Singleton in 1969 [38] and has proven to be an accurate 
and efficient FORTRAN version of the original Cooley-Tukey FFT 
algorithm. 

ft ft the computer implementation of the FFT, the input data are. 
assumed to be complex valued and hence formatted as two time series 
arrays — one the real part of the input data and the other the 
imaginary. A single real-valued time series x(ih) could be inserted 
into the real time series array while all points of the imaginary 
tims-. series array are set to zero. A more efficient method suggested 
by Singleton [38] is to insert odd-numbered data points of the real- 
valued time series into the real array and even-numbered data points 
into the imaginary array. Although a sorting algorithm must be used 


on the output of the FFT algorithm to yield interpretable results, 
storage space is halved and the computation time considerably speeded 
up . 

4.4 Properties of Spectral Density Functions 
Since spectral density functions are two-sided and symmetric 
about f = 0, and because negative frequencies have no physical meaning 
when analyzing real-life time events, the single-sided spectral density 
function is defined as: 

G(f ) =•- 2 X(f), 0 < f < + « . (4.12) 

Due to the squared-quantity nature of correlation functions, the 
Fourier transform of the auto-correlation is a "power-like” quantity. 
Hence, it is termed the "auto-power spectral density" or auto-PSD 
and is normally converted to single-sided form so that: 

G (f) * 2 FT{R (i) } . (4.13) 

XX XX 

Note that because the auto-correlation is a real and symmetric 
function of time, the resulting auto-spectral density function is 
real and symmetric [6, 45], Thus, the associated auto-PSD is real as 
indicated in Equation (4.13). This is a general result of the 
integral of the sine function with infinite limits. Equations (4.4) 
and (4.5), being equal to zero when real and even functions are 
transformed. 

4.5 Computer Implementation of Signal Analysis Functions 
This same result is obtained digitally on the computer by multi- 
plying the complex FFT result with its complex conjugate as denoted 
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by the *: 

C (f , Nh) =-• 2 X (f , Nil) X * (f , Nh) (A. 14) 

xx n ~n n --n n* 

- (A + jB) (A - jB) 

2 2 

- A + B , real , 

where A and B represent the respective real and imaginary frequency 
components of the complex FFT result. This yields the requisite 
power-like squared quantity that is obtained from the theoretical 
discrete finite Fourier transform of the auto-correlation. 

Equation (A. 13) indicates that the auto-correlation and the auto- 
PSD constitute a Fourier transform pair. This is known as the 
Weiner-Khintchine relationship [2, 32, 35] and implies that the 
inverse Fourier transform of the auto-PSD yields the auto- 
correlation: 


+00 


R 


xx 


<0=4:/ <f> 


j2irf t 


2tt ^ 'xx 

— 00 


df 


or in discrete finite form [2, 6]: 


(N-n) 


j2nni 

. »■> »*/ u 

W" h > = 2^ hG ™,i«V me 


(A. 15) 


(A. 16) 


where: n = 0, 1, 2, ..., (n-1) 

and the other variables are as defined previously. 

This is a very useful result as it allows the inverse FFT 
computer algorithm to compute correlation function values from spectral 


■*!* 


density values. It is much quicker and more economical than the 
straightforward point-by-point time domain computations. It also 
allows the correlation to be computed on data that have been 
•‘filtered" into frequency bands of interest. This "filtering" may 
be accomplished by editing the frequency spectrum of the data after 
they have been forward transformed by the FFT. This roundabout 
method of computing the correlation has been called "circular correla- 
tion" and is the method used to calculate the correlation values in 
the program CIRCXCOR listed in Appendix 11.4 [2, 6, 32]. 

4.6 Cross-Correlation and Cross-Correlation Coefficient 

The cross-correlation function is similar to the auto-correlation 
function in that two time series are compared. The cross-correlation 
is more general in that the two time series are of different origin — 
from two transducer locations, for instance. Thus, it is seen that 
the auto-correlation is merely a special case of the cross- 
correlation. 

Some of the more common characteristic cross-correlation curves 
are similar to the auto-correlation curves of references [2, 6, 19] 
except that the axis of symmetry no longer lies at t = 0 but rather 
at some other time delay. Nor is there necessarily an axis of 
symmetry because real-life phenomena and measurement systems may 
contain extraneous noise in data channels, nonlinearities such as 
dispersive media, etc. 

The cross-correlation is useful for detecting the propagation 
time delay between transducer locations, the existence of a signal 
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(not necessarily periodic) that is buried in noise, and the determina- 
tion of transmission paths [6, 35]. 

The cross-correlation of two continuous real signals, x(t) and 
y(t), as a function of the arbitrary delay or lag time t is defined 
as: 

, +T/2 

R (x) = limit ^ J x(t) y(t+r)dt . (4.17) 

Xy 1 -T/2 

As for the case of the auto-correlation function, physical 
limitations necessitate the use of finite time versions of (4.17). 
Therefore, the cross-correlation function is re-defined for continuous 
but finite data record lengths as: 

1 +T/2 

R (t) “ rr J x(t) y(t+T)dt (4.18) 

xy ~t/2 

and in the discrete finite-time form as: 

R _(nh> = FF l x(ih)y [ (i+n)h ] (4.19) 

, 11 ” n 

where: n = 0, 1, 2, .... r 

and the other variables are as defined previously. 

A normalized version of the cross-correlation is often more 
informative as to the "percent" of correlation between two signals. 
This is useful for comparing the results of one test with the 
results of another without regard as to absolute reference values or 
scaling of data. This normalized cross-correlation function is called 
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the "correlation coef f icient" p (i) and its value always will lie 

xy 

between +1 and -1. Provided the two signals x(t) and y(t) have zero 
mean, p (t) may be computed using the auto-correlation functions 

xy 

with t = 0 [2, 6, 32]: 


p xy <T) ‘ R xy <l)/[R xx (0>R yy (0)r ’ - p xy - +I ' (4 ’ 20) 


4 . 7 Cross-Spectral Density and Ordinary Coherence Functions 
In the same vein as for the single-channel auto-PSD the two- 
channel "cross-PSD" is obtained from the cross-correlation. But 
because the cross-correlation is in general not an even function, 
the cross-PSD is usually complex: 

G (f) = 2 FT{R C t ) } . (4.21) 

-xy xy 

In a simplified likeness of Equation (4.14), the cross-PSD as 
computed from the output of the FFT is given as: 


G (f) = 2X(f )Y*(f ) (4.22) 

-xy 

where the quantities are averaged for reasons given in Section 7.1. 

A normalized version of the cross-PSD that is void of any phase 
information is called the "ordinary coherence" or just "coherence." 
Maintaining the notation simplification used above with the additional 
stipulation that all following spectral quantities are assumed to be 
averaged and real or complex as necessary, the coherence is defined as 
a real number between 0 and 1: 


Y xy ( ° " G 


G (f) 
' xy 


(f) G (f) 
xx yy 


o < y < +1 
— xy — 


(4.23) 
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Whereas the cross-correlation coefficient indicates the degree 

of wide-banded dependency of one signal on another for a given time 

lag, the coherence indicates the degree of narrow-banded frequency 

dependency of the two signals on one another. It is analogous to 

computing the cross-correlation in narrowly filtered frequency bands. 

Note that the theoretical coherence is nonzero only if there is some 

degree of fixed phase information between x(t) and y(t). Therefore, 

2 • 

there would be no coherence (y = 0) between two signals, whether 

xy 

random or periodic, if they were totally uncorrelated over a long 
period of time. 


4.8 Partial Coherence Function 
In some instances a system may exist which has multiple input 
signals that combine to yield a single output signal. If these inputs 
are partially dependent on one another, then the ordinary coherence 
function cannot correctly yield the amount of relationship between 
any one single input and the output. 

Consider the hypothetical acoustic system of Figure 4.1. The 
two independent source mechanisms s^(t) and S 2 (t) partially combine 
with one another to yield emission signals x(t) and y(t). These 
emissions come from two separate locations and are sensed simultane- 
ously by the far-field microphone as z(t). It is desired to have 
knowledge of the amount of contribution each of the input noise 
source mechanisms s^(t) and S 2 (t) make to the output signal z(t) . 

Such knowledge is useful for locating the root cause of a particularly 
offensive sound from a machine or for determining the nature and 
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Figure 4.1 Hypothetical system demonstrating the need for analysis involving partial coherence 
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extent of a subtle noise sou-ce mechanism such as exists with the 
TE noise phenomena. 

2 

Suppose the ordinary coherence y is taken between x(t) and 

xy 

z(t). An appreciable amount of coherence would be indicated whether 
source mechanism Sj(t) were responsible or not. This is due to the 
partial dependence of x(t) on input mechanism 82 (f). 

The "partial coherence" is a more advanced type of coherence 
function that can overcome this problem in the following manner: 

If the phase coherent effects of y(t) were removed from both x(t) 
and z(t), the correct relationship between s^(t) and z(t) could be 
established (Figure 4.2). 

Thus, if all sound emission inputs of a system could be taken 
into account, then the actual degree of coherence between an acoustic 
source mechanism and the output sound at the far-field observer’s 
microphone can be computed. Often, however, the capability exists 
for analyzing only a limited number of channels. In that case the 
partial coherence will yield a result that, although a better estimate 
than the ordinary coherence, is still not the true value of the 
degree of fixed-phase relationship between a noise source mechanism 
and the received noise output. 

For a three-input/single-output system, the equations for 
computing the partial coherence from cross-spectral and auto-spectral 
values are given [3, 5, 19, 52]. They are so written as to yield the 
partial coherence between the input 1 and the output 4 with the phase- 
coherent effects of inputs 2 and 3 removed: 
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5. DATA COLLECTION 
5.1 Recording at NASA 

Test run data were collected at the NASA facility using the 
system shown in Figure 5.1. All data were recorded on the Honeywell 
model 96 14-track wideband FM tape recorder. It conformed to IRIG 
Group I specifications and featured flat response to 40 kHz when 
operated at the 60 inch-per-second speed used during the tests. 

A Lockheed model 417D FM tape recorder was used to bring data 
back to the North Carolina State University computing facility for 
analysis. Selected transducer combinations were dubbed from the 14 
tracks available on the NASA Honeywell machine onto the 4 tracks of 
the North Carolina State University Lockheed machine. The Lockheed 
machine had been modified to give flat record/playback response from 
dc - 10 kHz within ±0.5 dB with some sacrifice of signal- to-noise 
ratio when operated at 15 inches per second. Playback response of 
the Honeywell machine was also dc - 10 kHz at 15 inches per second. 
Dubbing speed for both machines was 15 inches per second, thus dc - 40 
kHz test run data response was possible after a simple 4:1 scale-up 
of frequency on computer analysis outputs. 

5 . 2 Data Collection and Analysis Facility 

To collect the data for digital analysis, it was necessary to: 

1) Filter the played-back data in a manner compatible with 
FFT and digital analysis requirements. 

2) Digitally sample the data at a rate compatible with FFT 
and digital analysis requirements. 

3) Qualify the data as suitable for subsequent analysis. 
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Figure 5.1 


Data collection chain of equipment used at the NASA Langley ANRL facility 
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4) Store the data in a mass storage/retrieval system for 
convenient use at a later date. 

The equipment used at the North Carolina State University 
facility for data collection and analysis is listed below: 

1) Lockheed FM Tape Recorder, Model 417D, Serial No. 800. 

2) Rockland Dual Hi/Lo Filter (two sets), Model 1042 F-03, 

Serial Nos. 1147709 and 1147710. 

3) Biomotion Waveform Recorder, Model 1015, Serial No. 3139. 

4) Digital Equipment Corporation PDP 11/40 Minicomputer, 

Serial No. WM 0106412. 

5) International Business Machines System /7 and System /370 
computer systems. 

Figure 5.2 illustrates the system configuration of the North 
Carolina State University data analysis facility described below. 

5 . 3 Digital Data Collection Procedures 

Various combinations of the four available data channels were 
chosen to be analyzed from the dubbed Lockheed recorder data. The 
selected data were then played back at the North Carolina State Univer- 
sity computing facility through the Rockland filters which were used 
primarily to prevent anti-aliasing when sampling (Section 6.1). A 
secondary function of the filters was to reject undesirable low 
frequency components up to 100 Hz (400 Hz in "NASA test-run time" 
due to the 4:1 recording speed change). The filters were also used 
on occasion to band-pass filter the data to verify the performance 
of the digital filter incorporated into the program CIRCXC0R (Appendix 
11.4). 

The data were then sampled by the Biomation analog-to-digital 
converter (ADC) which contained a 4096, 10-bit word buffer memory. 
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The waveform recorder could be set to collect onu, two, or four 
channels of data with corresponding memory allocations of 4096, 

2048, or 1024 data points (words) per channel. Individual adjustment 
for zero offset and input voltage scaling for each channel was avail” 
able so that the maximum use of the 1024 (2^) level dynamic range 
of the ADC in each channel could be achieved. The sampling rate 
was also set on the Biomation. 

A command from the PDF 11/40 minicomputer operator's console 
initialized the sampling procedure. Then the Biomation would sample 
data from its selected inputs until the 4096 word buffer memory would 
fill. 2he memory contents were then transferred to a disk memory 
connected to the PDP 11/40. The Biomation could then fill its memory 
again. This process was repeated until a specified number of Bioma- 
tion "dumps" were collected in the disk memory. 

The data could be recalled and plotted on a video graphics display 
screen. This was typically employed to insure that the data were of 
sufficient amplitude to yield good dynamic range results when later 
used with the signal analysis software. It was also used to detect 
gross collection errors such as signal drop-outs between recorded 
data runs, evidence of uad electrical connections, filter settings, 
overloads, etc. 

Additional operator console commands edited the Job Control 
Language (JCL) information which was then joined with the stored data 
dumps. The final command to the minicomputer sent out the JCL followed 
by the data as one large file to the IBM System/ 7 computer. The 
Systera/7 was used as a time-sharing data link between many on-campus 
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facilities at North Carolina State University and the Triangle 
Universities Computation Center, or "TUCC." TUCC is a large computing 
facility jointly operated by the three large universities in the 
region: North Carolina State University, Duke University, and the 

University of North Carolina. 

When the proper conditions permitted, the System/7 computer 
transmitted the JCL and data to TUCC. The JCL accompanying the data 
then directed the IBM System/370 computer at TUCC to identify and 
store the data for later use. 

Punched card versions of the computer programs listed in 
Appendices 11-2 - 11.4 were submitted to TUCC from on-campus computer 
terminals. These programs Included the JCL which identified the 
proper data to be analyzed. The results of the computations were 
then returned on paper printouts and in plotted form. 
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6. DIGITAL ANALYSIS PARAMETERS 

6. 1 Sampling and Frequency Parameters 

Analysis of data by digital means using the discrete finite 
forms of the functions outlined in Section 4 required sampling and 
conversion of the signals from analog to digital form. After storage, 
computer versions of the functions (Appendices 11.2 - 11.4) were 
utilized. The qualification of the data and the sampling requirements 
prior to analysis are herein discussed. 

Representative segments of a signal are used for data analysis 
due to time and storage space restrictions. These segments are known 
as time histories or "sample records." The sampling of the data at 
time intervals "h" for a total number of sample points "N" results in 
a sample record period of T = Nh seconds. It is this perior which is 
used by the discrete finite functions of Section 4. 

The resolution of the frequency bands which result from using the 
discrete finite Fourier transform (Equation 4.9) is then given by: 

BW = T = m ’ Hz * (6,1> 

The Shannon Sampling Theorem [32, 35] states that at least two 

samples are required per period of the highest frequency of interest 

"fjj" contained in an analog signal if the signal is. to be adequately 

represented by the sample record. Thus, all frequencies or Fourier 

harmonic components in a signal up to f^ are satisfactorily described 

by the samples, but those greater than f are not. Instead frequen- 

N 

cies greater than f N are mistaken for frequencies lower than f^ and 
are the basis for calling f^ the "Nyquist folding frequency." 
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To prevent this "aliasing," any frequency components in the 
analog data are filtered out before sampling takes place. The low- 
pass filter must therefore have a very sharp folloff above f^ and must 
minimize any alterations of the signal below f^, thus imposing strict 
performance requirements. The Butterworth-type filter set used in 
this study had a rolloff of 48 dB/octave. 

It follows from Shannon's Sampling Theorem that the sampling rate 
”f M must be at least twice that of the cutoff frequency f^, or: 


h = — <_ , sec /sample, 

s N 


( 6 . 2 ) 


Thus, for the 0-10 kHz analysis range which was most commonly 
used for this study the low-pass anti-aliasing filters were set for 
10 kHz and the sampling rate f = 20 kHz (or h = 0.05 msec). A 
higher sampling rate would have resulted in less chance for aliasing 
errors but the resulting loss of resolution would have increased any 
bias errors in the spectral density functions (Section 7.2). Also it 
was learned before any anal yd c took place at North Carolina State 
University that all discernible trailing-edge noise activity occurred 
below 6 kHz. 


6.2' Dynamic Range 

The data was digitized with a 10-bit binary analog-to-digital 
converter (ADC). This gave a dynamic range of ±512 levels or 1 part 
in 1024, Expressed in decibels the maximum signal- to-noise ratio or 
"S/N" possible if the total dynamic range were used is approximately 


[ 83 : 
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S/N ' 20 log [2(2 10 - ])] 

; 66 dB , (6.3) 

The fast Fourier transform subroutine used in this study (Sec- 
tions 4.3 and 11.5) had a dynamic range in excess of 100 dB [17]. The 
anti-aliasing filters used specified a S/N of at least 85 dB. The 
modified FM tape recorder employed had at best about a 45 dB S/N and 
was thus the limiting link in the entire data analysis chain. 

6 . 3 Data Qualification 

"Stationarity" of the data had to be justified before any finite 
Fourier transform functions could be computed; otherwise, results 
would have been in serious error [6, 11, 17], To be classified as 
"weakly" stationary, the mean value or "first moment" and the auto- 
correlation or "joint moment" of the data (Section 4.1) should be 
time invariant. That is to say the data points from any given sample 
record should produce similar values for the auto-correlation function 
and overall mean as the data contained in another sample record of , 
the same signal. If there are trends present in the mean and auto- 
correlation, the data is "nonstationary." 

"Strongly" stationary data are such that all possible moments and 
joint moments are time invariant. Since it is impossible to check all 
moments for trends, the establishment of weak stationarity usually 
justifies the assumption of strong stationarity. 

The data used in this study was assumed to be stationary based on 
knowledge of the wide-band random nature of the phenomena. Also, 
repeated visual checks for trends in the digital sample records were 


made throughout the data collecting process using a videl graphics 
terminal. Figures 6.1 - 6. A illustrate respective time and spectrum 
plots for representative near-field and far-field dynamic pressure 
signals. 

6. A Weighting Function 

Because the data involved were steady-state random noise, there 
was an abrupt truncation of values at the beginning and end of the 
time series sample records. These instantaneous changes of value from 
zero to nonzero are viewed as dirac delta functions by the finite 
Fourier transform. The resulting wide-band spectra associated with 
dirac delta functions (sometimes called ’’spike spectra") add to the 
valid data spectra, causing errors in both spectra amplitude and 
shape. 

To reduce these errors, a time-series weighting function was 

employed. Typically this weighting function has a smooth shape with 

tapering endpoints which are made to coincide with the endpoints 

of the time series. The effect of multiplying the weighting function 

and time series is to gradually reduce the data values to zero as the 

endpoints are approached. Figure 6.5 shows the use of the Hanning 
2 

(sin ) function which was employed for the computer programs used in 
this study. More detailed information on this and other weighting 
functions may be found in the literature [2, 6, 11, 17, 28, 32], 


GXX SPECTRUM LEVEL (RELATIVE DB) £ DIGITAL AMPLITUDE 


RUN 307. KUL.l (NEAR-FIELD SIGNAL), 0.05 MSEC SPLRT 
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SAMPLING INTERVAL. (SEC) 


0.08192 


;ure 6.1 Single representative sample record of a near-field 
pressure signal from a TE Kulite 
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Figure 6.2 Representative auto-PSD of the near-field pressure signal 
of Figure 6.1 after 64 such sample records were zero- 
meaned, Banned, Fourier transformed and averaged 
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7. ERROR ANALYSIS 

Statistical errors due to the use of discrete finite Fourier 
transform methods in signal analysis have been thoroughly reviewed and 
are well documented in the literature [2, 6, 28, 32 1. Such errors 
result from the inability of a finite length time history of an 
event to exactly represent the statistical behavior of the event for 
all time. Thus, some knowledge of the extent of the error is neces- 
sary . 

7. 1 Random Errors in Spectral Density Functions 

Spectrum ensemble averaging was used in this study to reduce the 
majority of statistical errors. The normalized standard (random) 
error as defined for the PSD functions is given in reference [6] as: 

er = (B T )“** (7.1) 

e t 

where: B e = W/T, effective line BW, (7.2) 

T = T x m, total time period, (7.3) 

m is the number of spectral ensemble averages, T is the time period 
of one sample record, and W is the factor relating the bandwidth of 
the time window weighting function employed to the bandwidth of a 
rectangular window. The use of the window weighting function is 
described in Section 6. A. 

Substituting relationships (7.2) and (7.3) into Equation (7.1) 
yields: 

er = (W x m) 


(7, A) 


The error in percent is stated as 100(er). The error in dB is 
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given by: 

er dB = 10 Iog 10 ^ 1 er ^ ’ (7.5) 

For this study Hanning weighting was employed (W = 1.5). The 
number of averages taken was either in = 32 or n = 64. Substituting 
W = 1.5 and m = 32 into the above equations yields a worst case error 
of 14.4 percent or +0.59, - 0.68 dB which was marginally acceptable 
for the purposes of this study. 

7 . 2 Bias Errors in Spectral Density Functions 

Whereas the normalized standard error discussed above depends 
solely upon data acquisition and analysis parameters, the amount of 
"bias'* error depends partially upon the nature of the data itself. 
Thus, it cannot be predicted before analysis. 

The nature of the bias error in the PSD estimates is such that 
the "peakedness" of the function is reduced. Thus, sharp "peaks" in 
the true spectra are underestimated or reduced, and narrow "valleys" 
are overestimated or increased when the PSD estimate is made [6]. 

The broadbanded character of the trailing-edge noise phenomena 
reduced the bias errors to a minimum. Because the standard and bias 
errors are most severe for the PSD estimates on which no corrections 
were made [6], no attempts were made at correcting the errors in the 
transfer function phase as computed by the program COHERENC (Appendix 


11 . 2 ). 
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7.3 Bias Errors in Coherence Functions 
The coherence function contains another type of bias for which 
corrections were made. By the nature of its definition the coherence 
function (Equation A. 22) always has a value of unity when only one 
spectrum ensemble average is taken. Increasing the number of averages 
reduces this positive bias and yields a better estimate. The bias 
corrected coherence is computed from [32]: 



» - v, 2 > 

2Wm 


(7.6) 


- 2 

where y q is the averaged uncorrected estimate for the coherence and W 
and m are as defined previously. 

The degree to which the corrected value of coherence at any one 
frequency can be relied upon may be stated in terms of a confidence 
interval. For a chosen percentage of reliability this interval indi- 
cates the upper and lower limits within which the coherence must fall. 

7.4 Confidence Interval for Coherence Functions 
The 95 percent confidence interval for the bias-corrected coher- 
ence estimate is computed by using a method developed by Benignus as 
outlined in reference [32]: 


A 2 

y upper ~ tanh “ b + s (za)] (7.7) 

y lower tanh “ b " s (m) 3 (7.8) 


where: 


z 

s 


= 0.5 fcn [(1 + y) / (1 - y) ]. 

= b^l - 0.004 (1 * 6y2 + °* 22) ] , 
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b = p / [2Wm - 3 p + 1], 
za = 1.96, 

p is the number of coherence function inputs, za is the one-sided 
100 a percentage point of the normal distribution of the data ^nd W 
and m are as defined previously. For 100 a = (100 - 95) = 5 percent, 
a/2 = 0.025 and thus za = 1.96 (as found in statistics tables). 

A plot of the 95 percent confidence limits as a function of the 
number of spectral averages and the original coherence estimate with 
W = 1.0 is given in Figure 7.1 [17]. 

7.5 Errors Due to Time Delay in Spectral Quantities 
' An additional error that affects all spectral quantities is the 
negative bias caused by signal path time delays. These delays which 
occur, for example, between separated transducers are best compensated 
for by time alignment of the data sample records prior to processing. 

Where time alignment is not possible, references [4, 23, 36] give 
detailed information on determining the amount of the error. Basic- 
ally, the PSD and coherence functions are in error by an amount: 


G(f) . 
G(f) - 



(f) 

(f) 



( 7 . 9 ) 


* - 2 

where G(f) and y are the computed averaged estimators with no time- 
xy 

2 

delay compensation, G(f) and y are the "true" values, T is the 
time length of one sample record, and x is the time delay. Although 
these simple equations are useful for obtaining approximate amplitude 


corrected values of G(f) and y yt they cannot compensate for the 

xy 
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increase in bias and standard (random) errors discussed previously. 

For determining the total extent of these errors, the reader should 
consult the above-mentioned literature [4, 23, 36 1. 

Another source of time-delay error arose from the refraction of 
the acoustic wave as it traveled through the shear layer of the open 
jet wind tunnel to the far-field microphone. Cross-correlation and 
cross-spectral measurements between a transducer at the trailing edge 
of the airfoil and the far-field microphone were thus potentially 
subject to this additional time delay (Figure 7.2, path A). 

A study of this effect and the derivation of equations for 
obtaining correction factors has been performed by Amiet [l]. The 
freestream Mach number of the flow in this study did not exceed 0.20 
(225 ft/sec), and the far-field measurement microphones were on a line 
perpendicular with the trailing edge and the air flow. Based on 
conclusions reached in Amiet.' s paper, it was determined that the 
acoustic pressure amplitude and phase corrections were not warranted. 
The amplitude corrections were less than 0.5 uB which is less than 
the normalized standard error. 

As an example, the worst case time delay error (not including 
refraction effects) is computed for 12-foot distant microphone M.3 
(Figure 7.2, path B) . 


t 


d _ 12.0 
c ~ 1125 


10.7 msec. 


(7.10) 


T 



1024 

20,000 


= 51.2 msec. 


(7.11) 


-•^) 2 “ (1 - 10. 7/51. 2) 2 - 0.628 


(1 


(7.12) 
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or expressed in decibel form: 

10 log j 0 (0.628) » -2.03 dB (7.13) 

r?here d is the spacing, c the speed of sound, f the sampling 
frequency, N the number of sample points per time record, and T and i 
are as defined previously. The quotient i/T indicates that the two 
sample records are effectively offset in time by 20.8 percent. 

Time alignment by data point shifting was not employed because 
of the limited capabilities of the North Carolina State University 
analysis facility. In order to maintain the maximum frequency resolu- 
tion possible out of the available fixed sample record length of 1024 
points per channel, no points could be sacrificed. Shifting the data 
would negate the negative bias errors due to time delay, but at the 
expense of an accompanying increase of standard (random) errors due 
to the lesser number of data points to be Fourier transformed. Also, 
comparisons between spectral density curves would become difficult 
if the frequency resolution were different among curves. 
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8. EXPERIMENTAL RESULTS 

Series I data brought back from the NASA test facility for 
analysis at North Carolina State University was only sufficient for 
a preliminary evaluation of the trailing edge noise phenomena. 
Subsequent analysis of data from Series II tests should yield substan- 
tially more detailed analysis results. Thus, the information present- 
ed herein primarily represents the methodology developed and the 
characteristics of the TE noise phenomena which were observed. 

8. 1 Microphone Calibration 

Calibrations for the one-half inch microphones were performed 
using a Briiel and Kjaer type 4220 Pistonphone producing 124.0 ±0.2 
dB SPL. 

Calibrations for the Kulites were performed after they were 
flush-mounted in the airfoil. A General Radio model 1562-A Sound 
Level Calibrator equipped with a seal was verified as producing 112.9 
dB SPL when pressed against a flat surface. It was then pressed 
over each Kulite in turn to yield an in-place "flush" calibration. 

All results in this thesis, however, are presented in "relative 
dB" meaning that the dB scale for each plot is not calibrated against 
any reference standard. Thus, pressure levels between plots may not 
be compared, but pressure level comparisons on any one spectrum plot 
are permissible. 

8.2 Validity of Far-Field Microphone Measurements 

Before analysis took place, a <<ijSck was made to insure that the 
closest "far-field" microphones were indeed in the acoustic far-field 
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of the radiated TE noise. An assumption was made (which was later 
justified) that a line source of Independent dipoles radiating wide- 
banded noise would have a SPL falloff with distance similar to that 
for a spherical wave source. Therefore a check could be made to see 
if Ka >>1, say 10 as a lower limit. For the 4-foot distant micro- 
phones, then a = 4 and: 

Ka >> 1, = 10 say, (8.1) 


K 


f 


10 = _10 
a 4 


2.5 ft" 1 . 


Kc = 2.5(1125) . , 
2tt 6.28 ~ 


the wavelength constant, 


Hz. 


( 8 . 2 ) 


Thus, due to the high-pass filtering at 400 Hz as set on the anti- 
aliasing filters, all observed spectral phenomena were assumed to be 
valid far-field measurements. 

The nearness of the 4-foot microphones to the jet did not appear 
to produce auto-spectra that were significantly different from those 
of the 12-foot distant microphone M.3. However, the noise floor of 
the FM tape recorder was reached in some instances, and so deviations 
between plots of 3-4 dB in frequencies above about 5 kHz were occasion- 
ally noticed in some auto-spectra. 


8.3 Initial Checks Using a Cylindrical Rod 
As an initial check on the test setup, the airfoil was removed 
and an 18-inch long steel rod of 0.375-inch diameter mounted in its 
place at the trailing edge position (Figure 8.1). The rod provided a 
dipole line source whose known radiation characteristics could be used 
to verify correct operation of the analysis system. 





f Figure 8.1 Initial check run using a rod mounted at the TE position 
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For a flow rate of 125 feet per second (500 rpm) , the peak 
frequency of the alternate vortex shedding for the rod was calculated 
using [25]: 

Su 

f = (8.3) 

where S is the subsonic Strouhal number for a cylinder, u^ is the 
freestreara velocity, and d is the rod diameter. Using the above 
values and a Strouhal number of 0.2, the peak frequency was calculated 
to be 800 Hz. 

Using the two 4-foot distant microphones, M.2 and M.5, the actual 
auto-spectra, cross-spectrum, phase and coherence were computed by the 
program COHERENC (Appendix 11.2) (Figures 8.2 - 8.6). The actual peak 
of 780 Hz was in good agreement with the above predicted 800 Hz. The 
180° dipole nature of the alternate vortex shedding is evident in the 
phase plot. The cross-spectrum and coherence also give evidence of the 
high ratio of fixed-phase common signal to uncorrelated signal present 
in such a phenomena. 

8.4 Test Conditions 

A generalized airfoil test configuration was used during the 
initial phases of analysis. No modifications such as the trailing 
edge extensions used in Series II tests were made to the basic airfoil 
as shown in Figure 2.2. 

The airfoil was tested at angles of attack a = 0°, ±5°, ±10°; 
but for any given airflow velocity no differences could be detected in 
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8.4 Cross-PSD of a tone generated by the 0.375 inch diameter 
rod mounted at the TE position as detected by M.2 and M.5 


w 500 RPM. RUN 21S , fl0A*0, ROD AT T.E.; X-MIC.2, Y-MIC.S; 6S AVGS. 

90.0 

55.0 
0 

-95.0 
-90.0 
-135.0 
-180.0 
-225.0 
-270.0 

0 1000 2000 3000 9000 5000 6000 7000 8000 9000 10000 

FREQUENCY (HZ) 




1 ' 1 '> *" l ‘ ”‘i " 1 

Returns to 0° on this 

and other phase plots caused by protective 






(La 


/ 

c< 

jaerlcal c 
atiputer pr 

ut-offs ir 
ogram 

the 




■ 







1 







■ 




















1 

1 



■ 

1 

■ 


■ 






1 

| 



■ 

■ 

II 



\ 

'1 




1 

■ 




■ 




■ 

■ 




■ 

■ 


■ 


Figure 8.5 


Phase between dynamic pressure field signals originating 
at M.2 and M.5 for the 0.375 inch diameter rod mounted 
at the TE position 
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the observable phenomena. Thus, all results reported in this section 
dea) with the airfoil mounted at a = 0°. 

Data were recorded for freestream velocities ranging from 52 to 
225 feet per second. However, a single data run. Run 307, was the 
single most widely analyzed collection of data. It consisted of 
recordings of the phenomena at the highest flow rate available (225 
feet per second, or 900 rpm) with the airfoil in the general test 
configuration. It was thought that at this freestream velocity the 
noise phenomena, if it existed, would be most strongly excited and 
therefore easily analyzed. It was assumed that the blunt leading 
edge radiated no noise. 

8 . 5 Auto-Power Spectral Density Measurements 
Figure 8.7 depicts the auto-spectrum of far-field microphone M.3. 
Immediately noticeable is the 4-5 dB hump in the 2.0 - 4.0 kHz band of 
frequencies. This auto-spectrum is characteristic of all far-field 
microphone positions for the airfoil general configuration. 

The wide-banded background noise floor or BNF (the downward 
sloping curve without the inclusion of the hump) is due to the aero- 
dynamic noise of the open jet and not the anechoic facility itself. 

Up to a frequency of about 5.0 - 6.0 kHz, the pressure level of the 
BNF can be predicted for the 225 foot-per-second freestream velocity 
by: 

p(f) = k(f)“ 1,125 (8.4) 




where k is a constant particular to each plot. Good agreement was 
found among the auto-spectra of all three far-field microphones (M.2 
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gure 8.6 


Coherence between dynamic pressure field signals 
originating at M.2 and M.5 for the 0.375 inch 
diameter rod mounted at the TE position 
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Figure 8.7 Auto-PSD of far-field dynamic pressure signal from 
12-foot distant microphone M.3; also representative 
of 4-foot distant microphones M.2 and M.5 
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M.3, M.5), The flattening out of the auto-spectral curves above 
5.0 - 6.0 kHz is possibly due to the TE phenomena, the noise floor 
of the FM recorder, or both. 

The auto-spectra of the TE Kulites (K.l, K.2, K.5) also display 
the characteristic "hump" in the 2.0 - 4.0 kHz band of frequencies 
(Figure 8.8). 

The auto-spectrum of the upstream Kulite (K.3) has the same flat 
BNF as the TE Kulites but without the hump (Figure 8.9). 

A check run was made to insure that the characteristic hump 
observed in all far-field and TE transducer positions was due solely 
because of the airfoil and not because of the sideplates or nozzle. 

The nozzle and sideplate assembly without airfoil were tested at the 
same freestream velocity of 225 feet per second as the standard 
configuration. The far-field auto-spectra appear the same as Figure 
8.7 with the same level of BNF except for the lack of a hump. 

Another check was made with the airfoil held in place over the 
nozzle but without the sideplates. The far-field auto-PSD is virtual- 
ly unchanged from the standard configuration, whereas the near-field 
auto-PSD of the TE Kulite shows a slightly more pronounced hump 
(Figures 8.10 and 8.11). 

8.6 Cross-Spectral Function Measurements 

To clarify the nature of the hump the coherence, cross-PSD, and 
averaged phase were computed for the symmetrically opposing equi- 
distant far-field microphones M.2 and M.5. The coherence indicates 
a strong relationship between the pressure fields which existed at 
the microphones in the range of hump frequencies (Figure 8.12). When 
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Figure 8.9 Auto-PSD of near-field dynamic pressure signal from 
upstream Kulite K.3 
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coupled with the knowledge that there is no visible hump in the 
far-field auto-PSD spectra when the airfoil was not present nor in 
the upstream Kulite K.5 signal when the airfoil was present, this 
leads to the conclusion that the phenomena responsible must be at or 
near the TE. 

The coherence also indicates a smaller secondary emission which 
extends from about 0.8 to 2.3 kHz. Both emissions average 180 degrees 
out-of-phase which indicates that the TE noise phenomena is "dipole- 
like'' ^Figure 8.13). 

The two humps are clearly visible in the cross-PSD plot of Figure 
8.14. Interesting to note are the values of the slopes on the upper 
frequency side of each hump. These slopes are frequency related to 
the pressure as was the far-field BNF and may be found by generalizing 
Equation (6.4): 

p(f) = k(f ) 3 (8.5) 

where k is a constant particular to each test-run plot and 8 is the 
power constant or slope of the 20 log (p) curve. Values of -1.5 for 
the lower hump and -3.0 for the upper hump were obtained. This results 
in a 2:1 slope ratio for the respective upper and lower humps. 

The coherence and cross-PSD curves for symmetrically opposing TE 
Kulites K.l and K.5 also show sigr-3 of a double-humped phenomena 
(Figures 8.15 and 8.16). However, the frequency range of the lower 
hump has shifted downward and the phase plot. Figure 8.17, shows that 
the lower hump is, on the average, comprised of in-phase components, 
contrary to the far-field results. These may indicate a totally 
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Figure 8.12 Coherence between signals from symmetrically opposing 
M.2 and M.5 
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Figure 8.14 Cross-PSD between signals from symmetrically opposing 
M. 2 and M.5 
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Figure 8.16 Cross-PSD between signals from symmetrically opposing 
TE Kulites, K.l and K.5 
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Figure 8.17 Averaged phase between signals from symmetrically 
opposing TE Kulites, K.l and K.5 


83 


different type of mechanism although the fairly low coherence reveals 
that even in the near-field close to the TE, the lower hump mechanism 
on the upper surface is not very strongly related to any that may 
exist on the lower surface. 

Again using Equation (6.5) to find the slopes 8 of the upper and 
lower frequency humps, it is learned that there is a 2: 1 slope ratio 
for the respective upper and lower humps (Figure 8.16). 

In both cases of symmetric transducer pairs, the averaged phase 
of the pressure-field components is random in frequencies above the 
characteristic primary (upper frequency) hump. 

A 2.0 - 4.0 kHz cross-correlograrn of the opposing TE Kulites K.l 
and K.5 as computed by the program CIRCXCOR (Appendix 11.4) also ex- 
hibits the out-of-phase structure of the pressure field at the TE 
(Figure 8.18). 

Indications of a strong relationship between the pressure fluc- 
tuations near the TE (K.l) and at a position 0.275 inches upstream 
(K.3) are given by the coherence in Figure 8.19. The averaged phase, 
Figure 8.20, exhibits the characteristics of a pure time delay. This 
distance-time information was used to obtain the mean convection 
velocity u c of the turbulent BL air flow by employing the equation: 


u 
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360d 
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( 8 . 6 ) 


where d is the transducer separation distance in feet and 8 is the 

slope of the phase in degrees/Hz. Taking 8 = 0.C678°/Hz and d = 

(0.275/12) feet yielded a convection speed of 122 feet per second for 

u = 225 feet per second. 
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Figure 8.18 Cross-correlogram of symmetrically opposing TE Kulite 
K.l and K.5 after digital band-pass filtering from 
2. 0-4.0 kHz 
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Figure 8.19 Coherence between signals from TE K.l and upstream K.3 




The amount of time t necessary for a nonacoustic disturbance to 


propagate from K.3 downstream to K.l is found using: 



where u is some mean airflow velocity and d is as defined previously. 
Using d = (0.275/12) feet and u c = 122 feet per second, a value of 
T = 0.188 milliseconds is obtained for the propagation time. 

The 2.0 - 4.0 kHz hump filtered cross-correlation coefficient 
plot for K.l and K.3 is shown in Figure 8.21. The timing of the 
first large positive peak at -0.150 milliseconds agrees well with 
T = 0.188 milliseconds computed above. It is likely that if finer 
resolution of the time delay axis were obtained (i.e., a faster 
sampling rate), a more accurate value approximately equal to 0.180 
milliseconds could be attained, as indicated by the dashed line. 

The coherence (Figure 8.19) and to a much lesser degree the 
cross-PSD (Figure 8.22) exhibit two strong humps in the 0.5 - 3.0 
and 3.0 - 6.0 kHz bands. This contrasts with the 2.0 - 4.0 kHz auto- 
PSD hump seen in both near- and far-field pressure spectra (not in- 
cluding K.3 which exhibits no hump at all). 

An interesting point to note at this time is how the center 
frequency of the characteristic 2-4 kHz hump mentioned above may be 
approximated by the use of the Strouhal relationship between the 
convection velocity u c and the characteristic blunt TE dimension of 
0.10 inches (see Figure 2.2) [37]. Substituting u c = 122 feet per 
second, d = (0.10/12) feet and S - 0.2 into Equation (8.3) results in 
a frequency of 2928 Hz. No justification is given for using 0.2 as 
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the proper Strouhal number for the geometries involved other than 
its acceptance as a value to use for most structures under general 
subsonic flow conditions. 

It seems likely that the suggested frequency selective Strouhal 
phenomena may be responsible for the "notch" in Figures 8.19 and 8.22. 
It also appears likely that were it not for the TE phenomena, whatever 
its true mechanism, the coherence would lack the "notch" and smoothly 
taper as denoted by the dashed line. The taper then gives evidence 
that the high frequency energy in a small turbulent volume of air 
convecting downstream dissipates much more rapidly than the low 
frequency energy in the volume. 

Thus it is possible to explain this "notch" on the basis of 
frequency-related elemental volume energies. If the TE noise radia- 
tion mechanism were selectively converting the turbulent kinetic 
energy of an elemental volume of air to acoustic ra4} *tion energy on 
the basis of frequency, then that energy transformation must have 
come from a change in the pressure/density structure of the volume of 
air. Thus the pressure field of the volume as it passes over K.l 
is not the same as it was when it passed over K.3, and hence the 
"notch" appears in the cross-spectral quantities. 

Spanwise communication of the turbulence is not nearly as great 
as in the strearawise direction. The random phase and low coherence 
as shown in Figures 8.23 and 8.24 indicate the lack of any definite 
relationship that exists between upper surface spanwise Kulites K.l 
and K.2 which are separated by 0.135 inches. 

The relationship between the pressure field at TE Kulite position 
K.5 and 12-foot-distant M.3 on the same side is typified by Figures 
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Figure 8.23 Averaged phase between signals from upper surface span 
wise TE K.l and K.2 
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8.25 - 8.27. Again the randomness of the phase and accompanying lack 
of any appreciable coherence reveal that the local near-field pressure 
and the far-field pressure are not intimately related. The cross-PSD 
plot does not yield any additional information. 

A similar result was obtained for the near-field/far-field 
transducer pair consisting of TE K.l and 4-foot-distant M.5 (Figures 
8.28 and 8.29). The cross-correlation coefficient plot for this pair 
(Figure 8.30) does, however, imply that: 

a) Some degree of dependency does indeed exist. 

b) There is some sound which travels the required 3.55 
milliseconds from the TE K.l to the 4-foot M.5 as 
computed using Equation (8.7) where u is replaced 
by c = 1125 feet per second. 

The cross-PSD plot for this pair is similar to Figure 8.27. 

8 . 7 Detection of Near-Field/Far-Field Relationships 
via Subtraction Schemes 

In an attempt to enhance the capability of the analysis system to 
establish a definite relationship between near-field Kulite and far- 
field microphone signals, two computational schemes were employed. 

The first involved the computer subtraction of signals from 
Kulite pairs via the program COHERENC after appropriate amplitude 
scaling. The intention was to take advantage of the 180° dipole-like 
nature of the pressure field in the 2.0 - 4.0 kHz hump frequency band 
and perform a point-by-point subtraction of the time data series of 
one TE Kulite from the time data series of its symmetrically opposite 
twin. 

A factor to amplitude-match the two channels to be subtracted 
was derived from the overall rms signal level as computed via COHERENC 
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Figure 8.25 Averaged phase between signals from near-field TE K.5 
and far-field M.3 
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Figure 8.28 Averaged phase between signals from near-field TE K.l 
and far-field M.5 
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in the 6.0 - 10.0 kHz band — i.e., those frequencies well away from 
the hump phenomena. 

The signal information in the out-of-phase hump frequencies was 
predicted to double in amplitude (6 dB) . The lower frequency in- 
phase information, which represented the insignificant life components 
of the pressure field, was expected to subtract out, and the higher 
frequency random phase information would increase by a factor of 1.414 
(3 dB). 

The first attempt at using this method involved K.l and K.5 
versus M.3. The results are shown in Figures 8.31 - 8.34, The auto- 
PSD for M.3 was identical to the general far-field spectrum given in 
Figure 8,7. 

The characteristic hump displayed in the auto-PSD of the subtract- 
ed TE Kulites seems to be enhanced and the relative level of the 
pressure in the lowest frequencies appears to be reduced as compared 
with an unmodified TE Kulite auto-PSD such as Figure 8.8. The cross- 
PSD has a more pronounced hump as compared with Figure 8.27, and 
there seems to exist a smaller secondary hump from 1.3 - 2.3 kHz, 
akin to that shown by the coherence between M.2 and M.5 (Figure 8.12). 
However, the averaged phase is still predominantly random and the 
coherence still does not indicate any near-f ield/far-f ield relation- 
ship. 

A second attempt at enhancing the ordinary coherence, cross-PSD, 
and phase functions involved not only the above-discussed Kulite 
subtraction, but also the subtraction of the symmetrically opposing 
4-foot microphones M.2 and M.5 (Figures 8.35 - 8.38). The subtracted 
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Kulite auto-PSD appears the same as Figure 8.31. The primary hump 
as seen in the auto-PSD of the subtracted far-field microphones 
(Figure 8.25) is markedly improved over that of 12-foot M.3 (Figure 
8.7). Also the smaller secondary hump is more clearly extablished, 
due in part to the better smoothing that results from the greater 
number of ensemble averages. 

The cross-PSD, phase and coherence (Figures 8.36 - 8.38) do not 
show evidence of improvement over any previous results. 

As a check to verify that the subtraction scheme was indeed 
enhancing any of the dipole-like phenomena occurring and that the 
amplitude scaling procedure used was proper, the identical informa- 
tion as in Figures 8.35 - 8.38 was again analyzed. This time, how- 
ever, the time data series of symmetrically opposite transducers 
were added point-by-point. The lack of any humps in either auto-PSD 
or in the cross-PSD (Figures 8.39 - 8.41) is proof that the method 
used was proper. There was no significant coherence and the phase 
again was entirely random. 

8.8 Detection of Near-Field/Far-Field Relationships 

via Partial Coherence 

The second scheme aimed at establishing some dependence of far- 
field sound on the pressure in the near-field involved the use of the 
partial coherence function as described in Section 4.8. 

It was believed on the basis of the span-wise TE cross-spectral 
measurements that the overall far-field sound due to the TE phenomena 
(i.e., the hump) was the result of emissions from many independent 
stream-wise "strips" of turbulent activity near the TE (Figure 8.42). 
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By taking into account more of these independent sources through the 
use of the partial coherence function, it was hoped that some near/far 
pressure-field connection would at last be established. 

In summary, the three-input/one-output program PARTIAL failed 
to show any improvement over the single-input/single-output results 
as computed by the program COHERENC. The results for K.l, K.2 and 
K.3 as inputs with M.5 as the output are shown in Figures 8. A3 - 8.45. 
Apparently not enough of the "strips” of acoustic emission as depicted 
in Figure 8.42 can be accounted for, even when u&ing the multi-input 
program PARTIAL. 


8.9 Verification of the Fifth-Power Hypothesis 
A fundamental precept of the theories of Howe and many others 
mentioned in the introduction is that the overall sound level in the 
acoustic far-field increases as the fifth power of the flow velocity 
[27], i.e.: 
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where: k = a factor dependent upon the "flyover" angle (constant 

for the case of the airfoil and microphone configura- 
tion used in this study), 

P o = the mean air density, 

= the mean turbulent fluctuation velocity of a turbulent 
volume or eddy near the TE, 

u c - the convection velocity. 


Regions of turbulence 
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a) coherent streamwise 

b) much less coherent 
spanwlse 

“=> serai-independent 
"strips” of turbulence 
generating many separate 
TE noise emissions 
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Figure 8.42 Depicting the overall far-field dynamic pressure field 
to be the result of TE noise emissions from many semi- 
independent strips of coherent turbulent flow 
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Figure 8.43 Partial coherence between signals from near-field K.l 
and far-field M.5 with the phase-coherent effects of 
K.2 and K.3 removed 
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Figure 8.44 Partial coherence between signals from near-field K.2 
and far-field M.5 with the phase-coherent effects of 
K.l and K.3 removed 
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Figure 8.45 Partial coherence between signals from near-field K.3 

and far-field M.5 with the phase-coherent effects of K.l 
and K.2 removed 
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M = the Mach number = u /c, 

0.1 

u B the freestream velocity, 
c ■ the speed of cound, 

L ■ the span of the airfoil TE which is "wetted" by the 
turbulent eddies, 

£ * the characteristic turbulence correlation scale 

(depicted in Figure 8.42), 

r = the distance from the TE to the observer position in 
the far-field. 

Since both and u are approximately linearly proportional to 

u , Equation (8.8) becomes: 

00 

p 2 : kp 2 u 5 [if] (8.9) 

r 

where k now has absorbed many other parameters. 

To verify this hypothesis, a series of far-field M.3 auto- 
spectra was computed as the freestream velocity was increased from 
123 to 225 feet per second. The spectrum curves were obtained for 
the test setup with airfoil and without airfoil. After proper ampli- 
tude scaling to insure equally pressure sensitive microphone record- 
ings, the line-by-line difference was calculated between the auto- 
spectrums. This resulted in spectrum values for the hump alone 
without background r- ise from which the overall SPL.of the hump was 
computed (Figure 8.46). 

The log of both the squared pressure and the velocity was taken 
and the values plotted (Figure 8.47). The slope of the line, 6, was 
computed from the data points using the method of least squares. 



Figure 8.47 Determination of the power 6 to which the far-field 
SPL scales 
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The far-field sound levei was found to scale on the freestream 

6 

velocity raised to a power of 8 “ 5.3. This is tending towards u 
which indicates why the dipole-like nature of the hump is observed. 
However, it is in fair agreement with the predicted power of 5.0. 

Hence, the trailing edge noise phenomena cannot be called purely 
dipole-like as the Strouhal relationship might suggest. 

8. 10 Measurement Difficulties and Their Relationships 

to Theory 

It is evident from Equation (8.9) that the acoustic power radiat- 

2 

ing from the TE falls off as (1/r ) as for a spherically diverging 
sound wave from a "compact" source. By "compact" it is meant that 
the size of the source is much smaller than the wavelength of sound 
it radiates. Compactness, then, was verified by both the correlation 
scale and coherence measurements in the streamwise and spanwise 
directions as covered in the previous pages and as depicted in Figure 
8.42. 

Based on the low spanwise coherence measurement (Figure 8.24) 
the spanwise correlation scale of the turbulence "£" must be equal to 
or less than the 0.135 inch spacing of the spanwise transducers K.l 
and K.2. Therefore, more than (L/i) = (18/0.135) = 133 sources may 
have existed along the "wetted" span of the TE. This serves to explain 
why none of the cross-spectral functions used in this study could 
successfully establish a near-f ield/far-f ield relationship between a 
single "source" and the sound in the far-field. When any single 
source region is monitored and compared against the overall far-field 
sound, all other sources are regarded as uncorrelated signal by the 
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analysis programs. Hence the great number of uncorrelated signals 
simply "swamps out" the noise emission from the single source and 
prevents the analysis programs from detecting any fixed-phase 
dependence. Thus, the lack of any definite near-f ield/far-f ield 
coherence indirectly verifies the existence of many independent, 
compact sources. 

Another factor which made the detection of TE noise difficult in 
the ANRL facility was the small size of the Cast airfoil relative to 
the dimensions of a full-scale aircraft wing. As an example, consider 
the wings of a Boeing 747 jumbo jet which have approximately 110 feet 
of wetted trailing edge per wing and a 30-foot chord. Since all flow 
and boundary layer characteristics of the test airfoil were approx- 
imately the same as that for a full-sized during landing approach 

2 

conditions, the factor (Li/r ) of Equation (8.9) implies that there 
are far fewer independent sources for the test airfoil than for the 
full-scale wing. 

If the far-field sound is dependent upon wetted span alone 
and not upon surface area, the sound power level difference may be 
expected to be on the order of 10 log^(l .5/220) = -22 dB. Hence, 
there is much less sound radiated from the test trailing edge as from 
the full-scale trailing edge. Noise from the open jet and side- 
plates is able to effectively mask the radiated TE noise (Figure 8.46) 
except at higher flow velocities causing the detection schemes to 


fail. 
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9. SUMMARY OF RESULTS AND CONCLUSIONS 

1. A miniature strain-gauge pressure transducer ("Kulite") was 
designed and manufactured with outer dimensions of 0.375" L x 0.125" W 
x 0.045" D and a pinhole opening of 0.0135" diameter. The resonant 
Helmholtz frequency of the pinhole and internal cavity was approxi- 
mately 63 kHz. Flat frequency response and phase shift of less than 
1.0° were achieved to well beyond the 10.0 kHz frequency limit of 
data analysis. The nominal sensitivity of the Kulite was 25 mV/psi. 
Several such Kulites were mounted flush with the airfoil surface near 
the TE where they satisfactorily detected the local dynamic pressure 
due to the TE noise phenomena. 

2. A companion two-stage amplifier and strain-gauge power supply 
for the Kulite was designed and constructed that met the operational 
requirements of the microphone. The final version as used at the NASA 
ANRL facility employed an Analog Device 521K operational amplifier 
and an Ithaco model 455 active filter for the respective first and 
second gain stages. The gain of the first stage was fixed at 50 V/V, 
while the gain of the second stage was allowed to vary according to 
signal level and recording requirements. A balanced direct-coupled 
input from the Kulite to the first stage and ac coupling between the 
first and second stages simultaneously provided for the rejection of 
electrical interference and prevented small dc current imbalances 

in the strain-gauge bridge from overloading the final stage output. 
Flat frequency response to beyond 100 kHz was achieved with this de- 


sign. 
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3. A pressure coupler was constructed which provided for the 
calibration and performance verification of the Kulite pressure 
transducers. It utilized a Bruel and Kjaer type 4134 one-half inch 
condenser microphone as a source and a Bruel and Kjaer type 4138 one- 
eighth inch condenser microphone as a reference. Interunit phase 
and amplitude matching of the Kulites was performed using the pres- 
sure coupler. 

4. Error analysis showed that because of the wide-band random 
nature of the signals analysed, the bias error of the spectral density 
functions was of no concern. The normalized standard error was within 
i0.7 dB for the number of spectrum ensemble averages used. The errors 
associated with the refraction of the sound passing through the open 
jet shear layer were insignificant. The method used for the bias 
correction of the ordinary coherence was shown. The error due to the 
time delay between trailing edge and far-field microphones was shown to 
give a constant negative bias of -2.03 dB for the worst case 12-foot 
microphone auto-PSD. 

5. All microphone data were assumed to be stationary and suit- 
able for analysis based on their widebanded random but steady nature 
as verified by visual inspection of the collected data sample records. 

6. Auto-PSD results from trailing-edge Kulites and far-field 
microphones exhibited a single characteristic 3-5 dB "hump’* in the 
2.0 - 4.0 kHz band. 

7. The frequency spectrum of the upstream Kulite exhibited no 
such hump as described in (6) above. It did, however, exhibit the 
same flat random noise spectrum as did the TE Kulites. 


1 16 


8. The slope of the far-field background noise floor spectrum 

_ 1.125 

was seen to vary as f 

9. Streamwise communication of the turbulent pressure field over 
a distance of 0.275 inches was quite significant as indicated by the 
coherence, phase, and correlation computations. A peak coherence of 
0.55 at 1.3 kHz was obtained. 

10. Rather than use the cross-correlation coefficient, the slope 
of the phase between the streamwise Kulites was used to clearly 
establish the turbulent boundary layer convection velocity at 122 feet 
per second for a freestream velocity of 225 feet per second. 

11. The coherence between the streamwise Kulites indicated that 
the coherent energy /Hz per unit volume decreases with increasing 
frequency for a given streamwise separation distance. 

12. The "notch" in the coherence at 3.1 kHz between streamwise 
Kulites coincided with the peak of the characteristic trailing edge 
hump and seems to indicate that the trailing edge noise mechanism 
selectively extracts energy in the "hump frequencies" out of the 
boundary layer turbulence for use in the generation of sound. 

13. Spanwise communication of the turbulent pressure field over 
a distance of 0.135 inches across the trailing edge appeared very 
slight as indicated by coherence and cross-correlation measurements. 

14. There existed a 2:1 ratio of the negative slopes between the 
primary and secondary humps as seen in the cross-PSD plots for either 
pair of symmetrically opposing microphones (K.l and K.5 or M.2 and 
M.5) . 

15. The peak values of the coherence for both the near-field 
symmetrically opposing TE Kulites (K.l and K.5) and the far-field 


m t * , 
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symmetrically opposing 4-foot-distant microphones (M.2 and M.5) 
coincided with the 3.1 kHz peak of the characteristic hump seen in 
the far-field microphone and near-field TE Kulite auto-PSD curves. 

Peak values for the coherence were approximately 0.59 for the far- 
field pair and 0.32 for the near-field pair. 

16. Symmetrically opposing microphones in both the near-field and 
the far-field indicate that, on the average , the upper and lower pres- 
sure fields were in-phase in the frequencies below the primary charac- 
teristic hump, out-of-phase in the hump frequencies, and generally 

of random phase above the hump frequencies. 

17. A computer method for amplitude scaling of the stored 
digital data was devised that allowed the computer to subtract the 
time series of symmetrically opposing transducer pairs to better 
detect the out-of-phase information contained within them. 

18. The Strouhal relationship was successfully used with the 
characteristic turbulent boundary layer convection speed of u £ = 122 
feet per second and characteristic trailing edge thickness of 0.10 
Inches to closely approximate the center frequency of the characteris- 
tic hump. A calculated value of 2.9 kHz was obtained versus the 
observed 3.1 kHz value. 

19. The pressure-field that existed in the vicinity of the trail- 
ing edge appears to be comprised of streamwise strips of coherent 
turbulence which were independent of one another in the span-wise 
direction. Item (13) indicates that the maximum width of such a strip 
is 0.135 inches. Therefore, greater than (18/0.135) = 133 such 
"sources” may have existed across the 18-inch span of the airfoil. 
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20. No spectral relationships or acoustic contributions could 

be established between the pressure fields at the near-field and far- 
field microphone positions, nor could the actual source-coherent area 
of location at or near the trailing edge be determined. 

21. The great number of "sources" which may have existed serves 
to explain why cross-spectral functions could not establish any 
relationship between a single "source" and the sound in the far 
field. In such a situation, the emissions from the other independent 
"sources" would be regarded as noise by the analysis program and 
their sheer number would prevent the program from resolving any fixed- 
phase dependence. 

22. On the basis of the ordinary and partial coherence results, 
it may never be practical to attempt the implementation of a partial 
coherence function with a sufficient number of inputs to adequately 
represent every "source" along the trailing edge. 

23. Much less TE noise is predicted to radiate from the 1.5-foot 
wetted span of the test airfoil as compared with that which radiates 
from the wetted TE span of a full-sized aircraft. For the similar 
flow conditions which exist during landing approach for a Boeing 747 
jumbo jet, the difference in sound power level may be on the order of 
22 dB. 

24. The overall far-field sound pressure level’ of the hump 

characteristic alone with background noise flor removed was seen to 
5.3 


vary as u 
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11. APPENDICES 

11.1 List of Symbols 

* linear separation distance 

* alternating current 

*= analog- to-digital converter 

■ angle of attack 

■> standard deviation coefficient for the coherence 
function confidence interval 

= effective bandwidth = WT * 

= boundary layer 

* background noise floor 

■ bandwidth = T * 

■ speed of sound z 1125 feet per second 

■ common mode rejection ratio 

■ linear separation distance 

* direct current 

* normalized standard error 

* frequency 

* discrete frequency line = n / (Nh) 

■ Nyquist frequency; highest frequency of interest 
analyzed <l h. f 

« sampling frequency or rate 
*= resonant peak frequency 

* lower half-power (-3 dB) frequency 
B upper half-power (-3 dB) frequency 
*= half-power (-3 dB) bandwidth 

* fast Fourier transform 


FT 

G, 

G 


_ G or 
xx w 


ir C 22 


G , 
xy 


or G 


12 


H 

i 

j 

k 


log 

K 

K.// or Kul .// 


I 

L 

m 

M 

M .// or Mic.// 

n 

N 

PSD 

Q 

r 
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* Fourier transform 

* real, one-sided auto-PSD 

■ complex, one-sided cross-PSD 

* time interval between consecutive data samples 

- f _1 ; height 
s 

* transfer function 
“ an index 

*= imaginary unit = /-T 
» a constant 
= base ten logarithm 
= wavelength constant = to/c 

= miniature strain-gauge pressure transducer 
("Kulite" microphone) located in the, acoustic 
near-field position // as per Figures 2.4 and 3.9 

*= length; characteristic turbulence correlation scale 

= wetted span of the airfoil trailing edge 

*= number of averages or ensembles 

*= Mach number = u/c 

* one-half inch B & K type 4133 condenser microphone 
located in the acoustic far-field at position // 

as per Figure 2.4 

*= an index 

= total number of consecutively sampled data points 
in any single sample record 

= power spectral density function 

= quality factor or sharpness of resonance 

= maximum time-lag number; radius; distance from the 
TE to the far-field observer position 
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rms 

rpm 




R 

yy 


s 

s(t) 

S 

S/N 

SPL 



u 



W 


x(t), y(t), z(t) 
X, Y 


*= root-mean-square 

* revolutions per minute 

* Reynolds number 

* real, even auto-correlation function coefficients 
of x(t), y(t) 

= real cross-correlation function coefficient of 
x(t) and y(t) 

» standard deviation 

■ piece-wise smooth continuous function of time 
= Strouhal number 
= signal-to-noise ratio 
*= sound pressure level 

* time 

= period of time; single sample record duration 
= total period of data sampling = mT 
e trailing edge 

* flow velocity (general) 

* convection velocity 

* mean turbulent fluctuation velocity 

* friction velocity 

* freestream velocity 

\ 

«= volume; voltage 

= sample record weighting-function bandwidth factor 
(relative to the bandwidth of a rectangular weight- 
ing function) 

= piece-wise smooth (continuous) functions of time 

- complex-valued two-sided spectral density functions 

= transformation for the distribution of sample 
coherence function; a general variable 


one-sided 100a percentage* point of the normal 
distribution of the data 

angle of attack; level of significance 

power constant; slope 

boundary layer thickness 

displacement thickness 

uncorrected, biased coherence estimate 

ordinary coherence function 

partial coherence function 

kinematic viscosity 

geometric angle; momentum thickness 

3.14159 

normalized cross-correlation coefficient 
mean air density 

lag or delay time; fluid shear stress 
phase angle; momentum thickness 
radian frequency = 2 tt f 
averaged quantity 
complex conjugate of [ ] 


complex quantity 
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11.2 COHERENC 


The two-channel program "COHERENC": 

1) can delay one time series sample record relative to another, one 
data point at a time (with the possible sacrifice of frequency 
resolution as points are deleted). 

2) can be used to find an amplitude scale factor. 


3) 


4) 

5) 

6 ) 

7) 

8 ) 
9) 
10 ) 

11 ) 

12 ) 


13) 


using the scale factor, can amplitude match and then subtract or 
add two time series on a point-by-point basis to yield a composite 
signal. A second composite signal obtained in the same manner 
from two more time series then makes up the pair needed for 
cross-spectral analysis. 


removes the overall mean from each of the two time series and 

2 

amplitude weights each time series with a Hanning (sin ) function, 
computes the spectral values for G , G , G , H , and <J> 


via the FFT. 


xx yy xy xy 


xy 


averages the above quantities by repeating steps 1, 3, 4 and 5. 

computes and bias corrects the ordinary coherence. 

computes the ordinary coherence and confidence interval. 

computes the coherent power output spectrum, G . 

can then scale the PSD quantities against a reference standard to 
obtain calibrated values if a reference and calibration level are 
supplied. 

computes the overall averaged G , G , G , $ , y G in 

^ 6 xx yy xy xy 'xy yy 

chosen frequency bands (method used for finding item 2 above) . 
outputs all above-mentioned averaged quantities on a printout, 
automatically will scale the plot axes and then plot the 

A A A A A ^ * 

G , G , G , <t t y . 

xx* yy xy T xy 'xy 
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TH- OATA OPCOFOING SPITO IN IN/SEC OR CM/SEC. C 

THF OATA R r PRODUC I NG SPEfO IN IN/SFT. OR CM/SEC. C 


If N P T APPLICAOLT, SET P"r = F^P=l .0 C 

c 

crcccr.rrccrrrcccr.ccr.cccc'-cccrccccrrcccccccccccccccccccccccccccccccccccc 


I OAT A — 2 
IPLOTAM 
I PLIiTO = I 


f FC=63.iT 

pro- is.n 


CCFCCCCCCCCCCCCCCCC' 


• N* = TNc-HH r 

N s T HE 

pos^hil'" in? a. 

01 LAYS1 . 


•IAVGS* If THF 
tcr>LPT» THP 


cccc.ccccccccccccccccccccccccccccccc.cccccccccccccccc 

c 

▼ wr MtiMnrr or oata points pff time: pftcoru; i.c.. c 


nr ppfctoal line?, pf esfntl y limited to 512 of c. 

Tine LEAVES ROOM F OF SHI r TINS POINTS l T I ME C 

c 

r 

NOMpro OF AVERAGE 5 PFP FOEOUFNCV LINE. C 

c 

SAMP1 f NO RAT* IN H? . C 

c 


ORIGINAL PAGE IS 
CK POOR QUALITY 
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i ns i 




H 1 r t nr ijfF*ri IN orc.RFrs. 


c 

c c 

rcrrccrccc.ccrr c rcrccrrcC rrrrcc':r cccccrcrcr*'crccrcrccc'-f cccccccccccccccc 
c 

Ns Ti I 2 
lAVGSr fift 

rVTrimo.o 

PSr-n.O 

Hcc.rcccrccccrcccccccccccrcrccrccrccccccccccccrcccccccccrcc ccccr.ccccccccc 


r 

C 

C 

c 

c 

r 

r 

c 

r 

c 

c 

c 

c 

r 

C 

c 

c 

c 

c 

c 


«rr2» . »u r V , AN) 'HKVOT ' ACT. tht Frf.-R r Nrr VALUES. 
M*T» = i s^rriri"* cpsinf*** weighting; = i r.ivrs n cxcao. 

• KMNF | • = 1 C ALOLA T rc TH^ CIINFIPFNCF INT c FVAl; r 0 Onrs NOT. 

• 1 *>L A Y * OfLAYS WH r»'»NN"L5 PY THE GIVEN NUMH*T nr FrjIMTS. 

•ilac.i* o r i. < :y$ ownnh i; •il*'’.?' delays channel 2 . 

*SM<-ap|« is t,;V -l IOING n~LAY rpr • I L AG I • IN MSFC. 

• c Mr aq? • is T Ml ! SI I PING PILAY CPP •ILAC..’’* IN MSFC. 

*KAOP' = 7 H c MUMH7F nr AVERAGE PfP fllr.MATlriM OMM-" 1 . IM'Pr ABF 

• )o ' IO-OIT WUFOS PI' DATA PCINTf PFp CHANNEL * yfFN EACH PUMP 

!*• pivnrn in MAi.r , n yitlos ?oap rata points available °fp 

C*<AMN'-|_. PCCM THIS. r.H- AVFBAGr PEP PUMP IwIYh ?.3AS POINTS PFP 
OMNN^L) MAY n- rrALIFCO, of TWO AVERAGES PCP PUM'» (WITH 1024 

POINTS orr CHANNFl » MAY p- PE *L I F FD » ETC. KAP’ = U324/NI 


c 

r 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

f. 

c 

c 

c 

c 

c 

c 

c 


rrrcrccrccccccrr.r.ccer.r.c'rrrccccsr ccccccccr cr.r.cccrccccrcccccccccccccccccc 

c 

r r I = I • T 

r r i . o 

PPT=SOPT (F=| *pr?l 

M»V> r r ="'OCT(Fr?/pT.l ) 


l*T=i 

FPNF 1= l 

lOLAY^T 
II AG 1 = 0 
II AC.2=0 


C 

SMPAD|=0.0 
SMP AP2=0 .0 
C 

KAOa=l024/N 

c 

ccccccccccccccccccccccccccrccccccccccccccccccccccccccccccccccccccccccccc 


c c 

C • C Al »A* = STL OF THF MIC. CAL. SCUFCE FXPFFSSFO IN OH FUF C 

C CIM r, CM2 I AND .1 r. 4 WHTN APPLICABLE; n THE P A I S E CAL»A = 3. 0). C 

c c 

C TA1 «T» * OP PFAOlNC. prr TM|S PIOGPAM IN THE FRFOUFNCY BANOS OF C 

C TH 1 - c»L. SOtlPCF W H C N MOATA* = I, C 


o o ^ r\ 
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= N r T nlFTF.F '■NC r IN C M U OF I ACTUAL DATA PUN KNUD 
SfTT INf.F )-IC*L . HUN KNUH SFTT INGS I rxPRTARFD IN f)H . 

•MPAS* = I CAI TUI A- r S THF OVHPAI L COM'frNrr AN 'I f.PFCTrU* t PVFLS 

tn up tp u chiirfn Fii.trPTO bands; - ■) onr s hot. 


*M“ NL'IV AND NMIH Filter CARDS C.O H^RC • t XA M'H .'T I PLOWS? NHlHxr>'> 

rrcrrrrccccrc.rcccc'trrcrr.ccrrr rcrccrcFccr.crrcr ccorccccrcrccccccrrcccccf c 


CALI A - o . 1 
CAl I A =0.0 
CAL? A=0 .0 
r At lArO . 0 
CM A A=0. n 

cal in-o.r> 
r \t ri - n . o 
r*\ io=o . n 
CAL AD - o « A 

a A IN 1 = 0.0 
C,A 1N? = 0 . 0 
C.A IN 1 = 0 . 0 
r.At,)A=o.O 


CA| | = 1 0. 0**< { C AH A-C AL I D- SAIN! 1/20.0 > 

CM. ?- I O.OAAM CAL2A-CAL?.1-',AlN2»/?rt.O > 

C«l.. 1 = 10. ■>*•♦( ( CAI.1A-CAL3D-r.A|N1)/20.0 ) 
r A(_4 = I n.0**l i AL A A— r Al. AD-SA ! N A » /?0 .0 t 

NPASr* 

M iw 1 1 1 = 42 
KM! U(l >=104 
NLPV I ? )= 10 A 
NMIHI 'I =?AA 
Nl OW C 1) = ir? 

NH IMI1I-1I ?. 

Ml |1W< « ) =?? 

NMIH (A )=*12 
KLPW|S>=1ST 
Nil I H C 5 I = A F 7 

Tur rnn CWINS CAROS Arr NOT CHANGED! THEY »S C FOP 
CALll'OA T I ON AND 4-CMANNFL OPFRATION. 

tf inATA.NF.,|l ,W TO |0 
I PL (IT 4 = 0 

I PI OTO=0 
Nr S | ? 

PS=0 .1 
in AV=0 

II AS 1=0 
II AS/'rO 


ORIGINAL PAGE IS 
OF POOR QUALITY 


r»r\nrn"»^r»r%n 


rto nn 


rwr^vpi =1.0 

r»*r*rs=o .0 
C*l 10.3 

c*l 20.0 
c«t?= 0 .o 

fHIrl.O 

KAr»r=| 
t.rA<- =5 
ni nwr i > = |-> 

NH?h{ | )=!/■, 

M.nv»i-> > =r •) 

Nll’lll I Ir5« 

cn Tr ?i 
io rru-iHJ - 

ir ( t n AT A • M r • a ) r,n TT 21 
K~ r 1 2 
IPLAVO 

n a 'i i = o 
1 L 1 r. ? = n 

SMT API =« .0 
'VJITPrO.n 

* *rv> = i 

SO Cf1N"|NU r 


COLL TPANS2 


'•rm M.'-.i) (T iri r, | j , i-| , ?r > j 

t»riT* i 3 * r ' 0 n t < t*tl»' « i > . i = i ,?o> 


call irnnw 


CC»’CC''Crrr':C ( TrC»"rCCCCrcCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCrCCCCCCCCCCCCCCCC 

ri QTTPR P'HITINr F1P nr)«rn 


ccccccrrrcccrrrcrrccrcrrcccr.ccccccccccccc.ccccccccccccccr.cccccccccccccccc 

c 

fr I lot ()T A.FO.O. AMO. in nrij.SO.O) c-0 Til o?rc 

C. 

•I MGTHrA.O 

*rrr=3.o 

xFir,HT=< crLPT/’.ojocnrc/ospf 

*!tNC = XPtf.HT/10.0 
V».Nr.TH=*S ..1 
VlMPsIO.O 
KfV-CV=l 
•HXPT^rN 
NITS l*N 
riNOU »=-i 
L I N p *=0 


AAA 


nnonn .•> r» n mn 
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PATA XL A IVT./ « Prr/j* , t IjrNf • . • V •/ 

pATA VL AGL I /• OXK • . • S P r C • . • T F'JM • . • l FV •»•(■!. t • * • PFL A • • • T I VC • # 
r. • pm » • , • • / 

DATA VLAf"_?"»;vv * . • • I * TPl»M» . • l.FV».*Fl « ♦ » • P‘ t A ♦ . • T I Vi= • » 

f, • DM J • . • • / 

CAT.\ VI AHL V7« | r.XV' , • J F"* , •FCTR* . *UM L'.TVL'.' < RE • . *L AT I • • 

t • vr O* . • B I • / 

PA*/. VLAHI H/'PHAS' .'I C-»» ,c '0 • ,«DrCi r »* ,»FF P*.»M4Sf*.» OFF*, 

r. • SFT I* I • • / 

data yl a,hl*' 7 * pr->j • ,• nary • . • con* ••r» fn* • »cr 

PAT A NPrex/r/.N*. If.T.V/l / .MAXCVS/I / • NCUP VS/ 1 / » N*’ TSi’/O/ .Not S3/0/ 

DATA NO* 54/0/ iNAT'.VO/ 

K t j^l ciTA.ro. ot on to to 


n or r a rjRr, rr.p channel » follcw this capo: 

YPAfc=T .3 

vTn r »=2o.o 

or m t = I.N 
xn^n Cl . i i = xr fco< i > 

Y'1^*>( I . 1 t=GXX< I ) 

120 CONTIN'JF 

ir c Yupnl i . n ,L f .J vtnp-s. iM on to no 

YTOn=YTOP*IO .J 
on to i?o 

I to FONT INI/' - 

VM4SF = YTCP-*.n .0 

TF ( YORO < I • I ) aLT • YHA ~T I YOFD ( I .1 ) = YO A “F 
100 CDNTINI/F 

CALI PJCSI Z< t*. 0. 7.0) 

TALI GRAFF IXLIGTH. XL CF T , X 0 1 GM T ,X I NC • NOECX . XL AO EL . XOFO. XOUMM Y, 
f. Y| NO’H, YflAr.r ,YTi)R. YINC ,Nn*CY« YLAOl 1 , YORD, YOUMMY.NSJNGV.MAXCVS. 
T,NCl|FVO .MAXPTC.NOTSl . N’T S ? « NPT S3 • NPT 5A.N'"*TS5»Kt NO »L INFS. TITLE) 


PLOT CAPPS F OR CHANNTL 2 FJILLG* THIS CAPO: 

yrtAorso. o 

V T 0P=20 .0 

nn i’oo ! = i,n 
xnRo 1 1 . t l=xn;roi i ) 

YOPO C I . T t=0YV( I ) 

?»(l CPNT INU r " 

?r fvopn.n ,T). LF.CVTrP-S.il) GO TO 210 
y»nr»=vTr>r»M o. .» 

GO TO 220 

2 r t continue 

YHASF = yTPP-60.0 




nnnnn no n -inn no 


if ivnor>ci.t).i T.vnA viimi i »i i = vpa 5 F 
?3<i ccnt Mur 
c 

CALL rMCMZC I A.O. 7.0) 

r 

call f.FArr i xlnotm. XL” r T , xi ir.MT ,x i*r . np=t x , xi Airi . kopo. xommm v, 

C VI NOTH. VPASF ,VTf«-». VINC ,NO=r V , VI A«l ,1 t Vl)PO . VOU’AM v , flS INGV .'A A XC VS. 
r.N ri iry., M AXP7 f • Nf'TT I ,NPTe.;»,N »7S I.NPTSA.NP7SS.ICIHD.L IN” S. T 1 TLK I 


*>i 'n cappo rcir |c,xv? rru.’iw this cafo: 
vmav--=o. o 

YTnr=-..?o.o 

nr mo I = | ,n 
xoroi i , t ) =* r F"n ( i > 

YM' C' I I . t J-rr.X VAHC I ) 

■»;tn (TN’lNUf 

I r I vn<?nil . I I .ir.l YTnn-r, ,<j | | co Tn 3 | 

VTOPsYTpP. I 0. O 
r.n in 3 ?i 
3io roMTjs'r 

YPAr^nYTOP-AO.O 

I* I Yr iu 0 (I « I ) .1 T »YP AS*" ! YOPDt 1.1 )=VHA5t 

mo Yn ip imup 

CALL "irSIZtl*. 1,7.0) 

• 

call r,p*rrixi noth, xi. fft . xi ioht . x iNC.Nnrr. x, xi aotl . xurra. xdunmv, 
r.YLN •.»« .VII Apr .VTIIP.V | Nr ,NOF<“V ,VL Af)l 3 .VOf i), V 0U N«Y • MS 1 N GV • M A XC VS. 
f.Nnir?vr<.MAXr»TS.N»1M .NP7 52. NPT P3. NPTSA.NO T S'). KINP.t I NFS .TITLE! 


PLOT r»U'DC Ftp. op AS r FOLLOW this cafo: 

vriAs e '=pr.PAsr 
YTnp='»r,TOP 
YINC=AS.O 

c 

OP 7110 1 = 1 .N 

xnupt j , I l-xrsrot T | 

YOPOU . I > = f»HAZP( t ) 

IF ( VITPOt J .1 ) .r.T.VTOPI VORDCI . I »=YTOP 
tr ( vrjp.ot 1 . 1 ) .1 T .Vt *5 r I V'lPOt l .1 IrYBASF 
7110 CONTINUE 

c 

cam. »|C«slZt IA. 3 .T. 3 ) 

c 

CALI r.PAI-r 1 XL NOTH. XL Fr T, XUGHT.XINC .NOFCX.XLAPHL.XriFO, XOUMMV, 

6 Y| NOTH. V BASF .VTfTP, V I NC .NOFC Y . VI. MU • • VORO • V0U4M Y • NS 1 NGV .NA XC YS. 
f.NrunvS,MAyoTF..N'»T«;i ,N?TS2,NPT5:i.NPTS«.NP1SS,Kt NP.L INKS. TITLE! 
c 


ft.ifi fvm i "> r» ,101 
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<o ccniMir 

If < l"»UhT«. r «:.3 ) Ml Til 1*3 

rim cappo pfr-* cph-i -mc* mu o*» this ca'>o: 

vf.A"«--n. o 

Y7P«r | .0 
YIN'"- - !.? 

OP ' 3 » 3 I = | i N 

xn'Jl-t i . ! ) -xrr :oi I ) 

vnr n ()•! »=Ciiicj < i > 

ir { VIIHOU a It .G T .vTi'nj vmol I a n =YTf)r» 
jr ( Yimru ,1 ) .L'.YnAsr » YOnr.< i , j »=Yr> A ?f 
too roNTjHijr 

CALI P!C r 1 Z( M.O.T.O) 

call r.rAfr t XLrjc.TM, xi r c r , *o i ,-.ht.x inc , nopc.x ,<L»nrL. xnr. n, x'>u amy, 
r. Yl NA’il, YHASr a Y Tfin a V I NC a NO”C. Y * VLAOl *>, YOCPa YA| MM Y , U5 I N GV a M A XC VSa 
r.NMjr VlaMAunTAaMil «; I aNPT^p . N"»T S H a NOT S A , N^T <5 5 ,K I N 3 aLINCSa TITLE ) 

O rpNTTNt/P 


CALL °KAI M OaOa'Tafll 


9rjOa fJITP 
CNO 



mi 
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nt»r»Rniir ls r ’ r\N -,3 

TNTrr,-p«? Itir IL . JOSY 

cnMr>L r x rrn .*-'p tp. sxy .*<* y. r,»ynp 

nj m-k«; irs 

01 '*'' nmi'n r.*x->?< si si , -,vvr> ■»< si v j ,<;xy r> 2 < mu 
cp*i**cin n.uvr, 1 ;. im ay, il agi . it ac,/,*- *nf>, iwt 
r n»>trN ktnp J »nmy ck , c rc • r, 'i' ,, ,PS* r r Tpo.r i',n»-;r 
C'immou 3 pi.pt , s« - »i | , ?mi ap-> ,dc | < .nxvnr r 

COIMiTI tsr II. < A vv. I , J'IM< C ?3 I . INK ( 70S I 

r n*MQN ox xt sis) .oyrt isi . ixyaMMsi ,hxvap<«- is I ,mw yr.ro cm 5 ) 
r T'Imon "HA/r <si s>, cnn->cs I s I. CO r r>f:C CO I 31 .C It wr< si 0 ) ,r it»rpj si 51 
crvvrs cmh ?oasi ,r nar’OA si . xFccocoio ) 

COMMO'I f T" I I T :1 SI .r FT?C SI SI .i'.XYI S 1 SI iHX vc f. is I 
r |IM»MN NT4' |IJil|H| | It.MnwtIOl 

T.'IMIIM IOATA ,C At. I.OAL ’.CAL J.CAl A . C AL X . 0 AL Y , C At. XV 
I 10 FnPMATI A 0 A.?) 

is? ’ o r ''4 a t c / / f « * sum’’ rn.'tiA.’ cyjTTro. now iavCiS«**iai 

I'l I l'rt'» T l' 0 ‘’P!NT MIM'irt* ,|A,» HAS A VALUO nr', 11 ?.* ON *J I OH A T I ON DtJM 

t r Ml) . * . t A | 

* 3,1 I ri|TM AT < * | • .////, ox, • I • , 7 K, «<;xx* ,sx, •r,VV».| 6 X, *GXV( PrAl) • , KX, 
r. *r.xv ( ! mac. I • . fx . • I • ,?x,»oxx» ovv • ,i.x . •r.xYcriiAl I • , 3 x , 

O'r.XYl I MAC I *.✓//. CSC I | 0 ,AF 12.1 I I I 


pi - s. t ai s^r? 

#A V*3 

JIVJMPr ■) 

NM'V)K = 0 
IN I - 3 

rrn ?o ir|,M 

rxx c 1 1 =o ,r> 

•tY Y | | I = 0.0 

C.XY ( | I =r M D L X I 3.3.3. n I 
to cohtinh*- 

rn S3 t = i,KAr» 

IMTii i = im 
TN!=-|N?FN + N 
S3 CONTI NUT 


t nr nomppp nr data points SHirTro itimf n r i,*Yri is calculatfo 

•>0 | OU'*P= ! OU *P* | 

| AM Mr API / r Lt) AT C I A VOS I ASPLFT f I 330 .3*FL0A T C JAV I 

I Ar,2=SMr«r2/rLOAT( 1 A V 0 S lA^UT/ I 030 . 0* FLOAT C JA V I 

IL T I - 1 LAG I ALAS1 

II T?=lLAr.?U.A 0 ? 

T*tr DATA IS F” A J (A AX PJOMATION f>UMP> AND CHFCKF3 FOR PRF()PS 
rr-AD (A. 1331 IINT ILC 1 I . t -I .A3CAI 

r 

CD I S3 l*t.A3?A 


fViri nftn ** nm 


ir C INF- IL< I ) .L T . I !’t .ANO. INF It. Cl I . Cr .0 > CO TO 1ST 
Nr»TI’*<= I 

WriT^ C3.I5EI I . INF ILC II • IOUNO 
1*40 f ONT I NOE 
C 

Jr f NOTOK . "O. 1 ) 'in TO ?? 3 

|H", «•= 1 \V, V K 4T" 

ten I to (.Til 1 --?) I O' IMP. I 4 VC * 

NnTC*=T 

C 

c.n TO TO 

c 

continof 

TM r r io^t ■>« i c A'-Ticoro t n chamnfl i; T *-r $rcnnr> to channel ? *'• 

f ( lOATA.ro. 4> CO To 23 E 

on ?3o !=i,;oaa 
CHI ( I ) = I NT IL C I > 

CH->( I > = I Nr J | C r*-?0AA) 

233 continot 
CAIXtCALI 
CA|.V=CM 2 

CALXY^snFT (CALI ACAI.2 ) 

CO TD 2.4E 

TME Tllll.nwlKr. SPECIAL CAP15 APE FOP 4-CHANNEL DATA— 

2JO FONT INUF 

prt 24-> F =1 .1 (04 

CM I t 1 )= I 0().0*< I NT JL( I >*CAL t-INr ILF !♦ I 024 >*CAL2 IFCAL2 
CM2C I I =1 OD.T* C INF IL C F »;>T4R> AC4L3- INF IC< t ♦ 30 721 *CAL4> /C ALA 
240 CON T INO r 
CALX-f Al 2 
CAL Y = CAl A 

CALXY^SOPT<CALX*CALV ) 

TF Nrcr c SAPY. TMF DATA POINTS AFC SHIFTED TO PROVIDE. TIME DELAYS - 

P4*i CONTINOr 

r=-o 

c 

122* CCNTINUF. 

NIK=I 

fval = 3 .n 

K - K 4 • 

jav=javc l 

JM0V = INPt <K> AIDLAY 

c 

no I 23T 1=1 .N 

AC I l=CHI ( I ♦! -!♦ JMflVf II. T | ) 

PdlsCHK I ♦ IaJMOVA II. T I 1 
FVAL*C VAL A AC I ) *Ft I I 
1230 CDNTIN'jr 
C 


nr>o n inn .■» d n n nn 


(TMANN'-t I Mr AN VAU>r CHPP-r.TlnN AND ►ANNIN'-. 
*VAt. r^VAL/CNANI 

nr i ? 3 i I-i.n 

ACM“(A( I l-TVAL ) 

rm«inc h-^vali 
ir.ciwr.Nr.il r,n t» 1211 

A I I ) - A ( I ) *? .OAISINC I | - I ♦*»l/«rLOATlN HI > **2 
r.»N< t • m-* t/ ifloa* cni m 1**2 

|?H rriNTINUr 

fin TO I -> 7 * 

I ? r it> FVA|_=0.1 

rn i-i.k 

AC n =CM,?I I A|- I ♦ JMIIVA TI.T? 1 
n» t »-CH,-»( T 4 | A.IMOVA II T?| 

r val - cr v At ♦ X( 1 1 *r ( i I 

I 2*5 CON - ' 1 NII- 

CHANNF L ? Mr AN VAl.l IT C««f-5CT|PM ANO HANN INfi 
fVAL=rvAl /(NANI 
nn i2Sfi t=i.n 

AH)-I»(II -rVAL * 

*«( I > = ( Ht I l-CVAL I 

ir ( r/T .Nr. i ) r.n t.» 135 *, 

A(l) = A(T>*?.n*<FlA 4 ((!-|)*- , I/( r LOAT(N|))M *2 

P( I »=M C I |*?. 0 *('> 1 N(( t » * n l/(FLnA'’(N» I 1 1 **2 
125 * ri NTlN'ir 


I 2 ^S r A(_L «TT( A.n.N.N.N.t > 
CAl L RFALTP ( A.P.N . I I 


ir cnix.fo.2i go rn iaoi 

g*x, r.yy, ano gxv A»r coM'urcr: and av b ragfo 

no 2001 1 = 1 *11 

Frr| (t|stMax(AU(,nmi 

fixxl 1 1 =fixx( T » *r fai (rcNjurnr i ( I n*rrTi ( i n 
201 1 CONT INI • 17 
f 

NtX=? 
rr tc 125.1 
c. 

Moo nn 2002 i = i.n 

rrT?n i=cmolx< ac 1 », rc 1 1 » 

0VV4 IlciiVVIl )A.F'!ALCCr>jr.Ccrf2Cn >*FF7?C 1 l> 
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GXVf I MC.’xvj ! M01Njr,<rrT ■»( ? ) |*FF1 U II 
21 V C«VJT IMIir 

1 r ALL TH r ’ Pri'r-, T r|t »V''PV'T5 <MAVr,G»> MAV« A C P N TAKEN. O'* 

rr 1 ► 4 '-. r AT A TFOM TH f TATA T.O'JPf r — — — 

*r < jAv.c r . i avg*.) r,n to *000 

ir all t nr •"»n r , r . ihl f avt^agt- <»kad<>*i apt in a stN.Lr ak 

ptination nunn, r.rT another pupa's worth of data —————— 

1 r ( K .<jr .K AGP > r,p TO T.i 

r.n to 1 »?r, 

rm» nr av' or.iw. »irn and CPc c s-snF.rTou« values ate roRRECtro 

FPH FTT PROGRAM rKiPR-,,— — 

A OO ■) CONTINUE 

0tinPM = 4 *N*M* 1 w -,s 
On 4 | .V) ! - 1 , N 
C»*l IPfiXyi! »/nr t 1A’M 

r.vv ( 1 > = r,vv< 1 > /nNU'PN 
OXY< I ) = GXVU ) /DNlir M 
r.K*r>?( I) = (;XX(I >«<CAl **•?) 

RYYPPM I ) =".YY( J )* (C ALT**"* I 
&xvn?( t » -r.KYI I >*|CAl KY*» 2 ) 

4103 C.ONT | NtJF 
C 

FF!T1- <3.43331 I I. GX XPJ>< I I . CV YP2 1 I ) . fiX V" 2< I * , I = I ,N > 

c 

aam prrtjON 

run 


: M$i 


nn r>nr»r*"*o.T n nor* 


FunFrUTiw: ynow 
1NTrr,FP*2 iNr IL • JUN* 

Complex rr T| ,Fr T?,c.xv,H*v 
CCs M «rN NiMvr.s. till \y. ilagi . It iwt 

crm:#rN konf i .not ok , r c .c nr- • rs tup, ns past 
C f’MMPN ,<;«f»r 7 ,iri .pcr.rci ,H>vrrr 

COMMON INF H. I « n-(. ) . JUNK ( 2 1 >. INIM25M 

common r.xxcsi s> .gy ycsi sj .gxy ahcsii. > .HxvAncsi r>> .Hxvomc pi 5 1 

CPMMnN PHA/r I •- 1 5 ) .CllHPI P | 5 I, CO SPEC ( *• | 5 1 . C II WP l PI r . I .C 1 UPH ( P I 5 I 

ccmmun c hi ippopi « rn? i?o< pi ♦ xrrrolPi p » 
rOMHPM FFTI l pi H| .FF’ Pt 51 PI ,GXVC PIP! ,MXY( PI PI 
COMMON NPA5 . N‘l I H C l .1 I ,M PWC I 0 I 

COMMON in AT A ,CAl 1 .CALF.rAL.I .C*L4 .CALX.CAI Y,C ALXY 

c 

i 4 f»i rrnMitcj' .//o' in ■'ht. Fc-or-Nev hand rPMM*,r«.i,» h ? rnsrs.i, 

f.* H7 T Hf nVDAt.L LEVELS ATT AS FCl LO WS 5 • * / / , 

T.* AXXl'iFT.I,' HO, PFISn.f,//, 
r.» r,yy:i ,F7.I ,• OH, F c : • . T«J .0 , //, 
r.» rjxv: • .pr.i , • ->r. p«: • ,fq.*,//« 
r. • cmi-c r NCF: • tFT.ji//, 

t* C PMC P 7 NT SPFfTPUX LEVELS • ,r 7.2 .* OR. PESSrr.bl 
15T7 FppvAT <• 0 •,///, 7X ,» I GXX GYY |GXV| |HXYj 

c nXYOPT PHAPF-HWY COHIERENCF CIUPR CILWR CO-SPECT RON* . 

r ,//// > 

1TP0 fpchay (PX.lA.n.i ,rrn..i ,ri i.i.r ia.a ,f n.t.r s.4,r9.i ) 

r. 

PRI=F C I**? 

BAI = P'.2l A? 

PS 1 **? 

H*YF SOsHXYr T • • ? 

e 

*ixxT = ').n 
C,YYT-=0 ,0 
GXYT=0.0 

c. 

nr. ?oo.-* i = i »n 


THE rtJFCFFNCF, CO-SPrCTRUM ANT PHAS C HETWFFN X t V AP* COMPUTED -- 

f»MIN=( I .O/Nl *».-> 

OPOD=GXX I I )* r.YY I I > 

IF ( ABSC ppno) .tr .PM I Nt PP0D=PPTN 

cfihpc r ism ai. cgxyc i > *ccnjgc gky cm i/ppod 
HXYC I l = GXY( I |/r,xxf l ) 

*!f t cnitrrrNrr to thfc \f eiti-fp gxx of gyv apt lfss than *cutoff*. 

TH’S IS A LC’W^r LIMIT CHETERMINTD PY THE NUMRTP OF A/D OITS ANO 
FPEO l!N:SI SE T TO "FrvrNT COMPUTER CPPOPS rPPM OPCOMING TOO 
GP 7 A T « THE PANT; OF TFF MAGNITUOT OF THF OA T A AT THIS POINT 
LI=S OFTwrCN 0**2 AND 512**2 CFOR 10 -HIT WORDS I . 

CUTnrrsl , 0 /N 

|F cr.xxc I I .LT .CUTOFF .rR.GVYC I |.t T.CUTOrr } COHP<l)r 3.0 

APPLY MIAS COPpFCTir'N ^fl COHERENCE EST I MATE — 



w pQQ|R QUALITY 


^ O A n OAA 
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c 

fcr I .0 

1 r f tWT.ro. I ) w -| .0 

rfiHP< i uniiiPi ii-ii .o-ci)Ht?i i > )/ir.o*w*rtoAi( iAvr.su 
ir icmtpm.Lf.'i.il c.n to isoi 
c 

* i magi mxv< 1 t » 
n CMK = P r At f MX V I T 1 > 

l r ( vi^tncur >.Lr.n't|Ni r)CMx=FviN 

r 

pvmrcu *=F,7.?'3P*A1 AN? ( T ANUP • OCHA I 

T »•=■ PHASF rj^FSTT TOP PLOTTING THr PHASE IS lNTPOrMIC r t> ----- 
PSTn”r'v; ♦ | TO . O 

psp A«; r.=pr -i oa .0 

ir 1 PH A 7 r C I ) « GT . os TOP I PHAZT < I »=PHA2S» I 1-3N0.0 
ir CPMAzri n .lt .p'-.nA?'-| puazf-c i i=pha xfc i ifsao. o 

GO TO 1 SO I 

I soo COMP Cl | = 0.0 
PH*7'-( I 1 = 0 .7 

»r < i .ro.i ) c.o to isoi 

PN\?rt 1 )=PH\7F< I-l ) 

i so i crttiT inuc 

CIU^P ( ! 1 - 0 . 7 
r ILWCC T 1=0.1 

ir CKONn.Fo.oi c.o t o ?ooo 

o^t^p^inf th~ cm-rr' c Nrr pit confiofncf inttsval 

7 As I .or, 

AVC.CsT AVC.S 
0=7 .*57< Avr.s- 1.3 1 
CPMr=COHR( I 1 
CnM=SO->TCCOHC) 

VAP I= r .or TCOI*CI..7-'3.374**<l.6*CPHf»3.?21) 

7= 0 • N* ALCIG I ( 1 . 0»CPH) /< | . 0-COHI ) 

rMH' = TANMt 7 — 0—7 A*V AB T ) 

tr CCOlN.l'-.O. 71 C M lN = *7.7 

CMtN=CMIM*C'7IN 

r^AX = C T7NHC 7 ->♦ 7.A*V AT I >)**2 

ClUOPC I 1 =*C M A X 

C1L *F I I 1--CMIN 

Till- TPANSrro FUNCTION, -|FT CFCSF SP r CT BUM AND X r. V CHANNFL 
SP r C T PUM? , A NO OTHFP VALUES AB* CCMPUTFO --- 

C.TTt l )=OXX( l >*( CAt X* ♦.?) 
svvi i ) = gvy C I) * ( CM v* *p > 
r.xvi 1 >=c-yvi I >* (cal xv**?> 

► XV | | )=HXYC 1 )*{ CA< V7C*l X 1 

c 


nno r% n nn-i 
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rrr.n'C c n r ow C l > «cnt-r ( i j 

|.vy(|pr,( ] |;(jv»| 1 )/,,> x( 1 J 

r.»y API i »=?cAi,r< <-.xy( i » » 

I )=r*HC(MXYj I I I 
p.XXT=i;*XT»>-.*X< l i 
<'.VVT-'-.YYT«r,YY( t > 
l»VTrr,xrT»r,|(YMII 1 ) 

y hp rr^cr'-i’nNJtNr, "anh '•=nitfg fFrorwc ap- c*'MPiiTr.'» 

Xl=l-I 

xfbpom J = ( xi i r*r -r i ro» 


20.13 CONTI NUT 


tr (NPts.i r.nt ■ii t tii ?*nn 

TM” rv-rxiL cnHrrrNC-. 'tc., **nr cnr.'s^N rands is COMF'lt^n ------ 

rn ?5’.i l=i ,N'H- 

c 

r,*¥TM=T .0 
<»vy rn= t • s 
gxvtr-o. o 
C 

KI. = NLD*( 1 ) 

KH = N»»IM< I | 

r. 

no py nn j=k». .kh 
r.XXTn= r .KXTn*GXX( J ) 
r,YYTP= -,VY' , «»PYY ( J t 

IK VT|lT-r,xV T n*';KY ADC J> 

2 r ^1 PMUTIHU" 

c 

CI'TOTO= I G» YTH*r,XYT 1 I /( PrX KTn*GVYT 0 I 

mon rm =r ot nr n * g y y T n 

c 

r,XKTO= I 0 . 0 YALP 1 I 01 '-,XKTp/C^| ) 
r.YVTfr I I.OAAt.nOI 01 GYY T R/P 32 I 

c,xvto=i o.n*ALnr,i o< gxyip/r.yi 

cnt »5 o”= n . o * alo r,j 3 c c /p s? i 

WFItr ( I.TtrOI xrf.FnlKL I . KFPEOIKH) .GXXTO .PEI .GYYTH.PE 2 .GXVTR.RE 3 . 
GCOTOTR , '-OPSDR. Pl?r> 
c 

2500 CIINTINU'* 

c 

2400 CCINTINUn 

c 

COTOT - ( vlXYT* GXYT 1 / < OXXT*GYYT ) 

C.ppr, p = c OTOT • r, Y V T 

c 

f.xxT*io.)»ALnr,nir.xxT/Fii I 
riYVT* I O • 0 A A 1 f Gl HOWT/PT?) 

GXVT=l 3 .J*AlllGlvMGXY 1 /P'-. 3 l 
roPSP=i o.o*ai nn 01 consiVFS?l 
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c 

r tmp n»ri , p«rr Ann imrr vu.u£R rnF cach hani 17 

C 

ni7Ks:rimTF*f m x* 
riiTv=rvi T iiT*r »i v**^ 
f 'll x. v = cu TnrF al xy • 
r 

HI Trt.'O lrr|,N 

c 

ir (GXXI M.LT.rt'TX) -.XX ( l * =Cl'T X 
* f (GYY ( ( » .L T .niTV * GYY(t)=CUTY 
II ( P.XYAIM U .LT.CU7XY I GXYAP ( U = rUT X Y 

ir < rfi r .."<rr ( n .1 t .nit y*i r our ( i » ♦ r . 1 7 1 > ) cn«5nfcil> = 

r.njiy t cnur ( 1 > t 7 . 70 I I 

t r (MxvAnii).Li, r 'ir’ i oiiY/pumi m x v ao i i * = fopt (city/ cut x * 
ir tHXYlll'OI 1 » .1 T.l fUTY/ri»TX II MXYOROI I >=(CU1 V/CUTX* 

c 

r.xv ( j )= i o.o* M (ir. i oi r,xx( ! »/psi * 

<-.yy i u -i o .. 7 * vi. nr, 1 3 < cvyi t )/Rs?? * 

o x v a n ( i ) = i o. 0 * ai or. i o< r. xyaou >/r *.7 * 

rni;pfr( ( >r|7 .7* ALIK. I 3(C')F-*FC( H/P6?> 

HXV AM< ! » -| .7.0*- ALHOI 0<NXY *17 C I IXHXYFSO » 

HXYUPOI | * = Al IIGI (7<MXVnrn< | )/(HXYF e 0**2> » 

r 

7<,oo coutinu 1- 

c 

c 

HXVAHt 1 »=o.n 
Hxvnrni i »=o.o 
PMAZM I ) =0.0 

c 

Wl-Iir (.7.3*671 XFPF 0 ( I ). X c RIT 0 ( N ) . T.X X 7 » R F I « G Y YT . RT «? . G X Y T . P£ 3 • 
f.rPTIIT.CIJPRO.P- ? 

WO ITF (7,76 70* 

WFITT (3,37,7(71 ( ( , XT RFO ( 1 I . GXX ( I ) , O.Y Y ( I > • GX Y AB ( 1 ) . HX YAH ( 1 * , 

r.MX YPC7 ( J * ,0(1 A,’T( I > ,rriHR( | ) ,ciur>r m ,C1 LWP( | * .COSPFC ( I * , I = t ,N» 

c 

CCM RFT'ITN 
FN0 
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11.3 PARTIAL 

The four-channel program ’’PARTIAL": 

1 ) removes the overall mean from each of the four time series sample 
records and amplitude weights each time series with a Hanning 

2 

(sin ) function. 

2 ) computes the spectral values for G^, G ^ , G 33 , G yy> ^y 2 ’ 

Gy 3 , G j2 , G 13> and G^, via the FFT . 

3) averages the above quantities by repeating steps 1 and 2 . 

4) computes the ordinary coherences. 

5) computes the partial coherences. 

6 ) can then scale the PSD quantities against a reference standard to 
obtain calibrated values if a reference and calibration level are 
supplied typically from the output of COHERENC, step 11). 

7) outputs all above-mentioned averaged quantities on a printout. 

8 ) automatically will scale the plot axes and then plot the 

p, p p -2-2-2 -2 -2 -2 

yl’ y2* y3* Y 12* Y 13’ Y 23* Y ly*23’ Y 2y-13’ Y 3y-12 * 


c 


lNTTCrptp INF |l 

rr>ifirx RFM .crT^.rrn.rFT* 

COMPLEX *YV, ?YR, SPY, fpo.S YNP. lYunV, GKYP.PCnM.G »VD 
C PMr>LF <t r.vv .01 1 .r.p-» « C-x.T , GY I ,GY? ,GV3 , Gl ?. G I 3, G? 3, GAG 
C 1 M r N 5 ’UN INriLIADCf. ),FFFD(5|5) 

r 1M«-NS | UN CHI < I 024 ) ,CH2< I 3 24 ) .CH3I 10 24 ) , CH4< 1024) 

M MEN", III N r,YYC515),GlM5l5).G;*.-M5l5)«G33<5l5)«r,Yl<'. tS),RY?I5t5» 
m M=hG ION r,Y3 (515) . G l?( r 15) . GI3< 5 I « I , G2 M 5 15 I .G*0< W ) , SV Yl 4 ) 
DI ur NG ION FFTI 1515 ) , FFT 2<5 I 51 ,FFT3( 5151 .FT T 4| 31 5) 

DIM^N? I UN 5VPI 0 ) ,5r>Yl 4 ) , SPP( 41 ,SYNP( 4 I »5YNPY<4 ),A(5I5).P(5I5> 
CI<t r N5ION GX YPI 4 I , PCGHI 3 ,5 1 5 I 

D1 M^NSirN XnRD< I , r I 5) , YORDI I .SIS I ,X DUMMY! 1.5 15) . V DUMMY <1 ,5151 
OIMrN5 ir N VL Al’F 2 | # * ■) ,YI. ARE3<23 ) 

O l W NS ION KlNOt 5) , M AOF|. | 30 ) , YL ADEL ( ? n ) , T I TL I < 20 ) . T I TL2 ( 2 0 ) 
TIMTUS irn TrTt.3<?,)>.TJTL4(?3),TITL5( J 3 ) , T I TL M 2 1 ) 


152 cnr-MAT ( ///, • OU'fF MO ,’,14,’ CEITTFD. N F W I4VGS-’.I4) 

I 55 FPC'UT ( • •’"•DINT NUMf» c F * . I 5. • HAS 4 VALUE PF • , II 2 , • Of, M I OM AT t ON HUM 
IP NO.’ ,14) 

?40 FOPM A T « 20A4 ) 

||.?1 rt'PM AT (4 P A3 ) 


I 7 IF* r PPYA” < • 1 • , ////,5X , • I • , | 5* . • GYV * , icx, « 31 l • , IPX , • 02,2 • , 19X . 

r. • G3 3 • . I ?X , ’GY I • . /✓, I ?X, ‘P~AL • ,?X. • IMAG* , 7 4 , • RF AL • ,7X. * I MAG* ,’X , 

T, ‘TR AL • ,TX , • 1 V45‘ ,7* , » FE AL* . 7X, • 1‘MAG* ,7X , ‘REAL* . 7K, • 1 MAG» , // ) 
irr rprM AT ( *1 ’,////, FX, ’I • ,iriX. , GY2 , ,l9X, , GY3' . I 9 X , • Gl 2 * , I 9 X , 

r. ‘Gt 3’ .1 °X , ‘C.2 I* .// . I 7X, »R" AL • ,7X , « I MAS* , 7X. 'RE AL • , 7X. • )MA5* ,7X . 

C 'T AL • . 7X , • l MAG* .-x. T’ AL‘ , 7X. M MAG* ,7X , • RFAL‘ . 7X, • IMAG‘ .// I 
) 375 r P PM AT 121X, I4,2F?3.?1 

?A50 rOTMAT ( /// ,’ CMINV rAIl. r 0 AT FFFQ. CDMT>.‘,I4,‘ of fhannel • . I 2, • 
r.CI=',FA.i,' 02=’ ,ro,i ,//) 

A 04 Y rnOM AT ( • 1 • ,////. RX.’MAGN I TU3T HF CF PS5- SOfC TP UM OF T WEEN CHANNEL | 

F, AMP THF miTP'l’ CHANNEL A i • , // , » X . 2 0/ 4 , / // , 2 3 X . • ! • , l 3X , • FRF OUFNC Y • 

G , 1 AX . • X-SPFC ’ ,/ ' ) 

404 1 FORMAT f • I • ,////. PX, ‘MAGNITUDE OF CFOS 5 - 5PFCTRUM HFTNECN CHANNEL 2 

r. AND Thf OUTPU T CHANNEL 41 • »//.PX. 20A4 , /// ,?3X, • 1 • , 1 3X,‘FRE0UENCY* 
r. ,1 4 X. » X-SPFC* . //) 

4042 FpPMAT c • | ♦,////, ox, ’MA5NJ TUDF CF CFOSS- SPECTRUM REThEEN CHANNEL 3 
r. AND THF OUTPUT CHANNEL A Z • , // • P X , 20 A 4 , / / / , ? 3X . • I • . I 2X , • FREOUENC Y • 
F. , I 4X ,* X-SPFC.* , // ) 

4053 rpPMAT | | ////.OX. •ORDINARY CCHTRENCE BCTWCEN CHANNELS t AND 21* 

f. . // ,5X . 2 0A4 ./// «?3X , • I • . | 3X . ‘FREOUENCY • . 12 X, • COHERENCE • ./ / > 

4051 rOTMAT I • 1 *,////, 5X, ‘ORDINARY CPHFP r NCE DETWEFN CHANN=LS I AND 3:* 
r. ,//, «X. 2 3 44, ///,23X. ‘ I * . t IX. ‘FRFOUENCY * . 12X. » C OHF PENCE • , //) 

405 2 FPPMA- (• I*. ////,8K. ‘ORDINARY CCHERENCF EF.TVEEN CHANNELS 2 AMD 3:' 

F..//.9X.2DA A . / // ,23X . • I ■ , l 3 X , ‘FRFGUrNCY * , 12X , • C OHF RENCE • , //) 

ADft!) F nr mat f • 1 *.//// .5X , ‘PART I AL COHERENCE 0 TTA EEN C HANN C L I AND THE 0 
CUTOUT CMANNFL 4 I • . //.AX.20A4 ./// ,23X , • I • . I 3X .• FRFOUENCY* . 
r. 12X . • COHEF F.N cr • ,// ) 

40M TORMAT l‘l‘ .////.AX, ‘PARTIAL CCHERFNCF OETVFF.N CHANNEL 2 ANO THE I) 
CUTOUT CHANNEL 4 t ‘ . // . P X . 2 3 A 4 . /// . ?3 X , ■ I • , I 3 X , • F R EDUF NC V • , 

F. I ?X , ‘COHERENCE* .//) 

A OF 2 format I • | •,////, 9X. ‘"ART! AL COHERENCE RTT KEEN CHANNEL J ANC THE O 
f.UTPUT CL- ANNFL 41 • .//.AX, 23A4 ,///,23X .• I • , MX. •FFEfYurNCV* , 
f. I ? X , ' CllHEPf: NCF * .//) 

A 2 T I FORMAT (AX. 14, I OF t | . | ) 
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ccrrrr rccrcc.ccr.c':r.crcrr.r.rrcrc'‘crcccrcccccccccc.cccc<.r.ccrr.ccrr.cr.r.ccccccccr, 


r. c 

c *toi.ota* i plots Thr MAGNiTUnr nr THr cross- spectpum hetween c 

C INPUT CHANNFLS I,?, r. T V.S. THT CUT PUT CHANNEL a; = O DOES NOT. c 

c c 

C * I n L CT B • = I PLOTS CRO INARY COHERENCE BCTWEEN INPUT CHANNELS C 

C I r. 2. | £ 3. AND 2 C 3; = 0 OOFS NOT • C 

c c 

C •|r*LCTr.' •= I PLOTS TM; O APT I AL '"OHEF ENCT: BETWEEN INPUT CHANNEL S C 

C 1.2. f. 3 V.S. TTT OUTPUT CHANNEL A *, = 0 COER NOT. C 

c c 

C ‘PEC* IS THE DATA fTCOPrilNG SPEFD IN IN/SEC OP CM/SPC. C 

C »rEP* IS THr OATA PpPPOPUCING SPFFp IN IN/SFC OR CM/SEC. C 

C TI NOT APPLI CABLE. SET RFC-PF.P=I.O C 

C C 


ccrrcccccrcccccr.cccccr.ccccccccccccccccccccccccccct.ccccccccccccccccccr.crc 

c 

I PL 0 T A “ I 
I °t HTS-I 
|Pt.OTC=l 

c 

RFC=NT .0 
PFR= I s. 0 

c 

ccrccccccrcccccccccccccccccccccc.cccccccccccccccccccc.ccc'-rccccccccccccccc 


c c 

c *N* = ONE-HALF THE NUMBER Or CATA POINTS PEP TJHF PFCORQI I.E.. C 

c n * the NU .in c P nr eptctral lines. c 

c c 

C • t AVGS* IS THE NUMPER OF AVERAGES PEP FREOUENCV L INF . C 

c c 

C * RPL PT • IS THE SAMPLING RATE IN HZ. C 

c c 

r »pEF* IS THF CEFFRrNCE VALUE FOP THF DFCIOEL SCALE. C 

c c 

c •RAPP* = THr NtJMHCr CF. AVEPAGES PEP BIOMAT ION DUMP. THERE APE c 

C AO-jr IT-BIT WOTOS op DATA POINTS PEP CHANNEL. AiHrN EACH OUMP C 

C IS Pivipro INTO FOURTHS. IT YtFLCS I 02 A DATA POINTS AVAILABLE PrPC 

c channel, fpom this, one avfragf ff.r. dump iwith 102 a points rep c 

C CMANN r L ) MAY OF PT AL I ZEO . OR TWO AVFPAGES PrR DUMP I W I TH 512 C 

C POINTS ppp CHANNFLI MAY BE PEALIZ r D, ETC. KACP * USI2/NI C 

C C 


crccccccccccccccccrcccccccccccccccccccccccccccccccccccccccccccccccccccr.c 

c 

NFS I 2 

|AVGS=T? 

eOLr.T-=S 000.0 

P.FF s I . 0 
r 


KAOPsSia/N 


cccccccccccccccccccccrccccrr.ccccccccccrccccccccccccccccccccccccccccccccc 
c c 

C • C AL •* * = SPL OF THF MIC . CAL. SGUPCF EXPRESSED IN OB FOP C 
C CHANNELS I.2.3. C A. C 
C c 


S' 1 Fm 

■ u8 


inn rt <o p n >o .'innecnno 


146 


• cai*m» = on nr/vpiNr, orr th" c*i ippation pfugoam in tit 

FP r OU 1 " NC V PAND3 fir THF CAL. SCUrCE t*HFN • I 0 A T A • = I. 


•GAIN** = NFT DITFEPTNCF IN GAIN Or (ACTUAL DATA PUN XNOH 
5 FTT FNGG )-( C AL • PUN PNOn SETTINGS! F*PP' 5 S e [) IN ON. 


ccrcccccrcccccccc.cccrcccccrcccccccccrccccccccccccccccccccccccccccccccec 


cal i a=o. n 

r AL ?A = I I 2.S 
CAL3A-I I 2 • R 
CAL4A= 124*0 

C*LIB=0.0 
CAL ?R-47.o 
C At 30-00 . 7 
CAL 4P=4 ( .n 

€* INI =0.0 
04 IN2 = 3 C .0 
CA !N3=3 r > .0 
04 I N4=3S . 0 


FAHelO .0»*( (CALI A-C A| I " -G A IN I I / 30 . 0 ! 

CAI. 2=1 0. 0** ( ICAL2A-CAL2H-GA I N2 1/20.0 ) 

CA|. 3=1 0 • 3 * * ( ( CAL 3 A — C *1. 3P— G A J N3I/20.0! 

C AL 4=10. 0**C ( CAL4A-CAL 4B-G AINA 1/20.0 > 

fin l?3) 1 =i.k 

GVY( I) =CN=LX<0.0.0.0 I 
Gl I ( l l=cnr>l XI 0. 0, 0.0 J 
r ,?2 ( | ) =CMraLX (O.O.n .0 * 

G33( I l=CM n LX(0«0.0.0) 

GY I ( T I =CMOLX ( 3 .3. 0.0 ) 

C,V?( I )=CNOLX(0. 0.0.0! * 

GV 3( I )=CMOLX( 0.0. 0.0 I 
C I 2 ( I 1=CM°LX (0 .0,0 .0 I 
G I 3( I 1=C MOL XI 0.0, 0.0! 

0? 3 ( l)=CNPLX( J.3.3.0! 

1200 CONTI NU~ 

P I = 3 . 1 A I SP27 
IPUMP-O 

notok=o 

ir>»s ipump= louppfi 

K= 0 

TH ,r OAT* is READ ( A 4K 0 I Cl MAT 1 ON OUNP! AND CMFCKFD FOR FPCOPS.- 

PF AO (4,11001 ( INF 1LI 1 ). 1= 1 .4 J9M 
C 

NO I GO I = 1 ,409 0 

IF ( |NF|L( 1I.LT.I024.AN0.INF IL( t I.GE.OI GOTO ISO 
NOTDK= 1 

WFITF (3.ISS! I . INFIL( I i.ioijmc 
IS' roNTiNur 






no no o o n oo non 


C 

ir (NnTriK.ro.O) tin To 12 3 

I »vr,?= i avgs - 1 

|Ol»MP=l DUMP- I 

VPITF t 3 • I <5? > IOUMP, 1»V« 

NOTCK=0 

r 

GO TO » r . I*> 

C 

123 CONTINUE 
J*riv = T 

iei#> crN*iNo= 

r.rr * i 

TMr r)r S T IK tr. ^r.niCN'-D TO CHI, *hE 5FCCN0 TO CH2, ETC, 

m ll i = i , i T 2 < 
rn|| I tsIKfU. | 1 )*CAL1 
<*H’ 1 1 1 = I Nr 1L I I * I 05>4 1 *t A) 2 
CHTI I I =lNr It II ♦‘‘OAOl *r AL 3 
CMAJ I ) = Itir ILI 1* T072I *C ALA 
11 CONTlMlir 


C V Al = 1 • 0 

no pp i = i»n 

All »=CMI 1 1*1* JPCV- 1 I 
D1 I I =CH 1 ( t*I*JMnV) 
eVA» =rv A1_»A1 1 1*01 t > 

?2 CONTI Ml IP. 

EVAL=EVAL/FLCAT1N*N) 

150 33 1=1 , N 

At 1 ) = ( At I l-CVAL) *2 ,0*1 SINK I -l »*P1/ 1 FLOAT IN) ) ) ) **2 

nil 1=101 1 1 -rvAL )*r,o*isiNi i i*f i > /ifloatcm i i i**2 
3 3 CONTI N'JF 

C»Ll FFT 1 A.n.H.N.H, I 1 
CAH OEALTPI A ,0.N. 1 I 

OO AA 1=1, N 

Ff T| ( 1 1 =CMPLX1 All), 0(11) 

Gill Ha'll 1 1 1 lACCNJOlFFTIl I l|«AFTim 
• A FONT 1 NO'" 


EVAl = 0.3 

c 

CO *>5 1--I.N 
All I =012 I 1*1 * JMC V* I ) 
Ol T )*CH?( 1 *1 ♦JHOVl 
CV Al *CVAL*A1 1)*F1 t ) 
«S CONTINUE 
C 


on n n no n or* 


FV*l -oval /float c N*N1 

nn /•/ i = i,n 

AC I > = ( AC 1 I-FVAL 1*?.0*< MNC C I - 1 1 *P I / C FL OA T ( N ) })) **2 
PC I >=ID< I l-rVAL»*?.n*(SIN( C 1*01 1/CFLOAT (N) } 1 |«»2 
*>/. continue 

CALI. rFTCA,n,N.N,N,l I 
CALL R r Al 7RC A.O.N. II 

CHI 77 1=1, N 

rrrzt i i=rMPi xi ac i i ,i>c n i 

C22C l \-r,Z?A\ iFCCNJGCFFTIMIl 1#FFT2C 1 I 
▼7 CONTINUF 


FVAL=9.0 

no no t= i ,n 

AC I 1 = CH.1C 1 4 I 4 JPOV- I ) 

PC 1 >=CH7 < 1 ♦ JMCVl 
F VAL=F VAL ♦ AC 1 1*BC I 1 
nn CONTINUE 

e V AL=E V AL /FLC A T C N * N 1 

DP 99 1=1 ,N 

AC t 1 = C A( | ) -rv AL 1*2 ,0*C SINC C I - 1 1 *P 1 / C FLOAT C N1 ) 1 1#*2 
BC I l=CO< I 1-F VAL 1*2.0 * C S 1 N ( I 1 • D 1 l/< FL O AT 1N11|1#*2 
99 CONTINUE 

C*LL FTTCA.n.N.N.N.l 1 
CALL Rr ALTPC A.B.N. 1 1 

DO 1010 1=1, N 

rrT3C i i=cmplxc ac 1 1 ,nc 1 1 1 

G33C 1 1=G 33 C IJ*CLNJGCFFT3CI) )*FFT3C 1 1 
1019 CONTINUE 


FVAL=0.0 

C 

no i ii i 1 = 1 ,n 
AC I 1 = CH*C I * I ♦ JMO V— I 1 
PCI 1=CH4< 1*|4 JM0V1 
EVAL= C VAL + AC 11*0111 
1111 CONTINUE 

c 

FVAL*FVAL/FLCpATC NfM 

c 

DO 1212 1=1, N 

ACt >*CAC 1 1— EVAL 1*2 ,0*C SINC C 1-1 1 * “1 /C FLnAT CN 111 1 *•? 
DC I l=Cni I1-FVAL1*2.0*CSINC C I*PI 1/ CFLOAT C N 1 1 » 1 * *2 
1212 CONTINUE 
C 

CALL FFTC A.P»N,N.N, I 1 


DC 


fpIv/lt'F 


mn o -1 o nn 
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CAUL P^ALTPJ A.H.N.I ) 

c 

OC 1313 f = I.N 

FfT 4 ( I >=CMO| MMtl.lMtll 

err < i > = gvvc I >*ccNjr.(rr ta ( I > )*fft At () 
1313 CONTINUE 


CO I A I A 1 = 1, N 

QVI ( ! IsGVIt! ItrrTM J )*CCNJG (FFT4 ( I 
GY2 ( J) = GVrC I ) *F FT2 ( 1 ) A CON J G t F F 7 A ( I 
GY3( I » = GV3f I )»rFT3 ( I )*COKJG(FFTA (I 
CI2C 1»=G12( I )*rFT2( I )*CPNJG(FFTI ( I 
C| Mil =GI 3(1 >*r-rTS( I ) *CCNJG IFF 3 l (1 
G?3( 1 l = r,?l(l I *f FT 3 1 I )*CON JG(FFT? ( I 
I 41 A CONTINUF 


IF ( fXIMP.Fn.I AVG5 I GO TO |7|7 
IF (K..~O.KAP'»l GO TC 1515 
ji'nvsjvnv*N»N 
GO Tp 1 6 If. 

I 7 l 7 CONTINUE 


FNP OF AVFBAGTNG. AUTO AND CROS S- SPECTRUM V AIDES ARC CORRECTED 
WITH FFT PROGRAM FACTORS.-— — ------- — — 

ONor m=a*n* iav,s*n 

DO 4|0O 1=1, N 

C.YYC 1 )=GVY( 1 I fONOPM 
Gill I ) = G I 1< t l/ONPRM 
02? ( I l = G?2( I I 'DHOPM 
C>31( I ll ~GT3 ( I I A) NO CM 
rr 1C 11= GY 1(1 l/DNORM 
GY’ C I )-=GY2C I l/ONCFM 
GY 31 I ) = GYM I ) AJNCPM 
Gl 2 ( I ) = Gl .’( I l/ONORM 
0) 31 n=G13< 1 l/ONORM 
623( I )=G23( I l/DNORM 
A I 03 CONTINUE 

WPITE (3,121(1 

wrirr (3,4201) C I • GVYC I » ,GI I C I »,G22C 1 1 ,G3?( I I, GY 1 C I | , 1 = | ,N) 
WFITE (1,1217) 

WFITF (.1,4201 1 ( l. GY2CI1 ,GV3( I >.G12( t > . G ) 3( I ) . G 231 I > . I » 1 *N ) 

PA RAM FT FRF FOR PLOTTING ARF INPUTFO — — -- — — — — — 

PF AO ( 1*200) (TIT(.KI). 1=1*20) 

PFAD (1 .200) (TITL2( 11,1=1, 20 ) 

PrAD (1.200) ( T tTLT( I ) *1*1 ,20) 

READ (1*230) (TITLA(I), 1=1*20) 

•TAD (l,?W) (TITL5(II.I*I*20) 


o«.“* rt n ft rtrtnnrt 


*P*n lt.2?01 ( T I tlo (I)«l = l» 20 I 

C 

XLNr.TH=6.3 

xlf:ft = o.o 

XF ICHT=( cni.rT/2 ,01 •( nrc/fi£P| 

X I NP = XP T GMT/ I 0 .3 
VI vir,TM = *-,.o 
v iNc-n. 

Nr>FCV=l 
MAX' 1 ! S=N 
KPT*; | = N 
«IND( 1 |ew| 

LIN=S=0 

C 

OATA XL AOTL / • T B r 0 * •' U r NT • • • V IH*.*7.1 •/ 

OATA VLAn“LX*r»nwr« , • r> S P • , • *?CTP • , • AL B* . »FNR I * . • TV 

PAT* VI M'CJVPAKT* •* 1 AL * • • COP*; • • • PFSC * • • E •/ ” 

FATA VL VPPV/ 'l>nni • , *NAPV •. • COH • . *IFB EU* . *CF •/ 

c 

OATA NOFCX/O/.N'MNOV/l / .MAXCVS/l / , NCUD V 5 / I / . NPT 52 / 0 / , NPT S 3 / 0 / 
OATA NPTSA/.T/.NFTS^/O/ 


pi ot capcs pop |crif FOLLO* THIF CAPO: 

VPA5F=0.0 

VTOP=?0.0 

on til i = i .n 

x i = i - 1 

rPPOl I »=l*Pl PT*PPC*XI |/( 2. 34PFP4FLOATIN > ) 

xrFpi i . i i=rm< ! ) 

VnPOl 1 . 1 I = , )0P T I PF AL (CChJGIGYtl 1 I I *GV 1(11)1 
vne nil , I 1 = 11.* ALOCI 01 VCBDI 1 . I l/FEF I 
ipo crNTiNu*" 

IF I VOPOI I I M.LF.I YTOP-5 •111 GO TO l|0 

VTO’ - '=VTnoT to. o 

60 TO 120 
»l« CONTINIir 

VnASE=YTOP-ftO.O 

IF (YriRPIl . I ) .LT.YOASEI YOROI l . I >=YRASE 
111 CON T I NUf 

XP1TF (3.40401 IT!TLl(ll,l=l,20} 

WPITF I3.107M I 1 , XPROl I . I I .VOPOII .1 1.1 = 1 .M 

IP I JPLUT A . r 0. 0 1 GO TO 10 

CALL P1C61 7. ( I 4.0. 7. 01 

CALL GRAFF I XLNOTH. XL FF T » XR IGFT « X INC .NOECX . XLARFL . XllRO. XDUHNV. VLNG 
VTM.vnASC .VTOP.VlNC .NPECV .VLAREL.VORn.VDURMV, NR INGV.M AXCVS.NCURVS.M 
TAXPTS.NOTSI ,N'»TS2.NPTS.V.NPTS4.NOTS5.KtN0,L!NER,TITLI » 

10 CIINTINUF 


nonAn f% n n nn 
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pi. ot rA» r ) , 5 fn>! |r.y?| follow this cafc: 

VHASEO.O 
YT0P*-»7. ? 

rn ? 2.? 1 - l • N 

xopdi i ,i )= ffeo< i ) 

’ YORDt I , I ) = '50PT(BrAL(eC6JG(GV2( I > >*GV2( t > ) > 

VOPD ( I , I >=IO.*ALCC,IO(VOKOI I , I )/PFF) 

220 CONTI NUF 

IF ( vnuri I . 1 ) .LE .(VTOP-5.3) » GO TO 210 
VTOB=YTnn»IO.O 

on to 2 .?o 

21 o CONTINUE 

vr*sr=YTOp-Ao.o 

|r (YnP0( J . I I .LT.YDASC > VnPDt I. I > = Y0A5F 
222 CONTINUE 

WF1TF (.1,4041) (T I TC ?< I)»I=I,20) 

WRITE (.1,10761 ( I.XOPOI I . ! > .YOFOd . I 1,1*1 ,Nl 

ir f IPLOTA.EC.1 > GO TO 20 
CALL PITS! ( 14.0,7.0) 

CALL 'IR AFF ( KLNGTH.Kt.2F T . XR I GMT , X I NC .NDFCX . XL ABEI. . XORO , XOUMHY.VCKG 
TTH, YU USE .YTOP.YINC, . NOrCV . YL AHFL . VOP 0 . VOUMM Y . NR INGV .M AXC VS.NCURVS »M 
7AKOT S.NPTSI , N^T S2. NPTS3 » NOT 5 A .NPTSf, • K I NO .LtNES.T ITL2 I 
20 CONTINUE 


FLOT CAPOS F C.P | GY 1 1 FOLLOW THIS CAPO* 

VOA^-O.O 

YTOP=2«.0 

or 333 1^1 .N 
XOPDd . I l=FRE0( I ) 

VORD ( I , I )-=SOPT(REAL(CCN JGIGY3I I ) >4GY3( I) ) I 
YOPO ( 1 , I )=I0.*ALCGI 0( YCRO( I .1 l/PEFI 
320 CONTINUE 

IF ( VORDC 1 .1 > .LG .( YTOP-S.O) ) GO TO 3|0 
vrnnaYTnp* to.o 
GC TO 320 
310 CONTINUE 

YBA YTOP— 60 » 0 

IF (VO PC (I .1) .LT .VPAS r ) VORDC I . I )=YBASF 
J3T CONTINUE 

WPITE (3.4042) I T ITL3( I). 1=1, 20) 

WFITe (3.1076) ( 1 , XCPDl I , I ) ,V0PD(1 .1 ) .1 = 1 ,N> 

IF ( IPLOTA.F.O.O) Or TC 30 
CALL PICSIZ (14.0,7.0) 

CALL OR AFT ( XLNGTH. XL EFT. XR I GMT , X INC .NDECX, XL ABC I. . XCIFO, XllUMHY. YLNG 
7TM, YOASP .VTOP.YI NC , NOFCY , VLABFL « YOPf) ,Y DUMMY . N5 I NGV , M A XCV S.NCURVS . M 
TAXnrS.N^TSl .NOTR2.NnTR3.NPTS4.NPTS5.KINU.LINES.TlTL* ) 


■OPS Tr; 


r»r»r» rv r» nnonnn n ft « rtonn«^'V nnn 


so continue 


or mo i=i . n 

"CHHC1 . I )=COK JG< Gl?( 1 ) )*GI2(1 >/<GI II ll*«?(lll 
PCMHI2 . 1 )=-CON JGIGI 3< 1) > «G I 3 t 1 1/ ( G 1 1 I I)*G33( II) 
PfO'Mj,! ) = CON JG ( G? 1 ( I ) )*G23( I )/(G22( I ) *G33< I ) ) 
300 CONTINUE 




Of n T C*°05 rcn TKtr ORDINARY CCHFPENCE BETVEFN CM| t CM2 

rnLi ow this capo: 


VBASr=0. 0 
ytdp=i .o 
V 1 NC = 0 . 2 


no io«i i*l . n 

V0PD( I . I I =PFAL( PCOMtl . 1) » 

Knon ( 1 . I )=FPFO< I ) 

O00 CrNTINUF 


WFITE <3, 4350) ( T 1 TL At I 1,1=1.20) 

fcrtTC (3, IOTA) ( I.XOPDII .1 ) » VOPO <I«I)»I=I«N) 


IF ( IPLTTO.CO.SI GO TO 43 
CALL PICSI* 114.0.7.0) 

CALL OR4FF ( XLNGTM.XL^FT.XRIGHT.XINC .NDFC X , XL AOF-L , XO BO , XOUMMV. VLNG 
7TH.VRASF .VTOP.V INC . NOECY , YL *EC3 . TOP 0 » V DUMMY . NS I NGV »M A XCV S . NCOP VS, M 
TAXPTS.NPTSt , NPTS2. NPT S3 « N°T S4 .NPTS5 » KI NO , LINES, TITL4 ) 

40 CONTINUE 


PLOT CAPOS FCP thf ORDINARY COHFRFNCE BETWEEN CHI C CH3 
FOLLOW THIS capo: 

DO I 9& ! I=I.N 
VORD ( I , t I^RE ALI °COH 12,1 ) > 

XOPOI I . I l=FREO( 1 ) 

I PA I CONTINUE 


NRITF (3,4051) (TITLSd ) .1*1.20) 

WRITE n.l«7F) ( I. XOPD( I • I ) .VOPfH I .1 >.!=!. N? 

IF I IPLOTP.ffi.O) GO TO 50 
CALL PICS1Z (14.9.7.3) 

CALL GR AFF ( XLNGTH , XLFFT . X F I GMT , X INC .NDECX , XLABEL, XORD. XOUMMV, VLNG 
7TH, VnASF, VTOP.V INC . NOrC V . VL AO e 3 . VOPD • VDUFMV . NS I NGV . M AXCV S • NCURV5 • M 
t»XOTS.N»TS| .NOTS2.NPTS3.NPTS4.NPTS5.KIN0.LINCS.T ITL5) 

50 CONTINUE 




oonfro rt n n'W 
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P| HT FAROS FOR TM 1 " OPOINARV COFFRFNCF B r Tt» F FN CM? t CM3 

rciuunw this caro: 

f>o ion? i = i , m 

VCRDI l . I ) -or ft. (PCQM( 3, I I ) 

XPODf 1 . I l=FPFO( I > 

»C«2 CONTINUE 

priTE (3.405?) (T)TL6( 1 > .1=1 .20) 

PFITF (3.1076) C I. X0P(>( I . I ). VORDU . I ) . 1 = I ,N| 

)f ( I "M.otp . ro .0 ) r.o to no 

CALC PICSI7 (14.3.7.0) 

FALL C.CArr ( XU NGTM i XL cr T • X P I GMT . X I NC .NDF CX . XL AOEL . XOFO . X DUMMY . VLNG 
7TH, VP ASF .YTOP, V I NF . NDFCV . VL APF. 3. VPPO . YDUMMY . NS 1 NG V , M A XCV S . NCUP V5 . M 
TAXPT S.NPTGI ,NI»TS2 , NPT S3 . N^T S4 .NPTS S , K I ND, L I NFS . T t TL 6 ) 

60 CONTINUE 


TMF CAROS FOR COMPUTING THE PARTIAL COHCPENCES FOUL OP THIS C4F0J 


C 

2000 

c 

c 

2100 


c 

c 

2200 


!P = -» 

00 7000 J--I.N 

1 F ( IP) 2 I 00, 2200. 2300 

GAF.( I l=GYY( J) 

C.AGI 2)=CCNJG( GY I ( J ) ) 

G A G ( 3 )=CCNJG(GY2( J ) > 
GAG ( 4 | =C0NJG(GY3( J) ) 
GAG(S)-GYMJ) 
£AG«6)=G1) (J) 

GAG( 7)=CCNJG(Gl 21 J) ) 
GAG( P > = FON J G ( G I 3 ( J I ) 
GAG( 9) =GY? ( J) 

GAG( 10 ) = GI2( J > 

G A G ( l I ) = G??( J) 

GAG( I 2 )=C r )N JG ( G23( J) ) 
GAG( 13 ) = GY 3 ( J I 
GAG ( I 4) =GI T( J ) 

GAG ( 15)=G?3( J ) 

GAG( 16 > = G33( J ) 

GO TO 2403 

GAG( 1 )=GYY( J) 

GAG( 2)=CCN.)G(GY?( J ) ) 
GAG ( 3 ) = CPU JG (GYI(J)) 
GAG( 4 > = CCNJF,( GV3( J ) ) 
GAG(5)=GV?( J) 
G4G(6)*G2?( J) 

G*G( 7 I = Gl 2 ( J ) 

G4G ( B) =CCN JG( G?3( J ) ) 
G4G( 9 )rGY I ( J ) 


w: 
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c 

c 

*300 


c 

c 

2A03 


c 

c 

c 

c 

c 


gag* i,i jii 

GAG* I I l«CI I ( J) 

GAG* 13 » kCONJG* C. I 3* J| » 
GAG* m«GY3* J> 
gag* ia)«g;>3* j» 

CAG* IM«G|3( J> 

GAG * I A I *G33 I J I 

GO to ?400 

GAG * I >nGYV< J I 
GAG I 2) «rciNjr,< GY I* J M 
r,AG< 3)nC('Njr,( r.Tp( j J > 
GAO.* 4 )®cr KJGIGV |< ,| > | 
GAG* *51 «3Y3* J> 

GAG* 

GAG *7|»G?3*3 I 
GAG* nlsC.tl* J» 

G»G* Q t J I 

GAG* I 0 I *0f IN J C * G P 1 1 J M 

GAC* I I >*0??* J> 

GAG* |?)tGI?< J l 
GAG* t 3>*GV* * J ) 

GAG* I 4 > = CnNJG«GM( J J » 
Gag* is ) jjCONJG* Gm* J1 I 
GAG* If. »«r, U* J> 

GO TO 

SW* l|nC AG* I ) 

«VY< ?l*fi AG* ? 1 
tVV* 3) *GAG* 5 > 

SVY* * |«f AG<6 | 

GPY* I IwGAG* 3) 

SPY* 2 )«GAG< A > 

5PV* H »GAG* T I 
SPY* AlrGAG* ni 
SYO* l ) * GAG* a » 

svotamr ag* to) 

SYP* 3 ) «G AG* 1 ‘1 1 
S YP * A >■ GAG* I A I 
SPP* 1 ) »GAG* I I I 
Sr>r>( ? ) m G A C> ( I ?1 
SPP*3>«GAG(I5J 
•P"t A | wG AG* IA1 

100 = 10*2 

CAU, C*1|NV*SP°. 2,Ot.D2 I 


10*01 .NP. 0 . 0 . AND. 02. Nr. 0.01 GO TO 2A70 
WPITFI 3 .54 AO » J. 1 PC .01. 0?. 

PCOM* J«C. J>*CMOLH« O.O.rt.OI 
CC TO 2C50 


-TRTpr 


nonruTo o o n(inrt«n 
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call rMPrrn( 5Y-*.‘?r*r’.r'VNP,?> 

CALL CMPFOC FVNP, SOY .5YNPV.21 

c 

nn ?soo i = i .4 

GXYO< 1 1 = SYV I 1 1-SYNPV I l » 

2500 CONTINUE 
C 

2553 CONT INUF 

GXVD*GXVP( 1 |«fiXYP( 4 1 

<ixYR*CAnsir.xvr»> 

c 

(S»"FS0=PFF*P cr F 

cFroP=prFso*PCFFo 

c 

IF ( GXYF.LT. RFFOP » r.o TO 2000 

c 

FCnMII°C , J> = < REALS GXVP< 3 1 1442 ♦ AIM A6< GXYPC 3» »4*2|/GXYO 
GO TO 2SS0 
C 

2 P 00 FCO'M IPC.J»=CMDLX<0.0.0.01 

c 

2° 53 CONTINUE 
TOM CONTINUE 

c 

»o= 1PM 

irnp.c.F. 2 ) c,n to moo 
c 

CO TO 2003 

c 

3100 CONTlNUr 


Cl OT CAPOS TCP TFE PAFTtAL COHERENCE OFTMEFN CHI t THE OUTPUT 
CHANNEL FOLLOW T HI S CARDS 

00 1 =53 1=1 ,N 

VOPDC l . I 1=PRAH CCOH1 I . I ) 1 
XOPOJ i . I >=FCE0< I I 
1C AT CONTINUE 

wrlTT C T .4060) (TITL1 < T 1 . 1=1 .20 1 

YriTc 13. 1 I , XCC0< 1 . t > .VOPDU , I > ,1 = 1 .N! 

IF I 1"LOTC.CQ.O ) GO TC 70 
CALL P I C F I 7 <14.0.7.01 

CALI craft < XLNGTH, XL TFT. XR IGHT. X INC ,NDEC X . XLABEL . XORO . XDUMM Y . YLNG 
7TH. Y0A5 r .YTOP.Y INC .NDFCV.YL ARr2.YCR0,Y DUMMY, NP 1 NGV , M AXCV 5, NCURVS. M 
7AXPTS.NPTFJ , NPTc; 2 , NOT 53. N^T S4 • NOT 55. K 1 NO .LINES. T I TLE 1 
70 CONTINUE 


RLOT CAROS FOR Tt-F PAFTIAL COHCPFNCF BFTXFFN CH2 f, THE OUTPUT 
CHANNEL FOLLOW THIS CAROS 


nrto r»nn r> r* 


156 


c 

no i op 4 ] = i.n 

vnr?o(i . n = rr •KPCdMj. i ) » 

XPOOC 1 . I ) = FPE<11 I I 
I ®8 A CONTINUE 

c 

WPITr (3.A0M) < T I Tl .'’(11.1 = 1.70) 

WFITT M,l«'76» < I. XOTOM . H.VOPOI I * I I . 1= I .N) 
r. 

IF ( ir>LOTC .FO.O ) GO TP *0 
C ALL P1CFJ7 (Ifl.D.-.O) 

CALL GP.AFT { XLNGTH.XL^FT.XPIGHT , xINC . NDF.CX « XL ADEL . *0 RO . X DUMMY . VLNG 
TTM. V»*ASC .YTftP. Y INC .NOECY , YL ARF2.VC.Rn. VOUMMY.NSINGV . M A XC V S .NCUPV5 • M 
7AXP7 S.NPTM ,NOTS2.NOTS3.N°T<54.NPTS5.K!N0,L!NrS. T ITl.? » 

AO CONT IMU~ 


CLOT CAnor rcn TFF partial cchffencf between CH3 t THE OUTPUT 
CH\NN r t nLLOW THIS f AR O I 

on i “="•? .n 

vnoni i . I )=PF ALI fcohi 3. i » ) 

xnmi i , I i^rrEnc 1 1 

I 9*5 CONT I >4U'- 

WFITT (3.406?) (TI TL3C ! » . 1 = 1 .20 I 

XP I T= ( 3 . 1976 ) ( I , xonol l . n . VCRP( ! • I ) .1 = 1 .N) 

l= c inLPTc.ro. o I on to «o 

CALI ® I C S I 7 (I4.0.T.C) 

CALL C.PAFF (XL NOTH. XL EFT. XP I C.FT , X I NC . NDFCX . XL ABFL . XORD . XDUMMV. VLNG 
7TH, YnASP .VTnP.YINC .NOFCY ,VLAB=? .YCRO.YOUFMY. NS INGV.M AXCVS.NCUPVS.M 
TAXPTS. NDTSI .NPTS2.NDT S 3 . NPT S A . NP T S5 • X I ND .L INES. T ITL3 I 
OO CONMNUF 


ir < IPLOT A.FC.O. AND. iPLOTn.eo.O.AND.IPLOTC.EO.OI GO TO 9009 
CALL PICS12I 0.0.0. 0» 


9999 STOP 
END 


SUHF'OUTlNr CMINW(SPP , Is ,01 ,02 » 

P1NFNOC.N SPP ( * I « A ( ♦ > ,B< 4) ,<:<0 ,F (4) ,G<4) 

fOMOLFX SPP 

no i o i=i, 4 
A ( I )="r*L(SPP( I > I 
m n = OMAG(spp( 1 1) 

P( I l = 4( I ) 

13 CONTI N'JP 

TALL MtNV( r ,h.n» 

r 1 =n 

CALL MPrO r 'tn ,p ,F ,N) 
call Mr>sr.n(p ,n ,o,n» 
on ?o l-i,4 
a» I ) = A( I i »r.( I > 

P.3 CP NT IMU<= 

CALL ‘AlNVIA.N.n) 

r>?=n 

CALI MOPOTt A.F.r.N) 

on 30 1=1,4 

ypr> ( j I-=r Mr>L K( A < 11,-00)1 
39 CONTINUE 
pcfupn 
Fin 


sunonur IN'- MINVCA.N.C) 
OIM^N^jrN Mt) iLU ) >M(4 I 
D« t .0 
W.*-N 

nn no Krj.N 

N« = NK AN 
t»K)=K 
MCK »=K 
KK=NKAK 

ntGA=A(KK) 
no 23 J “K . N 

I C = N*C J-I » 

OP ?0 I=K.N 
f J= I 7 A I 

10 in MiMnioAi- Ann caiijmi 

• 5 p ir,A*M t .1 > 

C «•*•* = 1 

M(K> =J 
23 CONTINUE 
J=LCK) 

lrt.l-K) 35.35.25 
25 risK-N 

OP 30 I-l.N 
K I = K I AN 
NOt.D=-A<Kl ) 

JI = K l-K A J 
A t K I ) = AC Jl ) 

30 AC J I l =Hr,t_r> 

35 t=Mt K I 

ir c i - k > 45.A5.in 

36 J® — N * ( 1-1 I 
no 4 0 J = I .u 
jK=Nka j 
JI=JPAj 
HOLO=-A( J*C> 

AC .IK l = AC .It I 

♦o A(jt) =nnu-> 

45 irmir.AI 40.44,45 

46 C*3.0 
PFT'JBN 

4a nn 55 i = i .n 

TFCl-K) 51,53,59 
50 1K = NK.AI 

AC !K )=AC I* I/C-R1GA > 

55 CONTINUE 
on f 5 1 = 1 , N 

!K=Ntr* j 

H0LD=AC IK) 

I J- I-N 

bn 65 j=i ,n 

1 J=l JAN 

JFCI-K) 60.65.60 
60 CFCJ-K) 62,65.62 
62 KJ*»J-!4K 

AC f J|*HCt 0*ACKJ)4 M I J) 

65 TOUTIN' CI- 


IS, 20. 20 


K J = K-N 

on - s j = i , n 

* J=K J*N 

IFIJ-Kl 70,75,70 
70 A1KJ1 = MKJ1/PIGA 
75 rnNTlNUr 
D=-n*m ga 
MKKI rl ,7/FUOA 
00 CONTtNU' 1 
K*N 

I JO K = < K-l > 

I r ( k i ifo.ino.lC 5 

I JO I=l(K) 

irii-Ki i po, i ,?o ,i on 

JOn JJfNMK-ll 
JP = N*( i -i ) 

on no j-=i ,n 

JK= JO* J 
WLl> = M JK ) 

Jl=JPK J 

AC JK )=-A( JI ) 

110 MJI> =HOL r > 

120 J=MCK> 

ire J- K l 100,100.125 
1 ?5 K1 =K-N 

DP 130 1=1. N 

KI=K1*N 
HCiLP=»<K 1 1 

JT=K T-K4 J 
MKI 1=-A< J1 1 
130 A( Jl 1 “HOLD 
GO TO 100 
150 P^TIIPN 
TND 
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pumpout ! N r :npro(».n,r.N| 
CJHrNStCN A( « ) , o ( , r (* | 
COWPLF* A.B.R 
|F = o 
IK* -N 

on 10 K=iiN 
i*=ik*n 
or to j=i,n 
IF* irM 
jf* J-N 
10=1 K 
PI IP 1 = 0 
on to i*i. n 

JI=J1»-N 

IP= i n* J 

10 PIIP|=BI |T > ♦ SI J1 l*H I I O' I 
PFTUPN 
fNO 


•UPBnnTi^ Mpont ia.o.b ,M 
MMrNSIfIN A( M «Rl 4 l .PI *1 
lF = 3 
IK®- N 

On I 0 K a I • N 

IK- IK+N 

on io j=i,n 
IP* jptl 
J! = J-M 
in* IK 
S| IP >=D 

or io i=i .n 

J 1 =JI»N 

|R~ 1 R+I 

10 PI IP >*PIIP»*M 
fftusn 

FNO 
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11.4 CIRCXCOR 


The two-channel program "CIRCXCOR": 

1) removes the overall mean from each of the two time series sample 
records and amplitude weights each time series with a Hanning 
2 

(sin ) function. 


2 ) 


computes the spectral values for the two-sided spectra G , 

G via the FFT. XX 

xy 



3) averages the above quantities by repeating steps 1 and 2. 


4) "edits" the spectral values in chosen frequency bands (digital 
filtering) . 


A 

5) computes the frequency-filtered values for R , 

the inverse FFT and the convolution theorem [3, 
cross -correlation . 


R , and R via 
yy xy 

4]; i.e., circular 


6) computes the correlation coefficient, p . 

7) outputs all above-mentioned averaged quantities on a printout. 

A 

8) automatically will scale the plot axes and then plot p 

xy 


uvouuuuwuuuuuuuuu 
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integer *2 infil.junk 
REAL44 msdat.mscup 
COMPLEX FFTI .FFT2.GXY.RXVCX 
DI MENS ION JUNK (201 

DIMENSION M5DAT (201 > MSCUR! 20 ) . DEL 1 99 1 , SHE IL ! 23 4 8 ) 

DIMENSION XDRD! 1,12011 .TOR 01 1,12011 , XDUWMY ( 1 .12011 .VDUMMVI 1 .12011 

DIMENSION XLA8EI.120I .VLABEL (211, TITLE! 29 I , K I N3 I SI 

DIMENSION RXXI2048) .RVY( 2048 ) .PXY 12048 > , GX V < 23 4 8 > , XF RE 01 204 SI 

DIMENSION T AU( 2 048 I.INFIL<4 096>,CHH294B),CH2!2048),INIT(100I 

DIMENSION AI 20481 ,B( 2048 ) .EFT I <2048 > , EFT 2(20481 

DIMENSION NLOWI I0I.NHIH! 1 0 1 . GXXC 204 8 I . G Y Y( 2048 1 .RXYCXI2048) 

c 

50 FORMAT !* I *.////. 40X .20A4) 

100 FORMAT ( 40 42 I 

155 EOPMATI'OPOINT NUMBER • « 16 » • NAS A VALUE OF • , 1 1 2 ) 

500 FORMAT I 20 A 4 I 

1058 FORMAT I • 1 • ,///, • CHANNEL I 3. • ZERO-MEAN INPUT :•,//.< I OF I 2. 1 1 1 
3000 FORMAT!* ITERATION NUMBER', 13.* EXPONFNT OF*.F6.2,* DEL= • . F9 . 1 .// I 
3311 FORMAT! • 1 «.///. • CHANNEL * , I 3 . » DATA AFTER DECAY CORRECTION:*.//, 
III9FI2.il) 

4300 FORMAT!* 1 IBAND. FREOUCNCY, GXX, GVY. GXY ! COMPLEX I *,// , 

1 12 ( I 1 1 ,F9.t . 4E13.I ) ) ) 

4345 FORMAT!* *.//,* ON FILTER PASS*, 13.* THE AVERAGE X-SPECTRAL VALUF., 

1 GXYAVE, IS* ,FB»2»* UNITS. 0F*,I5.* FILTERED* ./»• SPECTRAL LINES.* 
2.15.* WERE ABOVE THE DXYTST VALUE 0F*.*tb.2,* AND MERE ADJUSTED.*./. 
3* THE POST-WHI TENED X-SPECTRUM IS AS FOLLOWS! • .//. • IBAND FREOUF.NC 
4 Y . GXY I COMPLEX )*.//) 

4346 FORMAT I2II25.F9.I .2F10.I )) 

5100 FORMAT!* MSOAT ! * . I 2 . » ) **,FI4.i) 

5200 FORMAT!* THE MSV OF ALL DATA =*.F9. !.///> 

5500 FORMAT < * NSCUR! • .12. • I =».FI4.5I 

5600 FORMAT! • THE MSV OF THE EXP CURVE =*.F0.5.* THE REQUIRED GAIN IS 

» thus* .fs.ii 

6000 FORMAT! • 0* .//. • FILTER PASS * . I 3.F8. I * • HZ TO*.FB.l«* HZ*.//, 

216,* INITIAL POINTS ARE OMITTEO FROM EACH CHANNEL.*,//) 

6010 FORMAT 1*01. TAUIII-ICH1 LEAD IN MSEC). RXYlt). RXXll). AND RYYlI). 

2 1*1 TO 2*N* , // .1211 II. Flo. 2, 3F 13.4))) 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

XDELAY SPECIFIES THE POS. t NEG. DELAY I IN MSEC! FDR THE PLOT. 

IF JFLAG = 0 . NON-ESSENTIAL PRINT OUT IS DELETED IE. GXX.GYV... 

IOMPSK SPECIFIES DJMP TO BE SKIPPED DUE TO FILE ERROR. 

NOMIT INITIAL POINTS ARE OMITTED FROM EACH CHANNEL. 

IF IN9DW * 0, FILTERS WILL BE ABRUPT TRUNCATION. 

NLOWI JPAS) AND NHIHIJPAS) SET BAND-PASS FILTER CUT-OFF 
FREQUENCIES! NLOWI JPASI-l AND NMIM|JP4SI*N FOR NO FILTERING. 

SPECTRA ARE POST-WHITENED UNLESS !WHITE»0. 

DECAY CORRECTION: SPECIFY IDECAYsl 


nnrtnnnnnoonftrtonnn 


uuu 


C C 

cccccecccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 

XDFL AY=I .0 
JFLAG=0 
N=512 
I A VG 5=64 
1 DMPSK = 0 
N0NIT=3 
IN0OW=0 
3PLRT=20000.0 
Pt=3. 1415927 
C 

NPAS=2 
NLOPI 1 > ■= I 
NH I H ( l )=N 
NL0tM2> = l04 
NHI HI 2 1=206 
C 

I WH I TE=0 
IDECAY=0 
NPT$BL=125 
NBLOKS=B 
DELXPO= I .0 

N2H*N4N 
KAOP=I024/N 

data label capos go hfre: 

DATA XLABEl/*TAU *,*- CM* . • ANNE • . »L t • • ‘LEAD* • • IMS* 

OATA YL ABEL/ *COPR*»* FLAT • » • I ON • , • COEF • , • F IC l • . • ENT • 

C ” 

FE AD ( 1.500) (TITLril ) .1*1 .20) 

C 

OTAU=lOOO ./SPLRT 
JAV*0 
IOUMP=0 

c 

DO 70 1*1, N2H 
GXX ( 1 ) = 0.0 
GW I 1 ) =0 .0 

GXV ( t ) = C MPL X(O.O.O.O) 

70 CONTINUE 
C 

1NI«0 

DO SO 1*1 • K A DP 
In I T 1 1 )* Ins 
INI=INI4N«-N 
MO CONTINUE 

c 

90 REAO (4.100) (INFILI I), 1=1,4096) 

C 

NOTOK*0 

00 14B 1*1.4096 

1F( INFILI D.LT. 1024. AND. INFIL(I).GE.O) GO TO I 46 


X* 


uuu u u u uuu 


NPTOK* 1 

WFtTF!3.l55)I.INFlL! I) 

148 CONTINUE 

IF I NOT OK .EO. I ) GO TO 99 99 

220 IDUMP= I DUMP* I 

IFf IDUNP.EQ. IDMPSK ) GO TO 90 

K=0 

IFND*204 8-NOM I T 

00 240 1*1.1 END 
CHI I 1 >=INFIL! I ♦N0M IT ) 

CHS I n* INF IL I I 4 20484 NOMI T> 

240 CONTINUE 

1F( IDECAV.EO.O) GO TO 1225 

SET 2ERO— MEAN AND DECAY CORRECT! 

RMV1=0.0 
PMV2-0.0 
no 1221 1=1.1 END 

RMV1 = BMV1*-CM1 ( 1 I 
PMV2»BMV2«-CH2( I ) 

1221 CONTINUE 
PMV1*RMV1/IEN0 
PMV2=PHV2/ IEND 
DO 1222 1 = 1 . I END 
CHI { I ) *CHt t I I-RMVI 
CH2II l=CH2! I I-RMV2 

1222 CONTINUE 
C 

IOCM* 1 
NPO I T= I END 
C 

DO 1 056 I=i . NPO I T 
SNFILM > = CMIC I ) 

1056 CONTINUE 
C 

1057 SUMM 1*0.0 

c 

WRITE! 3. 10 58) IDCH, ( SNFIL! I».I*1 .NPOIT) 
C 

DO 1110 J*I .N0LOKS 

c 

SUM=0 .0 

DO 1010 1*1 .NPTSBL 
JRUN»NPT SOL* I J- 1 1 ♦ I 
SUM-SUMMSNFILI JRUN) )**2 
1010 CONTINUE 

c 

NSDAT ( J> =5UM/NPT*BL 
C 

WRITE! 3.5100) J.MSDATiJ) 
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c 

SUNMlsSUMMI+MSDATl J) /NllLOKS 

c 

1110 CONTINUE 

c 

WPITEI 3. 5200 ISUMMI 

C 

*PO=0.0 

ITNO* 1 

1150 SUNM2=0.0 

c 

00 1310 J=l t NOLOKS 

MSCURC J) *1 .0/5X°< 2. 0*XPO*(NPT*BL*1 J-l ) ♦NPT*OL/2 1/SPLRT 1 
C 

WRITE (3, 5500 ) J.MSCURC J ) 

c 

5UHM2=SUMM2*NSCUR1 J > /NBLbKS 
1310 CONTINUE 

GA IN=SUMN|/5UMM2 

c 

WPI TE<3 .5600 15UMM2 .GAIN 
C 

OFLC ITNO »=0 .0 
C 

OO 2000 I=l.N0L0KS 

DEL < ITNO l = DEL! IT NO I +ABS1MSDAT C Il-GA IN*MSCUR< II ) 

2000 CONTINUE 
C 

WRITE13.3000 > I TNO, XPO.DEH ITNO) 

C 

IF! ITNO.EQ.l ) GO TO 3100 
IF < 1 TN0.GE.99) GO TO 9999 
C 

CK0EL= 1.001 *DEL I ITNO-1 ) 

IF1DEH ITNO) .GT.CKDEL) GO TO 3111 

C 

3100 !TN0=ITN04.| 

XPO=XPO*OELXPO 
GO TO )150 
C 

3111 XPO=XPO-DELXPO 
C 

00 3211 I=l,NPOlT 

SNFILl I >*SNFHL1 1 )* EXP 1 XPOA I /S°LP T ) 

3211 CONTINUE 
C 

WRITF1 J.33I I )IOCH. (SNFILl I ). t*l.NPOIT) 

C 

1F< IDCH.E0.2 ) GO TO 6103 

C 

DO 6102 I=l,NPOIT 
SNFILl I )*CM2 1 I ) 

6102 CONTINUE 

C 

I DCH*2 
GO TO 1057 



m m 

A ' 


I? 
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6103 CONTINUE 

EN3 DECAY CORRECTION. 

225 JAV=JAV*I 
KaKM 

nix=i 

EVAl =0.0 

DO 1 230 1=1 .N 
A ( I >=CHI I I ) 

81 t 1=0.0 

EVAL=FVAL*A( J ) 

230 CONTINUE 

EVAL=EVAL/N 

CHANNEL I MEAN VALUF CORRECTION 

DO 1231 1 = 1 .N 
Am = (A< 1I-EVAL) 

231 CONTINUE 

DO 1232 1=1 • N 
I2H=C I+N 
A< I2H)*0 .0 
B( I2H)=0.0 

232 CONTINUE 

CO TO 1275 

1250 FVAL*0 • 0 

DO 1255 1=1. N 
Alt) =0 .0 
010 = 0.0 

1255 CONTINUE 

DO 1257 1=1. N 

12H»I+N 

AC I2H)=CH2(I) 

Bl 1 2 H) =0 .0 
EVAL*EVAL*ACI2M) 

1257 CONTINUE 

EVAL=EVAL/N 

CHANNEL 2 MP AM VALUE COPFSCTIOM 

00 1256 1=1. N 
I2H*ITN 

A«I2Ji>=AlI2H)-EVAL 

1256 CONTINUE 

1275 CALL FFTI A.B.N2H.N2H.N2H .1 ) 


nor* n non nnn 


c 

tFCNlX.E0.2l CO TO 1400 

c 

DO 2001 1 * l • N2H 

FFT1 (I >»CMP|_X( 41 II. 0( III 

GXX( I >«GXX< t >4REALIC0NJG<FFT11I » l*FFTI (111 

2001 C0NTINUF 

c 

HI X = 2 
CO TO 1250 
C 

1400 00 2002 1-I.N2H 

FFT21 I >=CMPLX( AC ll.eilll 

GVT 1 I)=GVYl n+REALtCONJGlFFT2H) I4FFT2CI )> 

C.XYf I»=CONJG(FFT 1(1) >4FFT2f I l»GXY( 1 ) 

2002 CONTINUE 

c 

IFf JAV.EQ.IAVGS) CO TO 4000 
IF(K.GF.KADP) GO TO 00 
GO TO 1225 
C 

4000 CONTINUE 

END AVERAGING 

ONQPN* 4 *N*N* I AV GS 
00 4100 1*1. N2H 
GXX( I ) * GXX ( I 1/DNORH 
GVYC I l*GYY( I l/ONORM 
GXV( I >*GXV( 1 »/DNORM 
XFREOl I )=< !-t )*SPLRT/1N*N) 

100 CONTINUE 

AUTO AND CROSS SPECTRA ARE COMPUTED. PUNCH 0* STORE OUTPUT. 
IF ( JFLAG.EO.O ) CO TO 4304 

VP I TE( 3* 4300 ) 1 I . XFREQ C II. GXX C I ) . 6VT ( I > . G XT ( I ) . t * I . N2H) 
BANDPASS FILTER THE TMO-SIOED SPECTRA... 

4304 CONTINUE 
,IPAS=0 

4305 CONTINUE 
JPA5* JPAS4-1 

C 

DO 4310 1*1. N2H 
RXXII )*0.0 
RTY( I )*0.0 

RXYCX1 I >*CMPLX(0 .0.0.0) 

4310 CONTINUE 
C 

ItOV-Nt.OH(JPAS) 

IHIH-Nl ilHI JPAS) 

C 

OO 4320 I-ILOW.IHIH 
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c 

N2I=N2H*2-I 
RXX< I (=GXXC I ) 

RXXIN2I I =GX X ( N? I I 

Rvvc I i=gvv< I ) 

RVVIN2I (=GYY IN21 > 

RXYC XI n *GXY< | ) 

PXVCXIN2I (=GXY(N2! I 

c 

IFINLOWC JPAS( .NE.II GO TO *320 

C 

NA1=N+I 

PXX(NA1 (=GXX(N* I I 
RVY (NAI > = GVY CNA1 I 
PX VC XI NAI I “GXY ( NAI I 
*320 CONTINUE 

c 

C POST-WHITE* THE BAND-LIMITED DATA... 

C 

IF ( (WHITE. EQ. 01 GO TO *3ST 
C 

NGAINS=0 

GXVAVE=0.0 

c 

DO *335 1*1. N2H 
GXVAVE*GXYAVE*CABS(PXVCX< I I » 

*335 CONTINUE 

c 

1DNM=?A( ! HIM- 1 LOW* I ) 

GXVAVE-GXYAVE/ IONM 
OX VT ST*GXYAVE*I 0.0 

c 

DO *3*0 1=1. N 
N2I=N2H*2-I 

CHECK=0.5*<CABS(RXVCX| I I (*C ABSIR XYCXIN2I ( ( ( 

C 

IFICHECX.LT. OXYTSTI GO TO *3*0 

c 

ngains=ngains*i 

GA I N=OX VTST /CHE CX 
C 

RXX ( 1 l-RXXII (AGAIN 
RXXIN2I (*PXX(N2I (AGAIN 
RVYf I ( *RYY 1 1 >*GAIN 
RVVC N2t (*PVYIN2I I AGAIN 
RXYCX( t (sRXYCXI IIAGA IN 
RXYCXtN? i >=RXVCXJN2? 

*3*0 CONTINUE 

C 

IONM* I DNM/2 

C 

WRITE! 3.43*5 I JPAS.GXYAVE. IDNM.NGAINS. DXVTST 
WPI TE I 3*43*6 (If. XFBF OC I I .PXYCXlI ) .1*1 .N2HI 
C 

C END POST-WMITEN. 

C 






non n nr»n 


4350 IF! INDOW.CO.OI GO TO 4360 

TAPER THE FILTER WINDOWS. 

IL0W=NL0Wt JP AS > 

I H t H=NH I HI JPAS ) 
OELF=!HIH-!LOW 

00 4355 I=ILOW,lHIH 
TAPER-S1NIPI *( I - ILOW ) /DELF ) 
N2I=N2H*2-I 
PXXI I )=RXXt I J4TAPEP 
PXXIN2I )=PXXIN?1 >*TAP C R 
RVYtl ) = RYV ( I I *T 4 PER 
»YVt N2I I=RVYIN2l ) 4 TAPER 
RX VCX I I l=PXVCX< I I FT A°ER 
RXYCXI N2 I I =RXYCX(N2 I ) • T A PE R 
4355 CONTINUE 

END TAPFR. 

4360 DO 4400 1=1 ,N2H 

A I I >=KEALtRXYCXI I > I 
Bt I )=-A I MAGI RXYCXI I ) I 
4400 CONTINUE 
C 

CALL FFT I A » B .N2H.N2H . N2H » ♦ I I 
C 

DO 4410 1*1 .N2H 
RXY(I)*Att) 

4410 CONTINUE 
C 

DO 4600 1=1 .N2H 
AI 1 1 =P YYt I ) 

Bt I ) =0 . 0 
4600 CONTINUE 
C 

CALL FFTI A.B.N2H.N2H.N2H.-1 ) 
C 

DO 4710 I=t.N2H 
RYYII I*A I I > 

471 0 CONTINUE 
C 

DO 4800 1=1 • N2H 
At I I =RXX II) 

BC 1 1=0.0 
4800 CONTINUE 
C 

CALL FFTI A.B.N2H.N2H.N2H.-I » 
C 

DO 4910 I * I » N2H 
RXX ( t I* At I I 
4910 CONTINUE 
C 

DNl-RXXtl) 

DN2-RYVI |> 
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ON3=SDRT(DNI*DN?> 

c 

on 5000 1 = 1 ,N 

FXX( I >=RXX< I 1/DN14FLOAT ( N> /(Ntl-II 
RY V 1 1 l=RYYt I )/r>N2*FLOAT( N|/(N* 1- 1 ) 

PXY ( I >-FXY< I ) /DN34FLOAT ( N)/| 

TAUI 11 = 1 l-N- I l*OTAU 
5000 CONTINUE 

c 

NIzMfl 

on 5010 I=N1 ,N2M 

»XX(!I = RXX(I )/DNI*FLOAT(NI/(I-N) 

»VV t I )=PYY< I |/ON2*FLnAT(N)/< I -Ml 
«XY< I)»PXY( 1 )/DN3*FLOAT(N)/(N»N*l-I ) 

TAUI I 1 = 1 l-N- 1 >*DTAU 
5010 CONTINUE 
C 

•PITF 13.501 (T1TLEIL).L*I.20I 

WR1 TE< 3.6000 I { JPAS. XFREOlNLOWl JPAS) I .XEREOINMI mI JPAS* I .NOME T» 
WPITEI 3.601 Oil I tTAUI I ) *RXV | I ) . RX XII) ,PYY< 11.1=1 . N2H) 

C 

IF I JPAS.LT.NPASI GO TO *305 

PLOTTER ROUTINE POP CIRCXCOR. 

MAXPTS*! 20 1 
X I NC = X DELAY / 4 .0 
XLEFT--XDELAY 
XP I GHT =X DELAY 
VTOP=0,0 
YBASF=0.0 

•N» x THE » OF POINTS PLOTTED 

M2-XDELAYM 1 I .0/12.0 1*1 SPLRT/1000.0 1 
INFG=N— M2 
IPOS=N4M2 
M=2*M2+ l 
J»0 

c 

DO 6050 1=1 MEG . I POS 

TORDI 1 . J»=RXY( I » 

XORDII. J»=TAUI II 

c 

IF I YORDI t . J) .GE.I .0) YORDI I . Jl«1.0 
IF I YORD (t .JI.LE •— 1.0) YORDI 1 . J 1 . 0 

c 

6020 CONTINUE 

IF IYORDCJ .Jl.LE.YTOP) GO TO 6030 
VT0P»YT0P*-0.25 
GO TO 6020 
C 

6030 CONTINUE 

IF | YORDI l.JI.GE.YBASEI GO TO 6040 
YBASE-VBASE-0.25 
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GO TD 6030 

c 

60*0 CONTINUE 
6050 CONTINUE 
C 

IF I-YBASE.GT .YTQP ) YTOP*-VBASE 
IF ( YTOR.GT.-YOASE) YBASE^-VTOP 
C 

VINC=YT0P/'5.0 
KINO I n=-l 
c 

CAUL PICSIM15. 0.9.0 > 

c. 

CALL GRAFF <9 .0. XLEFT . XR IGHT . X INC. 3, XL ABEL. X03D. XOUMMV . 

15.0. VBASC,VTOP.YINC.2.YLABEL.YORD.VDUMMY * 

21.1 , l.MAXPTS.M.O.O.O.O.K I NO .0. TITLE ) 

c 

CALL P1CSI MO. 0*0.0) 

c 

9999 STOP 
END 




173 





SUBROUTINE FFT (A, B. NTOT. N. NSPAN, ISN) 

C USING Ml XED-P ADI X FA5 T FOURIER TRANSFORM AIGDPITHM. 

C MULT IVARIATE COMOUF* FOUPIER TRANSFORM, COMPUTFD IN PL AC F 

C BY R. C. SINGLETON, STANFORD RESEARCH INSTITUTE, OCT. I96B 

C ARPAYS A AND 0 ORIGINALLY HOLD THF FFAL AND IMAGINARY 

C COMPONENTS OF THE DATA, AND RETURN THE REAL AND 

C IMAGINARY COMPONENTS OF T H c RESULTING FOUR IFF) COEFFICIENTS. 

C MULTIVARIATE OATA IS INDEXED ACCORDING TO THE FORTRAN 

C ARRAY FLEMENT SUCCESSOR FUNCTION, WITHOUT LIMIT 

C ON THE NUMBER of IMPLIED MULTIPLE SUBSCRIPTS. 

C THE SUBROUTINE IS CALLED ONCE FOR EACH VAPIATF. 

C THE CALLS FOR A MULTIVARIATF TRANSFORM MAY BE IN ANY OBOFR. 

C NTOT IS THE TOTAL NUMOFF OF COMPLFX DATA VALUES. 

C N IS THE OIMENSION OF THF CURRENT VARIABLE. 

C NSRAN/N IS THE SPACING OF CONSFCUTIVF OATA VALUES 

C WHIl E INOFX1NG THF CUPPTNT VAR|ABIF. 

C THE SIGN OF ISN DFTFPMINFS THE SIGN OF THE COMPLEX 

C EXPONENTIAL. AND THF MAGNITUDE OF ISN IS NORMALLY ONE. 

DIMENSION All), 0(1) 

C A TP I - VAR I AT E TRANSFORM WITH A (N I ,N? ,N3 I . B(NI.N2.N3I 

C IS CCMPUTFD BY 

C CALL FFT { A,B,N1*N2*N3 , Nt , Nl , | ) 

C CALL FFT< A.B.NI*N2*N3.N2.NI *N2 , I ) 

C CALL FFTIA.B.Nl AN2AN3.N3.Nt *N?*N3, I ) 

C FOR A SINGLF-VARI ATF TRANSFORM, 

C, NTOT a N = NSPAN = (NUMBER OF COMPLEX DATA VALUES), E.G. 

C CALL FFTCA.B.N.N.N.t ) 

C THE DATA MAY ALTERNATIVELY OR STORED IN A SINGLF COMPLEX 

C ARPAV A, THEN THE MAGNITUDE OF ISN CHANGED TO TWO TO 

C GIVE THF CORRFCT INDEXING INCREMENT AND AI2) USED TO 

C PASS THE INITIAL ADDRESS FOR THE SEQUENCE OF IMAGINARY 

C VALUES. E.G. 

C CALL FFTI A, A( 2) .NTOT.N .NSOAN.2 ) 

c ARRAYS ATIMAXF). CKCMAX C ), BT(MAXF), SKIMAXF). AND NPIMAXP) 

C APE USED FOP TEMPOBAPY STORAGE. IF THE AVAILABLE STORAGF 

C IS INSUFFICIENT. THE PROGRAM IS TERMINATED BY A STOP. 

C MAXF MUST BF .GE. THE MAXIMUM PRIME FACTOR OF N. 

C MAXP MUST RE .GT. THE NUMDER OF PRIME FACTORS OF N. 

C IN ADD I'. ION. IF THE SOUABE-FRFE PORTION K OF N HAS TWO OR 

C MORE PRIME FACTORS. THEN MAXP MUST HE .GE • K-t . 

C AFPAV STORAGE IN NPAC FOR A MAXIMUM OF It FACTORS OF N. 

C IT N HAS MORE THAN ONF 50UABE-FREF FACTOR, THE PPOOUCT OF THE 

C SQUARE-FREE FACTOPS MUST BE .LF. 1501 

dimension nf aci?i ) ,npi isoo ) 

c APR AY STORAGE FOR MAXIMUM PRIME FACTOR OF 373 

OIMENSION AT (373) .CM 373 ).BTI 3731.SK 1373) 

EQUIVALENCE II.II) 

C THE FOLLOWING TWO CONSTANTS SHOULD AGREE WITH THF ARRAY DIMENSIONS. 

MAXF=373 
MAXRatSOO 

IF(N .LT. 2) RFTUPN 
INC* ISN 

RAD*S.O*ATAN< | .0) 

S72*PA0/5.0 

C72»C0S(S72) 

S72aSIN(S7?| 



5l20®SORT< 0. 75» 

|F» ISN .cr. 0» GO TO to 

S72--S72 

SI?0=-SI2C 

U»t> = -p AO 

INC--INC 

10 NT® t NC4NTDT 
KS=INC*NSPAN 
KSOAN^K S 
NN=NT— INC 
JC-K S/N 

BiOT=B»n* r LO»T( JCM0.fi 
1-0 
JF = 0 

C OF T TPM INF THE FACTORS O r N 

M=0 
K«N 

GO TO 20 

|5 HrM«| 

NT AC (m®4 
K*K/ 16 

20 IF! K-U/16M |6 ,ro. 01 GO TO 15 
J-3 
JJ=*o 

GO TO 30 
25 M*MF| 

NFACIMMJ 

K*K/JJ 

30 IF CHOO(K.JJ) .EO, 01 GO TO 25 

J-J42 
JJ=J»«2 

I F ( JJ .LF. K» GO TO 30 
iroc .GT. 4) GO TO 40 
KT= M 

NFAC (MMMK 

|F(K ,NE. 1) M-NFl 

GO TO no 

40 1FCK-<K/4M4 .NH. 01 GO TO 50 

M«MM 
NFAC l MM2 
K»K/4 
50 KTsrM 

J-2 

50 IF|MOr»<K.J» .NF. 0» GO TO 70 

INMtl 
NFACIMMJ 
K»K/J 

70 J«<f JMI/2M2+! 

IFU .LF. K ) GO TO 50 
M |F|KT ,E0. 0) GO TO 100 

J-KT 

«0 *•«♦! 

NFAC«M|»NFAC< J> 

J- J- I 

IM J .NF. 0) GO TO 90 
C COMPUTE FOURIER TRANSFORM 
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too 


c 

2t0 


220 

230 


C 

320 


S0=PADF /FLOAT (KSPANt 
CD=?.0*SIN< SO ) **2 
SO=SIN(SD+SD» 

KK* I 
I* I ♦ 1 

IF I NF AC (I I .NF. 21 OO TO *00 

TRANSFORM FOR F*CTOP OF 2 (INCLUDING ROTATION FACTOR! 
KSPAN=KSPAN/2 
Kl=KSPAN*2 
K2sKK*KSPAN 
AK=A<K2) 

OK=HIK2> 

AIK2»=A<KK >-AK 
R|K2»=B(KK»-BK 
AIKK ) = A< KK | ♦ AK 
IIIKK )=B|KK ) A BK 
KK-K 2AKSPAN 
IF IKK .LF. NN> GO 
KK=KK-NN 

IFIKK .LF. JCI GO 
IFIKK .GT. KSPANt 
0=1 .0-CD 
M=SD 

K2=KKtKSPAN 
AK= A ( KK)*-AlK2) 

BK=FMKK >-fMK2 > 

AIKK >=A( KK) AA(K2) 

BIKK)=H<KKJAR(K2J 
A|K2 ) = C1 *AK-S1 *RK 
R(K2t=Sl*AK*Ct*BK 
KK=K2*KSPAN 
IFIKK .LT. NT | GO 
K2=KK- NT 
Cl=-Cl 
KK>K I-K2 

IF|KK .GT. K 2 ) GO 
AK=CI-(CDAC»ASD*Sl ) 

SI = I SO * C I — CD* S t ) AS 1 

THF FOLLOWING THREE STATEMENTS COMPENSATE FOR TRUNCATION 
EPPOP. IF POUNDED ARITHMETIC IS USED. SUBSTITUTE 
CI=AK 

CI=0.5/( AK*»2«-SI**2) *0.5 
St=C l*SI 
Cl — C 1 * AK 
KK*KKAJC 

IFIKK .LT. K2t GO TO 230 
K!*K|AINCAINC 
KK=IKI -KSPANI/2* JC 
IFIKK .LE. JC+JCI GO TO 220 
GO TO 100 

TPANSFORM FOP FACTOR OF 3 COPTIONAL CODE! 

K1*KKAKSPAN 
K2=K IAKSPAN 
AK«AIKK) 

RK=tMKK> 

A J=A(K1)4AIK2> 


TO 713 

TO 210 
GO TO 800 


TO 230 
TO 230 


r j=rik 1 1 ♦n ( k 2 1 

*( KK )=»K »A J 

«|KK )=BK*RJ 

AK= -0. *54 AJ 4AK 

UK*- 0. 5*BJ4BK 

AJ*I ACKi I-AIK2 > > *SI2Q 

EJ=I RlKt J-BIK2) >4SI20 

Al K I )=AK-BJ 

BlKl ) = RK * A J 

AIK? ) = AK4BJ 

01K2 )=RK-AJ 

KK=K?4KSBAN 

IF IKK .LT. NN) GO TO 320 
KK=KK-NN 

IF(KK .LE. KSPANI CO TO 3?0 
GO TO 700 

C TRANSFORM Frjp FACTOR OF 4 

400 ir(NFACII) .NF. 4) CO TO 600 
KSPNN=KSP AN 
K5PAN=KSPAN/4 
410 0=1.0 

51=0 

4?0 KI*KK*KPPAN 
K?=K 1 4KSPAN 
K 3=K 24KS PAN 
AK *>= A(KK)4A|K?) 

AK*=A|KK J-AIK2 ) 

AJP=A(KI I 4 A { K 3 ) 

AJMcMKl ) — A ( K 3 I 
AlKK ) = AK P4 A JP 
AJPxAKP-AJP 
BKP*BIKK)4B|K2) 

DKM=0lKK>-B|K2) 

«JP=BIK1 >40<K3> 

BJM=B(KI )-R(K3) 

Bl KK |=RKP40JP 
FJD»OKP-BJP 

IFIISN .LT. 0) GO TO 450 
AK P= AKM-BJM 

akmb akm+sjm 
bkp*bkmaajm 

£KM=RK M- A JM 

IF | S I .50. 0.01 GO TO 460 
430 AlKl I * AK P*C 1 — RKP* 5 I 
Bl K1 |=AKP*SI 4BKP*CI 
AIK2 » = AJP*C2-BJP*S2 
BlK2>=AJP*S2ABJP*C2 
BlKJI=AKM*S3Tr5K"FC3 
AIK 3I*AKM* C 3—BKN* S 3 
KK=K34KSPAN 

iriKK .LE. NTI GO TO 420 
♦40 C2*CI-ICn*C|4S0*5l » 

SI*ISO*Cl-CO*SI 1451 

C THE FOLLOWING THRrE STATEMENTS COMPENSATE FOR TRUNCATION 

C EPPOP. IF POUND:!) ARITHMETIC IS USED. SUBSTITUTE 

C CI*C2 




ei=0.S/’(C2**2fSt**2>*0.5 

S1*C l*SI 

CI»CI*C2 

C2»C1**2-S1**2 

S2*2.0*C1«'SI 

C 3*C2»C1 —52* S I 

S3»C?*SI ♦S2*C I 

kk«kk-nt*jc 

IFIKK .LF. K5PANJ GO TO 420 
KK*KK-KSPAN* INC 
IF IKK .LE» JD GO TO 4 10 
IFIKSP4N .FO. JC> GO TO POO 
GO TO 100 
450 AKP=AKM»BJM 
AK**~AKM-FIJM 
BKP=BKK- AjM 
BKM*HKMf AJM 

IFISI .NE. 3. 01 GO TO 433 
460 A(KI)=AKP 
WIKI J=0KO 
A I K2 ) - A JP 
BIK2 I = 8JP 
A I K 3 ) = AK N 
BIK J)=BKM 
KK=K3*KSPAN 

IF I KK .LE. NT I GO TO 420 
GO TO 440 

C TRANSFORM FOP FACTOP OF 5 (DPT IONAL CODE I 

510 C2-C 72** 2- S 72**2 
5?=2.0*C72*R72 
520 KI=KK*KSPAN 
K2=K l+KSPAN 
K3=K2*KSPAN 
K4=K 34KSPAN 
AKP^AIKI )*AlK4 ) 

AKM*A( K i )- At K4 ) 

BKP=«BIK I )+BtK4> 

BKM=BtKl )-BlK4> 

AJP*A(K2l*At K3> 

A JMa AI K2)— AIK3I 
BJP*BIK2I*8IK3> 

BJMtxP(K2 I-0IK3) 

AA&AIKK I 
BB^BIKK) 

A I KK )<*AA *AKP* A J° 

BCKK 1*00 tBKP*BJP 
AK>AKP*C72*A JP*C2*AA 
6K»BKP*C72*B J»*C2*RR 
AJ»AKM*S72*A JM*S2 
B J*BKM*S 7p *B JM *S2 
A I K I l*AK-BJ 
AIK4 )kAK»BJ 
BIKI | *BK *A J 
6IK4I=BK-A J 
AK~AKP*C2+AJP*C72» A A 
PK>nKP*C2+BJP*C72*BP 



i 



AJ=AKM*S2-AJM*S72 

ack 2 i=»K-nj 

AIK3J-AKFBJ 
BCK2 J=BK*AJ 
»<K3J=HK-AJ 
KK-K **KSPAN 

IF<KK .LT. NN) GO TO 520 
KK-rKK-NN 

IFCKK «LE » K5PAN» 50 TO 520 
OP TO 703 

C TPAN5FORM FOP OOD FACTORS 

600 K = NFACII) 

K«PNN=KSP»N 

KSPAN=K SPAN/K 

|F(K .EO. 3) GO TO 323 

IFCK .EO. 5> GO TO 510 

IFC* .'O. JFJ GO TO 6*0 

JF=K 

51=9 AO /FLOAT C K > 

CI=COS(Sl ) 

5I = S!N< SI > 

IFJjr .GT. MAXFI GO TO 996 
CM JF|»J .0 
5K( JF)*0.0 
J- I 

A30 CMJ>*CMK)*CI*5MIC>*SI 
SM J»»CMM*S|-SMkI*CI 
K*K-I 

CK(ttl*CK(JI 
SM* > = -SK( J> 

J*J*1 

IF t J .LT. K ) GO TO 630 
640 K | *KK 

K2*KK*KSPNN 
AA*A IKK J 
BB=B<KK> 

AK*AA 

W=BO 

J*t 

KI«X I4KSPAN 
650 K2*K2-KSPAN 
J*JT1 

ATI J)*A(KI >AA(K?) 

AK= AT ( J ) ♦ AK 

BTC J»=B(K | l+BCK?) 

BK*BTl J)*BK 

AT( JJ=A(K! J-A‘K2! 

BTC J»*6<KI»-B<K2> 
«l>Kt«KSPAN 

irtKI .LT. K2I GO TO 650 

ACKK>>AK 

BCKKl-BK 

m«kk 

K2-KK*KSPNN 
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J )T | 

66 0 Kt*KI«5P»N 

KZ=K2-KSPAN 
JJ=J 
AK*AA 
BK=BR 
AJ=0.0 
PJ=0 .0 
K*l 

670 K*K»I 

AK*AT( K1 *CK( JJ» *-AK 
AK*RT(K) *CM JJI FBK 
K*Kf 1 

AJ=AT(K)*SK( J J ) ♦ A J 
RJ=8TC K) *RKC JJ> »BJ 
J 

IF I J J .GT. JFI JJ-JJ-JF 

IF |K .LT. JF| GO TO 670 

K*= JF- J 

AIKI )=AK-RJ 

OIKI 1 = RK ♦ A J 

ACK?)=AK*BJ 

RIF? >=BK-A J 

JxJS'l 

IF tJ .LT. Kl GO TO 660 
IKK=KKFKSPNN 

IFCKK iLt. NNl GO TO 6*0 

ftf if x K K -o IJ FJ 

IFCKK .LE. KSrt'ANI GO TO 6*0 

C MULTIPLY BY POTATION FACTOP CEKCE®>T FOP FACTORS OF 2 AND *1 

700 IF C I .EO* Ml GO TO 600 
KK*JCf I 

710 C?= 1 .0“CD 

SI = S0 
720 C1=C2 

«j=5 | 

KK*KKFKSPAN 
730 AKsA(KK) 

AIKK l = C2*AK-S2*tl(KK ) 

HCKK) = S?*AK4-C2*BCKK| 

*K*KK+KSFNN 

IFCKK .LF. NT) GO TO 730 

AK»SI*S2 

S2"51*C?*CI *S2 

C2*C l*C2-AK 

KK«KK-NT*KSPAN 

IFfKK .LE. KSPNN) GO TO 730 

CZ>«Ct-ICO*CI FSD*SI ) 

SI*SI*CSO*CI-CO*SI ) 

C THE FOLLOWING THREE STATEMENTS COMPENSATE FOP TRUNCATION 

C FFPOR. IF FOUNDED ARITHMETIC IS USED. THEY MAY 

C BE DELETED. 

Cl »0.5/C C2**2*S 1 **2) *0.5 

SI«CI*SI 

C2«CI*C2 

KK»KK-KSPNN* JC 
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IF < KK «UF. KS°AN } GO TO V20 
KK*KK~KSP AN4 jc + INC 
IF(KK .LE. JC4JCI GO TO 710 
GO TO 100 

C PEFMUTE THE P^SULTS TC NORMAL OPOER- — DONE IN TWO STAGES 

c PFPMUT AT ION FOR SOUAPF FACTORS OF N 

800 ND(1 l=KS 

IFCKT ,EQ. 0) GO TO 890 
K*KT AKT+I 

IF(M »LT. K I K,=K-t 

J* 1 

NP(K*T > = JC 

•10 NPI JAI l*NP< JI/NFACC Jl 
NP(K t*NP(K»l >*NF ACC J > 

J«JM 

K«K-1 

I f c j .lt. k> go to eto 

K3*NO(KA 1 ) 

KSPAN*NP(2> 

KKwjCM 
K2-KSPAN* | 

J=£ 

IF ( N .NE. NTOTI GO TO 850 

C PERMUTATION FOP S INGLF-VAR I ATE TRANSFORM (OPTIONAL CODE I 

•20 AK= A (KK ) 

A(KKl=A(K2t 

4 ( K2 I = AK l 

BK»B(KK> 

R(KK)*R<K2> 

PCK2)*BK 

KKmKKMNC 

«2»KSPftNAK2 


830 

IFOC2 .LT. KS> 
K2*K2-NP( J) 

J= J«- 1 

K2=NP( J*1 )*K2 

GO 

TO 

829 


IF(K2 .GT. NP ( J ) ) 
J*l 

GO 

TO 839 

840 

IF(KK .LT. K2) 

KK»KKHNC 

K2>KSPANPK2 

GO 

TO 

820 


IFCK2 -LT. MS) 

GO 

ro 

840 


IF C KK .LT. KS) 
JC*K 3 

GO 

TO 

8 30 


GO TO 800 

C PEPMUTATION FOR MULT tVARI ATE TRANSFORM 

•SO KaKKPJC 
060 AK*A(KK> 

A(KK)>A(K2) 

A(K2)*AK 

BK-B'KK) 

0(KK |*B( K2) 
eiK2l*RK 
KK*KK+ INC 
K2»K2*INC 

IF(KK .LT. KI GO TO 860 
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KK-KK^KS-JC 

K2=K2»KE-JC 

IF IKK .LT. NT) GO TO 850 
K 2=«K 2- NT ♦KSf’AN 
KK*KK-NT*JC 

IF(K2 .LT. KS) GO TO 850 

ero k2*k2-np(j: 

J*J*1 

K2=NP( JFl ) *K2 

IF IK 2 .GT. NP(J)> GO TO 870 
J*=I 

MAO IFIKK .LT. K?) GO TO 850 

kk=kk*jc 

K2*KSPAN»-K2 

IF I K 2 .LT. KS) GO TO PPO 
IFIKK .LT. KS) GO TO 870 
JC*K3 

890 !F(?*KTM .GF. M) RETURN 
KSPNN=NPIKTT I ) 

C PFRMUT AT ION FOR SOUARF-FRFE FACTORS OF N 

JxM-KT 
NFACl Jft )=l 

900 NF AC I J ) = NFACI J)*NFACI JFl ) 

J-I 

IF! J .NE« KT) GO TO 900 

KT»KT ♦ 1 

NN*NFAC(KT)-i 

IF I NN .GT. MAXPI GO TO 998 

JJ=0 

J*0 

GO TO 906 
902 JJxJJ-K2 
l',2=KK 
l«*KM 
KK=NCAC(K) 

90* JJxKKFJJ 

IF I J J .G€. K 2 ) GO TO 902 
NPI J)=JJ 

906 K2XNFACCKT) 

K=KT ♦ I 
KKxNFAC(K) 

JxJM 

IF|J .LE. NN) GO TO 90* 

C DETERMINE THE PERMUTATION CVCLES OF LENGTH GREATER THAN t 

J*0 

GfJ TO 91* 

910 K = KK 

KKsNP(K) 

NPt K )»-KK 


IFIKK «NE * 
K3«KK 

J) 

GO 

TO 

910 

J*Jf I 
KK*NPt J) 
IFIKK .LT. 

0) 

GO 

TO 

91* 

>fjkk .NE. 

J ) 

GO 

TO 

910 


NP( J )■=- J 


|F|J .NE. NN1 CO TO 914 
NAXF=INC**AXF 

C REORDER * AND B. FOLLOWING THE PERMUTATION CYCLES 

CD TO «*50 
•2« J-J-l 

IF C NPI J ) «LT • 0) CO TO o?A 
jj~jr 

9 2ft KSPAN=JJ 

IFIJJ .r.T. MAXF| KSPANbMAXF 

JJ»JJ-KSPAN 

K-NPCJI 

KK*JC*K»U v JJ 

Kt*KK*K5PAN 

K2*0 

92* K2»K2*I 

4TCK2)*ACKI I 
BTC K 2 J *B CKI » 

KI*K l-INC 

IFIKI .Nf, KKl GO TO °2* 

0.1? F |*KK>KFPAN 

K2«K l-JC*CK*NPfK> | 

K*-NP(KI 

036 AC K I I = A C K? ) 

B(K| )*BCK2> 

K|*KI-INC 

F2*K2-INC 

|F<«| .NF. KK) GO TO «>36 
KK*K2 

IF C K .NE. Jl CO TO 932 
K|J-KK«-KSPAN 
K2>0 

040 K2=K 24 1 

At K I >*AT(K2> 

B(K|)»B7(K2) 

KI-K1-INC 

IFIKI .NE. KK I GO TO 940 
IFIJJ . Ntl r . 0) GO TO 926 
IF < J .NE. I) GO TO 924 
930 J~K3*8 

NT*NT-KSPNN 
I I*NT— I NC* I 

IF | NT .GE. 3) GO TO 924 
RFTIJRN 

C EPRPR FINISH. INSUFFICIENT ARRAY STORAGE 

996 ISN»0 
C PRINT 990 

999 FORMAT U4MC.*,fKV*.Y CSCCffOE© !sTTHIN SUBROUTINE FFT I 

RETURN 
END 


SURPOUTINF PEALTP(A«B«K»1 RN I 
DIMENSION A I I ) • R< I ) 

PEAL IM 
INC=IAHS(ISN> 

MK*N»!NC+? 

NH>NK / ;> 

SDs 2 • O * AT AN ( I .0)/FLOAT(N| 

CO=2.0*S!N< SO ) **2 
Sr»-SIN( SDTSD) 

SN=0,0 
CN=1 .0 
AINK-I )=A( I > 

RINK-l )=BI 1 ) 

13 OO 70 J — I • NH • I NC 
K=NK-J 
AA = A( J ) ♦ A( K ) 

AH* A I J I — A< K I 
RA=B( J »*R<K ) 

FB=RI J» -R(K» 

PE=CN*BA»SN* AM 
!M=SN*RA-CN* AO 
PI*>=IM-RR 
Bt J1=I M*RR 
A<K J = AA-RE 
A< JI^AA+PE 
AA=CN- tCr>*CN*SO*SN» 

SN* { SD*CN-CD*SN I FRN 

C THE FOLLOWING THREE STATEMENTS COMPENSATE FOR TRUNCATION 

CN*0.*»rl AA**2*SN*A2M0.S 
SN=CN*SN 
70 CN=CN*AA 
RETURN 
END 
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